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Finite  element  analysis  of  linear  structures  is  an  increasingly  popular  tool  of 

choice  for  analysis  by  engineers.  Only  recently  has  this  tool  been  applied  to  the 

case  of  nonlinear  cable  structures  with  a  sufficient  degree  of  veritable  accuracy. 

Work  in  this  dissertation  provides  a  finite  element  approach  to  the  analysis  of 

cable  structures  with  verification  of  the  approach  done  for  a  traffic  signal  cable 

structure.  There  has  been  a  preponderance  of  work  done  on  the  dynamic  analysis 

of  single  cables  and  cable  networks  but  very  little  is  known  about  dynamic 

analysis  of  non-grid-like  cable  elements  in  conjunction  with  other  elements  such  as 

beams  and  trusses.  The  approach  of  analysis  presented  in  this  dissertation 

xi 


provides  a  valuable  bridge  in  available  theories  by  creating  a  program  that 
dynamically  analyzes  a  cable  in  context  of  other  types  of  elements  also  present 
within  the  structure.  Using  the  principle  of  virtual  work,  the  cable  element 
developed  in  this  dissertation  is  based  upon  a  total  Lagrangian  formulation.  The 
cable  element  models  geometric  nonlinear  capability.  Also  modeled  in  the 
program  developed  is  a  nonlinear  loading  capability.  A  realistic  wind  load 
generation  was  accomplished  by  the  use  of  auto-regressive  numerical  filters.  With 
a  complete  finite  element  model  sufficiently  exhibiting  the  elastic  and  kinematic 
properties  of  the  true  structure,  dynamic  analysis  was  accomplished  by  numerical 
integration  of  the  equation  of  motion.  An  incremental  form  of  the  general  equation 
of  motion  is  used  because  of  a  piecewise  linear  assumption.  The  numerical 
integration  scheme  presented  in  this  dissertation  is  the  Newmark's  Beta  method. 
The  Newmark's  Beta  method  reduces  the  equation  of  motion  to  a  set  of  nonlinear 
algebraic  equations.  These  equations  are  then  solved  using  a  modified  Newton- 
Raphson.  The  finite  element  approach  presented  in  this  dissertation  is 
programmed  into  a  software,  DYNASS.  Verification  of  the  program  DYNASS 
and  subsequently  the  solution  approach  is  done  by  comparing  predictions  from 
DYNASS  to  test  results  obtained  from  recent  traffic  signal  experiments  performed 
at  the  Civil  Engineering  Department,  University  of  Florida. 
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CHAPTER  1 
INTRODUCTION 


Problem  Definition 


The  use  of  cables  as  load-carrying  elements  is  common  for  structures  such 
as  bridges,  cable  reinforced  membrane  roofs,  traffic  signal  systems,  spatial  cable 
nets  and  mooring  lines.  The  set  of  forces  seen  in  each  cable  application  is  different 
but  there  is  a  common  trend  in  the  behavior  of  these  structural  systems.  This 
behavior  is  nonlinear.  One  source  of  nonlinearity  is  the  large  displacement  the 
cables  undergo.  This  type  of  nonlinear  behavior  is  known  as  geometric 
nonlinearity.  Several  theories  have  been  proposed  to  handle  this  problem.  One 
such  theory  is  presented  in  this  report. 

The  Structure 

The  structure  investigated  here  is  the  traffic  signal  system.  The  choice  of 
this  type  of  cable  stayed  structure  was  a  result  of  the  recent  work  done  for  the 
Florida  Department  of  Transportation.  The  Civil  Engineering  Department  at  the 


University  of  Florida  developed  the  software  program,  ATLAS,  to  perform  a  static 
analysis  of  a  two-point  connection  traffic  signal  system  as  shown  in  Fig.  1.1.  The 
restrictions  encountered  during  developing  this  program  spurred  the  need  for 
another  analysis  approach.  The  typical  system  considered  consists  of  a  catenary 
cable  (slack  cable  configuration),  a  messenger  cable  (taut  cable  configuration),  the 
signal  heads  and  signs  (points  of  dead  and  wind  load  application),  and  the 
prestressed  concrete  poles  (as  the  support  system). 


Catenary  Cable  Prestressed  concrete  pole 

/  Connection  hanger 


Fig.  1.1    Traffic  Signal  System 


The  traffic  signal  system  considered  includes  that  shown  in  Fig.  1 . 1  but  is 
not  limited  to  the  two-point  connection  traffic  signal  system  shown.  The  method 
presented  in  this  research  is  capable  of  handling  virtually  any  system  of  cables  and 
frame  elements,  for  example,  the  one-point  connection  traffic  signal  system  and 
the  suspended  box  system. 

The  Theory 

The  approach  presented  here  uses  the  theories  of  structural  dynamics  to 
model  the  behavior  of  cables  subjected  to  both  dynamic  wind  loads  and  static 
loads.  The  static  load  case  can  be  considered  a  special  case  of  the  general  dynamic 
model.  The  static  load  is  modeled  as  a  constant  load  with  duration  equal  to  the  full 
time  of  analysis.  The  structural  response  often  of  interest  in  the  static  load  case  is 
the  damped  out  response,  i.e.  the  stresses,  moments,  forces  and  deflections  at  die- 
down.  For  this  response  the  initial  rate  of  application  of  load  is  inconsequential.  It 
is,  however,  good  practice  to  maintain  a  slow  ascent  to  full  load  in  order  to 
minimize  impact  effects  of  the  load  and  to  obtain  a  dynamic  response  more 
reflective  of  a  stationary  load. 
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Fig.  1.2  -  Static  load  model 


The  wind  load  is  modeled  as  both  static  and  dynamic.  The  equivalent  static 
wind  load  as  determined  by  AASHTO  (American  Association  of  State  Highway 
Officials,  "Standard  Specifications  for  Structural  Supports  for  Highway  Signs, 
Luminaries  and  Traffic  signals")  is  applied  as  a  horizontal  force  (called  the  drag 
force)  and  a  vertical  force  (called  the  lift  force).  An  aspect  of  nonlinear  behavior 
of  traffic  signal  cable  is  the  nonlinear  wind  loading.  The  research  work  presented 


provides  a  dynamic  wind  model  comprising  a  steady  wind  pressure  and  a  variable 
wind  component.  The  variable  wind  component  models  wind  gust  effects.  A 
wind  gust  is  a  short  duration  wind  turbulence  often  characterized  by  its  frequency 
content  and  its  duration.  Fig.  1.3  shows  a  typical  wind  gust  representation. 


Fig.  1.3  -  Mean  and  time  varying  components  of  wind  speed 


The  dynamic  theory  presented  uses  geometric  and  time  discretization.  The 
geometric  discretization  or  discrete  element  technique  employed  is  the  finite 
element  method.  Numerical  integration  is  used  in  the  solution  of  the  equations  of 
motion  in  the  discrete  time  domain.  Generally  a  nonlinear  response  analysis 
involves  a  significant  amount  of  computation  as  it  was  in  this  case.  Consequently, 
a  computer  program  was  written. 

The  Computer  Program 

The  resulting  program  is  called  DYNASS.  This  program  performs  a 
nonlinear  dynamic  analysis  of  the  traffic  signal  cable  system.  This  program  is 
written  in  Fortran  with  some  wind  generation  routines  in  C.  The  results  of  this 
program  are  compared  with  results  from  other  programs  where  possible  and  with 
test  results  from  full  size  traffic  signal  system  tests. 


Overview 


Each  chapter  in  this  report  discusses  topics  or  themes  built  upon  in 
subsequent  chapters.  Topics  discussed  start  from  the  theory  proposed  to  the 
program  written,  to  the  comparison  of  the  results  from  the  program  with  test 


results  and  finally  to  the  typical  applications  and  conclusion.  Figures  and  formulas 
are  numbered  consecutively  throughout  each  chapter. 

Content 

Chapter  two  presents  a  literature  background.  This  includes  an  overview  of 
existing  work  done  on  this  topic  and  an  explanation  of  the  fundamentals  of 
structural  dynamic  theories  and  cable  analysis  theories. 

Chapter  three  contains  a  detailed  presentation  of  the  theories  and  methods 
developed  in  this  research.  These  include  the  various  assumptions  made  and  the 
equations  derived.  This  chapter  is  essentially  subdivided  into  two  parts.  The  first 
part  is  the  theoretical  development  just  mentioned,  and  the  second  part  discusses 
the  algorithm  used  in  the  program,  DYNASS. 

Chapters  four  discusses  the  wind  load  models  developed  plus  a  general 
description  of  the  other  loads  on  the  traffic  signal  system.  There  are  two  wind 
models  formulated,  the  variable  area  wind  load  and  the  variable  velocity  wind 
load.  Both  models  are  discussed  in  some  detail  in  chapter  four.  Finally,  chapter 
four  presents  a  brief  discussion  of  the  wind  load  data  obtained  from  dynamic 
experiments  performed.  These  experiments  were  conducted  previously  for  a 
Florida  Department  of  Transportation  sponsored  research. 


Chapter  five  and  six  present  comparisons  of  the  results  from  DYNASS  to 
other  software,  where  available,  and  to  test  data  obtained  from  experiments  done. 
Chapter  five  compares  results  for  static  load  cases  only  and  chapter  six  deals  with 
the  comparison  of  results  for  dynamic  wind  load  cases.  Each  chapter  is  concluded 
with  a  discussion  of  strengths  and  limitations  of  either  the  theory  proposed  and/or 
the  program  developed,  for  each  load  case. 

The  final  chapter,  chapter  seven,  deals  with  conclusions  observed  from  the 
research  and  recommendations.  Areas  needing  further  research  are  also  identified 
and  noted. 

Also  included  are  Appendices  A  and  B.  Appendix  A  shows  a  detailed 
derivation  of  the  2-node  cable  element  developed  in  this  research.  Appendix  B 
gives  a  user's  manual  of  the  computer  program  developed.  In  appendix  B,  defaults 
assumed  in  the  program  are  explained  and  for  portions  of  the  program  not 
developed  by  the  author,  references  are  provided. 

Scope 

The  scope  of  the  research  presented  is  limited  by  the  inherent  assumptions 
made  in  the  cable  analysis.  The  assumptions  made,  common  to  other  types  of 
cable  analysis,  are  as  follows: 
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♦  The  cables  are  perfectly  flexible  and  possess  no  bending  stiffness,  only 

axial  stiffness. 

♦  Cable  materials  are  of  hookean  materials.  Ferrous  materials  are  commonly 
used  for  cables. 

♦  Only  small  strains  are  involved,  i.e.  no  yielding  of  cables  ever  reached. 
By  the  nature  of  the  type  of  structure  investigated  (traffic  signal  systems), 

the  loads  are  essentially  predefined.  For  this  research,  loads  examined  are 
restricted  to  the  dead  load  from  the  weight  of  the  signal  heads  and  signs,  and  wind 
loads.  Other  assumptions  discussed  in  chapter  four  further  define  the  application 
of  the  theory  and  subsequently  the  program. 


CHAPTER  2 
LITERATURE  BACKGROUND 


Introduction 

Several  authors  have  dealt  with  the  problem  of  the  geometric  nonhnearity 
of  cables  and  have  provided  some  solutions  to  the  problem.  A  brief  summary  of 
pertinent  results  of  their  work  is  given  in  this  chapter  as  they  provide  a  foundation 
to  the  author's  work.  Also  presented  is  a  general  theoretical  background  of  cable 
analysis  and  dynamic  analysis.  The  work  covered  here  reflects  only  the  solutions 
provided  for  cable  dynamics.  It  is,  however,  congruent  to  consider  the  static  cable 
problem  as  a  special  case  of  the  cable  dynamics  problem  hence  extending  these 
solutions  to  cover  static  cable  analysis.  Other  methods  strictly  for  static  solutions 
include  the  force  density  method  [1]. 

Fundamentals  of  Cable  Analysis 

Consider  a  uniform  cable  which  hangs  in  static  equilibrium  in  a  vertical 
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plane  suspended  between  two  supports.  The  forces  acting  on  a  differential 
element  of  the  cable  is  dependent  of  the  profile  assumed  by  the  cable  at-rest 
position  and  the  cable  final  deformed  position.  Let  any  point  on  this  cable  be 
described  in  the  3-D  Cartesian  coordinates  shown  in  Fig.  2.1. 


(u) 


Fig.  2.1  Cartesian  Coordinate 


Let  the  cable  ordinate  after  deformation  be  denoted  x,  y,  z  as  shown  in  Fig.  2.2. 
x  =  x  +  u 

y  =  y  +  v 
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z  =  z  +  w 


The  cable  segment  shown  in  Fig.  2.2  below  can  be  described  by  a  vector  along  the 
chord  joining  both  ends  of  the  cable  segment. 


-xLi  +  IjJllLj  +  llIjLk 


ds 


ds 


ds 


The  terms  *j  -  *,-      y}  -  >,    and    *j  -  z,-      are  referred  to  as  the  direction  cosines  for 

ds  ds  ds 


the  element. 


wl£ 


(x,>  y» z.) 


Fig.  2.2  Displacement  of  a  cable  segment 


The  relationship  between  the  nodal  forces  and  displacements  at  any  point  along 
the  cable  relative  to  the  next  is  defined  as  follows, 
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Uj^Ui+J7  (2-1) 


This  says  that  the  displacement  at  nodey  is  the  sum  of  the  displacement  at  node  i 
and  the  differential  strain  in  the  cable  multiplied  by  the  arc  length  of  cable. 
Similarly  for  the  displacements  in  the  y  and  z  direction,  the  relationship  is 

vJ  =  v'  +  T7ds  (2-2) 


d  w   , 
wj  =  wi+—ds  (2-3) 


The  above  displacement-strain  relationship  is  referred  to  as  the  linear  strain  - 
displacement  relationship.  The  relationship  between  the  forces  at  node  i  andy 
both  along  the  cable  can  similarly  be  summarized  as  follows, 
Tj  =  Tt+dT  (2-4a) 

Fxj  =  Fxi+dFx  (2-4b) 

Fyj      =     Fyi      +     dFy  (2-4C) 

Fzj=Fzi+dFz  (2"4d) 

For  static  analysis  the  forces  acting  on  a  finite  length  of  a  cable  are 
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T         =         Cable  tension 

Fx,  Fy,  Fz=      equivalent  nodal  loads  due  to  for  example  dead  weight  of  cable 
and/or  cladding  or  applied  nodal  load.  In  the  "unloaded 
configuration"  otherwise  referred  to  as  the  null  state,  the  applied 
load  is  taken  as  zero. 

The  tension  in  the  undeformed  slack  cable  is  computed  from  the  cable  profile. 

The  basic  shape  of  a  slack  cable  subjected  to  its  own  weight  is  a  catenary  profile. 

However,  if  the  slope  of  the  cable  is  everywhere  small,  the  profile  adopted  by  the 

cable  can  be  accurately  described  by  a  parabola.  It  is  advantageous  as  a  means  of 

simplification  to  describe  cable  shape  in  two  dimensions.  Here  the  x-z  plane  is 

used. 


T: 


Fig.  2.3     Shape  of  cable  due  to  self  wt. 
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According  to  Ref.  [2]  the  parabolic  profile  is  given  by  Eqn.  2-5. 


4hx 
z(x)  =  x  tan(0) T~(L~  x) 

L2 


(2-5) 


The  tension  in  the  cable  at  any  point  x  along  the  span  of  the  cable  is  as  given  in 
Eqn.  2-6. 


T(x)  = 


w  L 


8/; 


1  + 


\2 

,A.     Ah     8hx 

tan(0) +  ^— 

L        ,2 


(2-6) 


S         =         length  of  the  cable  along  the  sag  (Fig  2.3). 

z  (x)    =         vertical  distance  of  any  point  along  the  cable  relative  the  starting 

node, 
tan  (0)  =        slope  of  the  chord  joining  both  ends  of  the  cable, 
h  =        distance  from  the  chord  at  mid-span  of  the  cable.  This  is  used  to 

define  cable  sag. 


w 
L 


weight  of  the  cable  per  horizontal  distance, 
span  of  the  cable. 


T  (x)   =         tension  in  the  cable  (null  state)  at  any  point  along  the  cable. 


sag 


Comparing  the  above  profile  with  the  catenary  profile  in  Ref.  [3],  the  following 
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equations  form  the  catenary  profile  equivalent  of  Eqns  2-5  and  2-6  above.  Note 
that  the  catenary  shape  is  more  difficult  to  define  and  indeed  is  only  defined  by  an 
iterative  process  of  trial  and  error. 


Z(X)  = 


H 


w 


cosh  ( e )  -  cosh(  2(f) e) 


(2-7) 


£  =  sinh 


-i 


<p  tan  (0) 
sinh  (<j> ) 


+  0 


(2-8) 


4>  = 


w  L 
2~H 


(2-9) 


T(x)=H     1  +  sinh 


sinh 


-i  2 


wL      wx 


2H       H 


(2-10) 


H 


horizontal  component  of  the  tension  in  the  cable.  This  is  invariant 
along  the  cable. 


The  other  variables  retain  their  definitions  given  on  previous  page. 

There  are  two  unknowns  in  the  above  equations,  z(x)  and  H.  The  value  of  H  is 
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often  obtained  by  approximately  setting  z(x)  to  sag  x  span  (L)  at  x  =  L/2.  The  rest 

of  the  configuration  is  obtained  by  using  the  H  computed.  It  is  apparent  that  it  is 

more  convenient  to  use  the  parabolic  approximation.  The  error  introduced  by  the 

approximation  is  of  the  order  of  1%  -  2%  for  slack  cables  with  small  sag.  For 

example,  if  the  tension  at  the  first  third  point  and  the  z  coordinate  were  required  of 

the  cable  with  the  following  properties: 

sag  =  10%,  w  =  0.12  plf,  L  =  3500  ft,  take  tan(0)  as  zero  i.e.  both  cable  ends  are 

on  the  same  level.  The  results  will  be  the  following. 

Assuming  the  cable  shape  is  parabolic  - 

z  (x  =  VHL)  =  3 1 1 . 1 1 1  ft.  And  T  (x  =  VaL)  x  529.646  lbf . 

Using  the  catenary  equations  above  the  results  - 

z  (x  =  V&L)  =  31 1.557  ft.  And  T  (x  =  VaL)  =  536.839  lbf. 

An  error  of  just  1.4%.  For  the  remainder  of  this  dissertation,  the  parabolic  (null) 

profile  will  serve  as  the  reference  state  for  the  slack  catenary  cable.    The  null  state 

shape  is  important  in  the  derivation  of  the  stiffness  matrix  of  the  cables  as  well  as 

in  determining  the  initial  tensions  in  the  cable.  Any  assumption  about  this  initial 

configuration  could  significantly  influence  analysis  results. 

Cables,  like  trusses,  are  idealized  as  not  having  any  bending  stiffness 

( in  reality  they  possess  some  bending  stiffness).  This  introduces  equation  of 

conditions  in  the  equilibrium  analysis  of  cables  making  them  inherently  unstable. 
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In  addition,  they  can  undergo  large  displacements  and  cannot  resist  compression. 

These  complexities  invalidate  the  more  conventional  linear  theories  often  based 

on  the  assumption  of  small  displacements  and  unconditional  stability.  To  solve 

the  problem  of  inherent  cable  instability  it  is  necessary  to  introduce  additional 

resisting  forces  that  can  keep  the  cable  in  equilibrium.  The  problem  of  large 

displacement  is  best  dealt  with  by  applying  the  load  in  small  enough  increments  as 

to  achieve  small  enough  displacements  in  each  increment.  Lastly,  the  inability  of 

cables  to  sustain  compression  can  be  taken  into  account  by  solving  the  problem  in 

steps  and  keeping  track  of  the  stresses  in  the  cable.  The  cable  elements  are 

rendered  inactive  when  in  compression  and  active  when  in  tension.  A  very  logical 

method  of  solving  a  cable  problem,  therefore,  is  to  use  the  theories  of  structural 

dynamics.  It  is  advantageous  to  use  this  technique  because  it  introduces  inertia 

and  damping  forces  that  can  be  used  to  stabilize  the  cable  and  it  provides  a  time 

step  solution  which  allows  the  application  of  load  in  small  increments  and 

supplies  a  framework  for  keeping  track  of  cable  stresses. 

Cable  Dynamics 

The  solutions  to  the  cable  dynamics  problem  can  be  discussed  in  two  major 
categories,  closed  form  solutions  otherwise  referred  to  as  the  "exact"  solution,  and 
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the  finite  element  solutions  i.e.  explicit  solutions  such  as  the  Euler's  method, 

Central  Difference  method  and  Runge-Kutta  method  and  the  more  common 

implicit  methods  such  as  the  Wilson  Theta  method  and  the  Newmark's  Beta 

method. 


Closed  Form  Solutions 

Most  closed  form  solutions,  for  taut  cables  and  cables  with  small  sag  to 
span  ratios,  are  based  on  a  linear  theory  (i.e.  small  displacement  theory).  Irvine 
and  Caughey  [4],  Uwe  Starossek  [5]  and  Veletsos  and  Darbre  [6]  have  all 
presented  semi-empirical  solutions  based  on  this  linear  theory.  Of  importance  is 
the  experimental  results  of  Irvine  and  Caughey  [4].  Their  work  revealed  the 
observed  dynamic  behavior  of  cables  with  small  sag  subjected  to  free  vibration. 
When  a  cable  vibrates  it  does  so  in  two  major  motions.  The  transverse  horizontal 
motion  and  the  in-plane  motion.  Both  motions  are  observed  to  be  uncoupled.  The 
transverse  motion  is  essentially  a  swinging  motion  as  the  cable  rocks  in  and  out  of 
its  at-rest  plane.  The  in-plane  motion  consists  of  two  components,  the  symmetric 
and  anti-symmetric  components.  The  anti-symmetric  components  consist  of  in- 
line motion  similar  to  a  wave  traveling  along  the  length  of  the  cable  and  in-plane 
transverse  motions  which  are  unsymmetrical  in-plane  deflections  of  the  cable. 
The  symmetrical  components  are  analogous  to  simple  harmonic  motions  occurring 
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along  the  cable.  For  small  sag  to  span  ratios  (<  1:8)  these  components  of  motion 

constitute  the  modal  shapes  that  could  then  be  analyzed.  The  results  of  most  of 

the  aforementioned  authors  establish  the  mode  shapes  of  cable  vibration  and  their 

corresponding  frequencies.  The  limitation  of  the  linear  theory  is  that  for  catenary 

cables  with  large  sag  the  theory  ceases  to  be  accurate    The  theory  also  becomes 

problematic  when  the  cable  is  considered  as  part  of  a  structure  and  not  standing 

alone. 

The  equation  of  motion  from  linear  theory  is  given  in  Eqn.  2-12. 

From  Uwe  Starossek  [6], 


H 


m    g   I2  (2-1D 

8    / 


H        =         Horizontal  component  of  the  cable  tension 

m        =         mass  of  the  cable  per  unit  length 

g         =         acceleration  due  to  gravity 

/  =         cable  span 

f  =         sag  of  the  cable 

The  in-plane  equation  of  motion  is  for  sag  ^  1/20; 

„  d2w  d2z  d2w         dw  (2-12) 

ti  y  +  nr =  m —  +  c 

dx  ax  dt  dt 
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T  +  dT 


Fig.  2.4  Notations  For  Eqn.  2-12  for  a 
cable  segment 


c  =         damping  force  per  unit  length 

u,  w     =         displacement  in  the  x,  z  directions  respectively. 


K  =  T 


dx 
ds 


dynamic  part  of  the  total  cable  tension.  This  is  invariant  along  the 
cable. 


An  alternative  to  the  exact  solution  is  the  discrete  idealization  of  the 
continuum.  There  are  two  major  techniques  of  discrete  idealization  of  a  structural 
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continuum,  the  lumped  mass  method  and  the  finite  element  method.  The  finite 

element  method  is  more  commonly  used  because,  among  other  reasons,  the 

numerical  modeling  has  a  clearer  physical  interpretation.  With  the  finite  method, 

the  dynamic  solution  for  cable  vibrations  are  either  explicit  solutions  or  implicit 

solutions. 


Finite  Element  Method 


Due  to  the  nonlinear  behavior  of  cable  systems,  direct  superposition  of 
loads  and  displacements  as  would  be  necessary  with  the  typical  modal 
superposition  analysis  is  not  strictly  valid.  However,  according  to  Leonard  [7] 
and  others,  if  an  incremental  load  approach  is  adopted  and  the  load  increment  is 
small  enough,  then  within  each  step  it  is  reasonable  to  assume  linearity  of  the 
superposed  dynamic  response.  An  additional  condition  for  a  linearized  solution  is 
that  the  prestressed  cable  configuration  prior  to  load  application  be  known  and 
that  the  cable  system  possess  a  significant  stiffness.  Using  finite  element  analysis 
the  traffic  signal  system  is  idealized  as  follows;  the  strain  poles  as  beam  elements, 
the  cables  as  nonlinear  cable  elements  and  the  connectors  as  large  displacement 
beam  elements.  A  common  practice  is  the  use  of  steel  bracing  for  the  poles. 
These  braces  could  be  modeled  as  linear  truss  elements. 


23 


Formulation  of  Element  Properties 

The  properties  of  the  linear  beam  elements  are  given  below,  their 
derivations  can  be  found  in  most  structural  analysis  books.  The  typical 
geometrically  nonlinear  truss  element  often  used  for  cable  analysis  is  also  given  in 
this  chapter.  The  cable  element  adopted  for  this  dissertation,  however,  is  derived 
in  chapter  3. 

Beam  Properties 

According  to  Ref.  [8]  the  stiffness  matrix  for  the  3-D  beam  element  is  given  in 
Eqn.  2-13. 
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Truss  Element  Properties 


The  degree  of  freedom  for  a  typical  truss  element  is  shown  in  Fig.  2.6. 


Fig.  2.6  Degrees  of  freedom  for  a  truss  element 


The  stiffness  matrix  for  the  truss  is  given  in  Eqn  2-15  [10].  This  stiffness 
matrix  is  in  the  element  local  coordinate. 


[K]    = 


1 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

EAX 

0 
-1 

0 
0 

0 

0 

0 

1 

0 

0 

0 

I 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(2-15) 


In  the  structure's  global  coordinate,  the  stiffness  matrix  is  given  in  Eqn.  2-16. 
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Fig.  2.7  -  Truss  element  showing  direction  cosines 
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a 


Ax 


L         C7 


A  z 


(2-17) 


V(Ax)2    +   (Ay)2    +   (AZ)2 


(2-18) 


The  mass  matrix  in  global  coordinates  is  also  given  below  [9]. 


[M]    = 


2 

0 

0 

-1 

0 

0 

0 

2 

0 

0 

-1 

0 

ml 
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0 

2 

0 

0 

-1 

6 

-1 

0 

0 

2 

0 

0 

0 

-1 
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0 

2 
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0 

0 

-1 

0 

0 

2 

(2-19) 


m 


mass  per  unit  length  for  the  element 
length  of  the  element 


Cable  Element  Properties 

A  straight  truss  element  capable  of  large  displacement  is  commonly  used 
for  cables.  The  derivation  [3]  of  the  stiffness  matrix  for  this  element  is  briefly 
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discussed  here  so  as  to  be  differentiated  from  that  in  chapter  3.  Fig.  2.8  presents 
the  straight  element  nonlinear  truss  in  its  local  axis. 


Fig.  2.8  Straight  nonlinear  truss  element  used  for  cables 


From  Eqn.  2-1,  the  component  of  horizontal  displacement  is, 

du 


du  , 

U:   + OS  -  U: 

1       ds  ' 


ds 


-ds 


(2-20) 


Similarly,  the  component  of  displacements  in  the  y  and  z  directions  are 

(2-21) 


dv  dv  J 

v-  +—-ds  -  v.-      =     — -  ds 
1      ds  '  ds 


dw  . 

W:  +— —  dS  -  W; 

1       ds  ' 


dw 

77 


ds 


(2-22) 


The  resultant  axial  strain  along  the  member  due  to  the  above  displacements  can  be 
computed  by  determining  the  deformed  element  length  ds.  The  deformed  length 
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due  to  displacement  in  the  x  direction  is  simply, 


ds        =     ds  +  — ds 

x  ds 


(2-23) 


The  deformed  length  due  to  displacements  in  the  y  and  z  directions  are  shown  in 
Fig.  2.9. 


[<*>"♦(>)] 


[(*>»♦(**-.)] 


Fig.  2.9  Deformed  length  due  to  displacements 


The  axial  strain  from  the  above  in  the  x,  y  and  z  directions  is  given  in  Eqns  2-24, 
2-25  and  2-26. 


Applying  the  binomial  expansion 


^    =    |+ifl2-V^ 


2  8 
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£x    = 


C7.S  <9w 


(2-24) 


<fr2  + 


(dv 
ds 


.2\X 


ds 


ds 


ds 


1  + 


2  \ 


-1 


V         *       '    J 


(2-25) 


<is2  + 


<  ^ 


-<fc 


ds 


ds 


(        /.     n2> 


1  + 


V  / 


(2-26) 


and  neglecting  the  higher  order  terms  Eqns  2-25  and  2-26  become 


h    =    l  +  2 


-i  = 


fa^ 


ydSj 


(2-27) 
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2     ds 


-  1     = 


'  dw 
ds 


(2-28) 


The  total  axial  strain  in  the  element  is  the  summation  of  the  strain 
expressions  in  Eqns  2-26,  2-27  and  2-28.  This  summation  is  referred  to  as  the 


large  deflection  strain-displacement  equation. 


du      \_ 

ds  +  2 


dv]2  +  ifdwV 

k  d  s  J  2  y  d  s  / 


(2-29) 


If  the  nodal  displacements  are  known  for  a  finite  length  L  of  the  nonlinear  truss 
element,  Eqn  2-29  becomes, 


e    - 


uJ-ui  ,  1 

+  — 

L  2 


(vj-vA 

{     L     J 

2      1 

+  — 
2 

(  Wj    -  Wj  > 

I       L       ) 

(2-30) 


Using  Castigliano's  theorem  for  the  force-displacement  relationship,  the 
expression  for  the  strain  energy  in  the  element  is  given  in  Eqn  2-31. 

U    ~    i/I-A      =    lfle2ds  (2-31) 

Substituting  the  expression  for  strain  and  expanding  the  brackets  we  obtain 


u    = 


AE 


fdu\2      a../-a..\2 


ydS; 


du 

ds 


dv 

yds; 


du 
ds 


ydSj 


(2-32) 


1    dv 


4{ds 


(dv\2f^..\2 


^ds 


dw 

yds  j 


(dw^ 


ds  t 


ds 


34 


Ignoring  the  fourth  order  terms  and  substituting  the  notation  used  in  Eqn  2-30,  we 
obtain  a  new  expression  for  U. 


V    =    ^J 

2   ' 


fuj-u^2 


L      ) 


Uj-Uj 


VJ     vl 


v     L     J 


Uj  -  «,- 


(yti  -  wt  ^2 


I     £     ) 


ds 


(2-33) 


Note  that, 


T     =     —(«,-«,) 


(2-34) 


Where  T  is  the  tension  force  in  the  truss  element.  Solving  the  integral  we  get  Eqn 
2-35. 


2       AE 


"    -    £("!->)    +#("/  -»Xv-»)    *%bl-lfa-lf       (M5) 


Differentiating  the  Eqn  2-35  with  respect  to  the  nodal  displacements  gives  the 
stiffness  coefficient  associated  with  that  displacement.  The  expression  for 
tension,  from  Eqn  2-34  above  is  substituted  into  Eqn  2-35.  The  assumption  is 
made  that  the  tension  in  the  cable  is  approximately  a  constant. 
du  AE 


dUj 


("'■  -  uj) 


(2-36) 


dU  T  (  v 


(2-37) 
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du 

d\Vj 


(w,  -  Wj) 


(2-38) 


dU  AE  i  s 

^7   ■   — (»;-«/) 


(2-39) 


du 


fh-,) 


(2-40) 


dw : 


=     j(Wj-Wi) 


(2-41) 


In  matrix  form,  the  geometric  stiffness  of  a  nonlinear  truss  can  be  written  as  given 
in  Eqn  2-42.  The  first  matrix  term  below  is  the  conventional  elastic  stiffness  of  a 
truss  element.  The  second  term  is  referred  to  as  the  geometric  stiffness  term.  This 
geometric  stiffness  accounts  only  for  initial  cable  stress  and  not  for  nonlinear 
strain  -  displacement  relationship. 
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Connector  Element 

The  load  on  the  traffic  signal  system  is  applied  by  the  traffic  signal  heads 
and  the  traffic  signs.  These  are  attached  to  the  system  via  connectors.  The 
connectors  have  translational  and  rotational  degrees  of  freedom.  In  this  they  are 
identical  to  the  beam  element  discussed  previously.  The  connectors,  however, 
undergo  large  displacements.  These  displacements  affect  how  the  stiffness  of  the 
connector  is  added  to  the  structure's  overall  stiffness.  The  connector  stiffness 
matrix  presented  next  is  a  beam  element  stiffness  accounting  for  large 
displacements.  The  element  presented  here  and  adopted  for  this  dissertation  was 
derived  by  Chen  and  Agar  [11].  The  element  stiffness  is  broken  into  the  elastic 
component  and  the  geometric  component.  The  elastic  stiffness  is  as  given  in  the 
earlier  described  beam  element,  the  geometric  stiffness  is  set  up  in  terms  of  the 
natural  displacements  (deformation)  of  the  beam  shown  in  Fig  2.1 1.  The 
derivation  of  the  geometric  stiffness  is  given  in  detail  in  Ref.  [11].  The  derivation 
presented  in  Ref.  [11]  rely  predominantly  on  work  done  by  Bathe  and  Bolourchi 
[12].  In  this  dissertation,  only  the  definition  of  the  coefficients  of  the  geometric 
stiffness  and  its  application  is  described.  The  degrees  of  freedom  for  the  large 
displacement  beam  element  used  as  connectors  is  identical  to  that  of  regular 
beams.  The  resulting  beam  end  forces  corresponding  to  these  degrees  of  freedom 
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are  shown  in  Fig.  2.10.  To  accurately  determine  the  end  forces  another  set  of 


forces  referred  to  as  the  natural  forces  or  basic  beam  forces  must  be  defined.  The 
natural  beam  forces  are  the  forces  developed  in  the  beam  due  to  relative 
deformation  of  the  beam  ignoring  rigid  body  motions.  These  forces  are  shown  in 
Fig.  2.11. 


Deformed  Position 


/     ui 


Original  Position 


Fig.  2.10        Joint  displacements  and 
corresponding  end  forces 
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Weak  Axis  Bending 


Strong  Axis  Bending 


Fig.  2.11        Relative  deformations  and  associated 
forces 


From  Fig.  2.10,  the  nodal  degrees  of  freedom  are, 
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(2-43) 


The  nodal  displacements  in  Eqn  2-43  consist  of  the  two  major  components  the 
rigid  body  component  and  the  natural  displacements  (i.e.  relative  deformation) 
component.  From  Ref.  [13]  six  rigid  body  modes  and  six  natural  displacement 
modes  completely  defines  any  displacement  that  may  be  seen  in  a  large 
displacement  beam.  The  identified  modes  are  as  given  in  Fig  2.12. 
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Rigid  Body  Modes       Sequence  of  modes        Natural  Modes 
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Fig.  2.12  Rigid  body  and  natural  modes  for  large 
displacement  beams 
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From  Fig.  2-12  it  can  be  seen  that  natural  mode  1  can  be  defined  by  the  change  in 

length  of  the  member  i.e.  axial  deformation.  The  rest  of  the  modes  of  deformation 

are  bending  deformations.  These  can  be  defined  by  the  end  rotations.  The  natural 

displacements  (i.e.  relative  deformations)  that  are  necessary  for  the  computation  of 

the  end  forces  can  therefore,  be  summarized  as  given  in  Eqn  2-44.  Notations  are 

shown  in  Fig.  2.11. 
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The  computation  of  {dN  }  in  Eqn.  2-44  is  not  intuitively  obvious.  For  illustration 
of  the  process  of  computing  the  relative  rotations,  an  assumption  is  made  here  that 
the  connector  is  initially  straight  as  shown  in  Fig.  2.13.  The  relative  rotational 
angles  are  the  total  nodal  rotations  obtained  from  analysis,  i.e.  static  or  dynamic 
analysis  minus  any  member  rigid  body  chord  rotations.  In  Fig.  2.13,  it  can  be 
seen  that  the  rigid  body  rotations  about  the  y,  x  and  z  global  axes  are  as  follows: 
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Deformed  chord  Position 


■'  Initial  chord  Position 


X  -  Z  plane  projection 


Fig.  2.13    Rigid  Body  Rotations 
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Using  the  cartesian  coordinate  system  shown  in  Fig  2.1,  clockwise  rotation  about 

the  y  axis  is  considered  positive  and  counter-clockwise  rotation  about  the  x  axis  is 

positive.  Therefore  when  dx  is  positive  (|>y  gets  a  positive  sign  but  when  dy  is 

positive  (J)x  must  be  assigned  a  negative  sign.    The  relative  rotational  deformation 

is  then  computed  as  shown  in  Eqn.  2-45  below. 
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Eqn  2-45  above  assumes  that  the  nodal  rotational  displacements  (with  the  nodal 
designation  i  or  j)  have  been  resolved  to  the  local  coordinate  system  but  the  chord 
rotations  ((J))  are  in  global  reference.  Nodal  rotations  are  measured  relative  to  the 
original  orientation  of  the  beam.  Relative  rotations  are  measured  relative  to  the 
chord  of  the  beam.    /'  is  the  deformed  length  of  the  beam  and  /  is  the  original 
length  of  the  beam.  To  check  the  validity  of  Eqn.  2-45,  consider  the  chord 
rotation  and  bending  deformation  of  the  next  connector  example.  In  Fig.  2.14,  for 
ease  of  illustration  the  assumption  is  made  that  deformation  exist  only  in  the  2  - 


dimensional,  x  -  z  plane. 
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Chord  of  deformed  beam 


Undeformed  Position 


y  axis 


4>r      Rigid  body  chord  rotation 

9'      Relative  deformation 

©       Nodal  rotation  1 


Deformed  beam 


Fig.  2.14     Relative  Nodal  Rotations 


Note  that  the  nodal  rotation  at  node  (i)  is  clockwise  hence  positive,  at  node  (/')  the 
nodal  rotation  is  counter-clockwise  hence  negative.    The  chord  rotation  is 
positive  at  both  ends.  The  resultant  relative  rotation  at  node  (i)  is  smaller  in 
magnitude  but  with  the  same  sign  as  the  nodal  rotation,  i.e.  clockwise.    The 
resultant  relative  rotation  at  node  (j)  is  bigger  in  magnitude  but  with  the  same  sign 
as  the  nodal  rotation,  i.e.  counter-clockwise.  This  conforms  with  results  from 
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applying  Eqn.  2-45.  The  natural  end  forces  of  the  beam  element,  as  shown  in  Fig 
2-11,  can  be  obtained  by  multiplying  the  above  relative  deformation  vector  by  the 


basic  beam  stiffness,  [kNE]. 
K}    =    KfilK  } 


(2-46) 


Substituting  the  matrix  terms  into  Eqn.  2-46,  Eqn  2-47  is  obtained. 
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The  end  force  vector  computed  from  Eqn  2-47  can  be  seen  in  Fig.  2.1 1. 
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The  coefficients  of  geometric  stiffness  are  obtained  from  the  basic  end  forces. 
The  geometric  stiffness  matrix  for  a  large  displacement  beam  according  to  Chen 
and  Agar  [1 1]  is  shown  in  Eqn.  2-49  below. 
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The  terms  in  the  matrix  above  are  as  follows: 


a  =  6/7  5/, 


b=  M     II,        c=  Mz   //,        d=P\l_. 


(V'JA-    e=M'/1-    f=p'xo 


8  *  [Mzi  +  Mzj)/6'        h  -  [Myi  +  Myj)/6>        '  =  2  W/15,        j  =  ^jf^i 


k=  M    .11, 


y,ll 


l=M7.jl,        m=Mt/2,         n=  Pl/30, 


(2-50) 


46 
Note  that  /  is  the  length  of  the  beam  and  the  other  variables  are  as  defined  in  Eqn 

2-48.  The  final  expression  for  the  stiffness  of  the  connector  element  with 

reference  to  the  element  local  axes  is, 

W  -  M*[*J  (2"51) 

The  above  stiffness  in  global  coordinates  can  be  obtained  by  applying  the  standard 
coordinate  transformation. 


Transformation  Matrix 

The  translational  transformation  matrix  can  be  obtained  by  using  the  k  -  node 
method.  This  method  is  used  in  the  SSTAN  program  developed  at  the  civil 
engineering  department,  University  of  Florida.  The  transformation  matrix  is  a  3x3 
matrix.  The  terms  of  the  transformation  matrix  [T]  is  obtained  as  follows.  The 
first  1  x  3  row  entries  represent  a  unit  vector  lying  along  the  member  local  axis, 
dx   dy   dz 

T'T'T 

/  =         length  of  the  member 

dx,  dy  and  dz  are  similar  to  Ax,  Ay,  Az  in  Fig.  2.7  and  similar  to  those  given  on 
page  40.  Another  unit  vector  must  be  defined  along  a  node  in  the  plane  of  the 
element  but  not  along  the  member.  This  node  is  referred  to  as  the  k-node. 
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dxk    dyk    dzk 


The  second  1  x  3  row  in  the  transformation  matrix  is  obtained  from  the  cross 
product  of  the  k-node  unit  vector  and  the  member  unit  vector.  The  third  1  x  3  row 
is  obtained  by  getting  the  cross  product  of  the  previous  member  unit  vector  (in  the 
first  row  of  the  transformation  matrix)  and  the  unit  vector  gotten  in  the  first  cross 
product  (now  in  the  second  row  of  the  transformation  matrix).  The  resultant  3x3 
matrix  [T]  is  used  to  transform  nodal  displacements  from  global  to  local 
coordinate  systems. 

For  the  case  of  the  beam  with  large  displacement,  in  addition  to  the  initial 
transformation  matrix  [T]  described  above,  it  is  often  advantageous  to  use  another 
transformation  matrix  [R]  based  on  Euler  angles.  This  is  because  this 
transformation  matrix  formulation  can  better  keep  track  of  the  twisting  of  the 
element  than  the  former.  A  combination  of  both  transformation  matrices  could  be 
adopted  with  the  first  transformation  [T]  used  to  transform  element  properties 
from  the  local  coordinate  system  to  the  initial  member  global  position.  The 
second  transformation  [R]  can  then  be  used  to  transform  from  the  initial  global 
position  to  the  member's  deformed  position.    The  transformation  matrix,  [R], 
used  in  the  context  just  described  is  described  next.  The  Euler  transformation 
matrix  [R]  is  also  explained  in  the  context  of  a  numerical  time  step  solution, 
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which  is  reflective  of  its  usage  in  chapter  3. 


Fig.  2.15  Euler  angles  for  transformation 


u 


w 


original  length  of  the  beam, 
axial  relative  displacement, 
relative  displacement  in  the  y  direction, 
relative  displacement  in  the  z  direction. 


The  Euler  angles  in  Fig.  2.15  can  be  computed  as  shown  in  Eqn.  2-52  and  Eqn.  2- 
53.  Note  that  u,  v  and  w  above  are  computed  using  the  cumulative  nodal 
displacements  in  the  case  of  dynamic  time  step  solution. 


a  =  tan 
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7TY  (2-52) 


a  -1  w 

P     =     tan = (2_53) 


j{l  +  u)2  +v: 


The  final  Euler  angle  that  is  needed  is  the  angle  of  twist,  y.  The  angle  of  twist  can 
be  obtained  from  the  rotational  displacement  about  the  local  x  axis.  It  is  assumed 
that  the  nodal  incremental  displacements  have  been  transformed  to  the  element's 
initial  local  coordinate  system  using  just  the  initial  k-node  transformation  [T]. 
The  angle  of  twist  is  cumulated  at  each  iteration  using  the  change  in  nodal 
rotations. 

y'    =    /~    '  +-]cosj8cosaf  AG   .+A0    .  )  +  cos/Jsinaf  A9   .+A0    .)  +  smB\A0  .+A0    .11 
2[  \      xi  xj)  \      yi  yjj         ry      u  ZJjj 

(2-54) 
y'  =         the  angle  of  twist  at  time  t 

y'  +a'  =         the  angle  of  twist  at  time  t  +  At 

A0(axisHnode)     =         change  in  nodal  rotation  about  the  axis  {x,  y  or  z}  at  node  {i 

orj}. 
a  and  P  are  as  computed  from  Eqn.  2-52  and  Eqn.  2-53. 
The  final  transformation  matrix  [R]  is  presented  in  Eqn.  2-55. 
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(2-55) 


The  final  transformation  matrix  for  the  large  displacement  beam  is  therefore: 

[A]    =    [R][T] 

[A]      =  final  transformation. 

Applying  the  transformation  matrix,  the  stiffness  of  the  member  in  terms  of  the 

global  nodal  displacements  is 
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(2-56) 


The  last  property  that  needs  to  be  defined  is  the  damping  matrix.  The 
damping  matrix  for  a  structural  system  is  easier  defined  proportional  to  either  the 
mass,  the  stiffness  or  both.  There  are  two  classical  damping  matrices.  The 
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Rayleigh  damping  and  the  Caughey  damping  [14].  The  Rayleigh  damping  matrix 

is  described  below.  With  Rayleigh  damping,  the  damping  matrix  is  defined  as 

follows. 

[C]    =    a0[M]    +    ax[K]  (2-57) 

ao  and  a,  are  damping  constants  depending  on  the  controlling  modal  frequencies 
and  the  corresponding  specified  modal  damping  ratio. 

°°  =  ^7^"  ai  =  *Wj  (2"58) 


£  =         modal  damping  ratio 

C3j        =         modal  frequency  corresponding  to  mode  i. 

GJj        =         modal  frequency  corresponding  to  mode  /. 

Controlling  modes  are  often  the  higher  modes  corresponding  to  the  number  of 

degrees  of  freedom  for  the  system. 

For  nonlinear  analysis  the  eigenvalues  are  never  obtained  as  part  of  the 
solution.  There  needs  to  be  a  way  to  estimate  the  natural  frequency  of  the  system 
prior  to  a  rigorous  nonlinear  analysis.  For  the  structural  system  defined,  an 
estimate  of  the  natural  frequency  of  the  system  can  be  taken  as  either  the  natural 
frequency  of  the  cables  on  infinitely  rigid  supports  or  the  natural  frequency  of  the 
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poles  without  the  cables.  Both  estimates  are  given  below. 

For  the  catenary  cable  subjected  to  self  weight,  no  initial  pre-tension  and  a 

specified  sag,  the  natural  frequency  [3]  is 


tan     =     nnj^-  n=l,2,3,...  (2-59) 


n         =  modes  to  investigate 

g         =  acceleration  due  to  gravity 

f  =  sag  at  mid-span  of  the  cable 

/  =  span  of  the  cable 

For  cables  where  the  sag  is  assumed  to  be  zero  the  above  equation  ceases  to  be  a 

good  estimate.  Therefore  for  a  messenger  cable  with  initial  pre-tension  and  a  near 

zero  sag  ,  the  natural  frequency  can  be  estimated  as 


/     V   m 


For  the  poles,  taken  as  cantilever  beams,  the  natural  frequency  can  be 
estimated  as  follows  [9]: 

..    =    2*L^W  (2-61) 
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<»-^ 
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Matrix  representation  of  the  element  properties  have  now  been  developed. 

The  solution  techniques  that  will  be  used  to  determine  the  time  history  response  of 

the  discretized  structure  is  considered  next. 


Dynamic  Solution  Methods 

The  equation  of  motion  [14]  for  a  typical  structure  is 

[M]{d]  +  [C]{d]  +  [K  ]{d]  =  {P}  (2-62) 

[  M  ]  =         assembled  mass  matrix  of  the  structure. 

[  C  ]  =         assembled  damping  matrix  of  the  structure. 

[  K  ]  =         assembled  stiffness  matrix. 

|^j  =         acceleration  vector. 

id)  =         velocity  vector. 

{d)  =         displacement  vector. 

{  P  }  =  applied  load  vector. 

For  nonlinear  analysis  direct  time  integration  is  used.  This  means  that  the  solution 
to  the  equation  of  motion  is  determined  at  discrete  times  ft  i=0,l,2...).  The  initial 
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conditions  at  the  beginning  of  the  time  steps  are  always  known  and  it  is  assumed 

that  some  variation  in,  and  self  consistent  relationship  between,  the  displacements, 

velocities  and  accelerations  occur  during  the  time  step.  The  methods  of  solution 

available  for  dynamic  analysis  vary  depending  on  the  assumptions  made  for  the 

variation  in  displacements,  velocities  and  accelerations  within  the  time  step.  As 

earlier  mentioned  these  methods  are  classified  as  either  explicit  or  implicit. 

Explicit  Methods 

A  disadvantage  of  using  explicit  method  is  that  they  are  conditionally 
stable.  The  more  common  explicit  solutions  are  Central  Difference  method  and 
Runge-Kunta  method.  In  explicit  methods  it  is  assumed  that  the  displacement  and 
velocity  at  time  t  +  At  are  independent  of  the  acceleration  at  t  +  At . 
Consequently,  the  displacement  at  the  end  of  the  time  step  is  projected  from  the 
equation  of  motion  at  the  beginning  of  the  time  step.  In  particular,  the  elastic 
forces  and  damping  forces  within  the  time  step  are  computed  using  displacements 
and  velocities  at  the  beginning  of  the  current  time  step  and  the  beginning  of  the 
previous  time  step.  For  details  of  the  available  explicit  methods,  the  writer 
recommends  Ref.  [7],  [9],  [14].  For  the  cable  analysis  it  is  important  to  apply  the 
load  in  small  increments,  it  is,  therefore,  more  convenient  to  perform  an 
incremental  analysis.  Implicit  methods  of  analysis  allows  us  to  rewrite  the 
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equation  of  motion  in  terms  of  the  change  in  the  displacement,  velocity  and 

acceleration. 


Implicit  Methods 

The  incremental  steps  are  measured  relative  to  a  reference  state  at  the 
beginning  of  the  time  step.  The  equation  of  motion  given  in  Eqn  2-62  can  then  be 
rewritten. 

[M]{Ad(tj\  +  [C]{Ad  (t)}  +  [K]{Ad(t)}  =  {AP(t)}  (2-63) 

The  new  shape  of  the  cable  at  any  time  (t)  can  be  obtained  as  follows 

d(t  +  At)  =  d{t)  +  Ad{t)  (2-64) 

The  solution  of  Equation  2-63  above  can  be  obtained  by  a  number  of  different 
implicit  methods.  The  more  common  methods  are  the  Adam-Stoermer  method, 
Newmark's  Beta  method  and  Wilson  Theta  method.  These  methods  actually  all 
have  the  same  form  and  can  be  easily  inter-changed  with  a  manipulation  of  the 
parameters  associated  with  each  method.  The  Newmark's  Beta  is  used  by  the 
writer  and  is  discussed  next. 
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Newmark  Beta  Method 

Newmark  Beta  method  [14]  is  based  on  the  following  equations  written  in 
terms  of  incremental  parameters. 


{Ad}=  At{d(t)}+   (I  -^}  {<*(*)}  +p{d(t  +  A*)} 


(2-65) 


{Arf}      =      y-[p(r)+i'(r+AO}] 


(2-66) 


Solve  Eqn.  2-65  for  {</(;  +  At)}  and  Eqn.  2-67  is  obtained. 

{<(<♦*<>]  -  jfrM-fem-^rm  (2"67> 

From  Eqn.  2-67  the  change  in  acceleration  within  the  time  step  can  be  expressed 
as  follows: 

{Ad}     =     {d(t+At)}-{d{t)}     =     -^{Ad}-^-{d(t)}^{d(t)}      (2-68) 

Substituting  Eqn.  2-67  into  Eqn.  2-66  and  the  expression  for  the  change  in 
velocity  can  be  rewritten  as, 

M  -  3^o-^w}-i^.{*w}  (2"69> 
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Substituting  Eqn.  2-65,  2-66  and  2-68  into  2-63  and  rearranging  terms  we  have 

the  final  equations  for  the  Newmark's  Beta  method. 

[K~(t)]{Ad}    =    [a  f(t)}  (2-70) 


[*«]    -    J^[M]+2j^[Cit)]  +  [K{')]  (2-71> 


{AP(t)}     =     {AP(t)}  +  [tAJ-{d(t)}  +  ±-{d{t)} 


pAt*        J       2/3 


[c(0] 


'^m+^m 


(2-72) 
The  solution  procedure  using  Newmark's  Beta  method  may  then  be  summarized 
as  follows: 

•  From  initial  conditions,  calculate  effective  stiffness  (Eqn.  2-71)  and 
effective  incremental  load  vector  (Eqn.  2-72). 

•  Solve  for  incremental  displacement  (Eqn.  2-70)  and  cumulative 
displacement  (Eqn.  2-64). 

•  Solve  for  incremental  velocity  (Eqn.  2-69)  and  velocity  at  t  +  At. 

•  Impose  total  equilibrium  at  the  end  of  the  time  step  to  obtain  the 
acceleration  at  the  end  of  the  time  step. 

{d(t  +  At)}  =  {n'  +  ")}-{M<^)} -{/.«+*)}  (2.73) 
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fd         =         damping  force  vector 

fe         =         elastic  force  vector. 
•  Repeat  the  process  for  the  next  time  step. 

The  parameter  (3  reflects  the  variation  of  the  acceleration  within  the  time 
step.  There  are  various  values  that  could  be  used  for  p.  For  P  =  1/6  we  have  the 
linear  acceleration  method.  For  P  =  1/4  we  have  the  constant  acceleration  method. 
And  finally  for  P  =  1/8,  the  acceleration  is  constant,  equal  to  the  beginning  value 
over  the  first  half  of  the  time  step  and  constant,  equal  to  the  ending  value  over  the 
last  half  of  the  time  step.  When  p  is  taken  as  zero  the  Newmark's  Beta  method 
degenerates  to  the  Euler  integration  scheme.  Similarly  if  p  is  taken  as  Vz 
Newmark's  method  becomes  identical  to  the  central  difference  integration  and  to 
backward  difference  when  P  is  unity.  The  constant  acceleration  method  is 
unconditionally  stable  .  This  means  that  the  solution  will  not  become  unbounded 
as  the  integration  progresses.  With  the  numerical  integration,  there  are 
inaccuracies  due  to  round-off  or  due  to  time  increment  being  too  large  to  represent 
the  forcing  function  or  the  structural  response.  The  errors  can  be  redistributed  at 
each  time  step  by  using  the  Newton-Raphson  method. 


CHAPTER  3 
THEORETICAL  DEVELOPMENT 


Introduction 


In  chapter  2,  the  typical  cable  element  commonly  found  in  literature  was 
discussed.  The  formulation  of  a  different  2  node  straight  cable  element  adopted 
for  this  dissertation  is  presented  in  this  chapter.  This  formulation  includes  the 
equation  of  motion  and  solution  to  the  equation  for  the  traffic  signal  system. 

The  equations  derived  in  this  chapter  are  based  on  the  Cartesian  coordinate 
system  shown  in  Fig.  2.1.  The  deformation  parameters  u  in  the  x  direction,  v  in  the 
y  direction  and  w  in  the  z  direction  are  deformations  relative  to  the  initial  (null) 
state.  The  initial  states  for  the  cables  are  taken  as  either  taut  or  slack.  The  taut 
cable  has  an  initial  pretensioned  cable  state.  The  load  acting  on  this  cable  at  the 
initial  configuration  (time  t  =  0)  is  the  self  weight  of  the  cable  and  the  initial 
tension.  The  cable  with  this  null  state  (condition  of  zero  deformation),  in  the  case 
of  the  traffic  signal  system,  is  referred  to  as  the  messenger  cable.  The  messenger 
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cable  null  state  has  an  initial  axial  strain  but  an  assumed  zero  sag.  The  catenary 

cable,  in  the  case  of  the  traffic  signal  system,  has  a  slack  null  state.  The  catenary 

cable  in  this  case  is  presumed  to  have  zero  pretension  (except  for  tension  from  self 

weight)  and  a  specified  sag.  The  default  sag  if  none  is  specified  is  5  %.  Two 

additional  reference  configurations  are  of  importance.  These  are  the  present 

configuration  at  time  t  and  the  next  configuration  at  time  t  +  At.    The  incremental 

displacement  at  the  nodes,  therefore  is  the  difference  between  the  coordinates  of 

the  nodes  at  time  t  and  at  time  t  +  At.  The  cumulative  displacement  at  the 

beginning  of  the  time  step  is  the  difference  between  the  coordinates  at  time  t  and 

the  null  state.  Similarly  the  cumulative  displacement  at  the  end  of  the  time  step  is 

the  difference  between  the  coordinates  at  time  t  +  At  and  the  null  state.  The  null 

configuration  is  the  base  reference  configuration.  Both  the  messenger  and 

catenary  cable  are  idealized  by  the  2-node  non-linear  cable  element  derived  in  the 

following  section.     The  main  difference  between  this  cable  element  and  that 

derived  in  chapter  2  is  the  inherent  assumption  in  the  strain  -  displacement 

relationship. 


Two  Node  Nonlinear  Cable  Element 


The  shape  function  of  the  element  adopted  is  given  in  Eqn.  3-2. 


^1 


Fig.  3.1  Cable  Element 
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The  element  properties  for  this  cable  element  are  derived  directly  in  the 
global  coordinate  system.  An  explicit  derivation  of  a  2  node  cable  element  from 
first  principles  using  some  novel  manipulations  of  typical  energy  expression  is 
presented  next.  The  derivation  presented  is  based  on  the  work  presented  by 
Henghold  and  Russell  [16]  and  Haritos  and  He  [17].     Ref.  [17]  presented  a  total 
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Lagrangian  formulation  of  the  incremental  equation  of  motion  for  cable 

movement.  This  writer  has  adopted  this  approach  with  a  deviation  in  that  the 

cable  tension  used  for  the  geometric  stiffness  is  cumulated  in  each  time  step  in  an 

updated  Lagrangian  fashion.  This  is  observed  to  be  more  convenient  in  terms  of 

numerical  implementation  of  the  element.  Ref.  [16],  however,  uses  the  method  of 

perturbation  of  the  virtual  work  expression  for  a  cable  element  to  derive  the 

equation  of  motion  of  the  cable.  Though  this  is  very  different  from  the  approach 

adopted  in  this  dissertation,  the  writer  has  chosen  to  adopt  the  reference's 

approach  of  defining  the  element's  displacement  field  (shape  function)  in  such  a 

way  as  to  derive  the  element  stiffness  directly  in  terms  of  the  deformed  global 

coordinates  as  opposed  to  the  more  classical  approach  of  using  a  local  coordinate 

system  with  a  local  to  global  transformation.  The  shape  function  is  defined  to 

relate  the  coordinates  of  any  generic  point  along  the  cable  element  to  the  nodal 

coordinates. 
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(3-1) 


In  terms  of  a  convenient  notation  Eqn.  3-1  can  be  rewritten  as 
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shown  in  Eqn.  3-2. 

M'WM 


(3-2) 


N 


coordinate  of  an  arbitrary  point  s  along  the  cable 

coordinates  of  the  nodal  points,  n  is  the  degrees  of  freedom  shown 

in  Fig.  3.1. 

shape  function 


The  interpolation  function  for  a  linear  displacement  field  is  given  below  [18]. 
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(3-3) 


Applying  the  principle  of  virtual  work,  the  expression  for  total  work  done  for  a 
single  element  at  the  time  t  +  At  is  as  follows. 

8W=\Ls-j-T—ds+     £    8—T—^5—T—ds_^A—T-ds=Q 


(3-4) 
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6W      =         sum  of  the  internal  and  external  virtual  work  done  on  the  element. 

qs        =         gravity  and  inertia  forces  acting  along  the  arc  of  the  cable. 

Fj        =         applied  nodal  load  -  j  is  the  degrees  of  freedom,  (j  =  l...n). 

f  d       =         damping  force  vector. 

5ds      =         vector  of  small  virtual  variations  in  coordinates. 

Note  that,  though  for  convenience  the  writer  has  chosen  to  omit  the  brackets  { } 

used  as  vector  notation,  the  displacement  and  forces  in  Eqn.  3-4  and  subsequent 

equations  are  vectors.  Note  also  that  the  straight  line  overstrike  ," ,  indicate  values 

or  properties  at  time  t  +  At.  The  superscript ,  T,  indicate  the  transpose  of  the 

vector. 

The  last  term  in  the  Eqn.  3-4  is  the  equivalent  internal  virtual  work  due  to 

stress  o  at  time  t  +  At.  This  is  referred  to  as  the  virtual  strain  energy.  The 

variation  in  the  coordinate  at  time  t  +  At  is  equivalent  to  the  variation  in  the 

difference  between  the  coordinates  at  time  t  and  time  t  +  At.  It  is,  therefore  useful 

to  express  the  coordinate  at  time,  t  +  At,  as  follows. 

KLai-Wi  +  Kl  (3-5) 

{Ads}     =     [N]{Adn]  (3-6) 

Adn      =         change  in  nodal  displacement  vector. 
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Ads      a         change  in  displacement  vector  of  any  generic  point  along  the  cable. 
Substituting  Eqn.  3-2  into  Eqn.  3-4  yields  Eqn.  3-7. 


SW  =  ^[N]T8AdnTqsds+    ^[NfdAdJ  F.  -  £[N]T  SAd/fjds-fi  A  8eT  a  ds 


(3-7) 


Kirchhoff  stress. 


Lagrangian  strain  of  deformed  cable. 


A  nonlinear  strain-displacement  relationship  is  adopted,  therefore  the  strain  is  as 
given  in  Eqn.  3-8. 


s      2 


<** 


Kdsoj 


(3-8) 


as 

ds„ 


stretched  length  of  a  cable  segment  at  time  t  +  At 
unstressed  length  of  a  cable  segment  at  time  t  =  0 


In  terms  of  the  shape  function  expression  in  Eqn.  3-3,  Eqn.  3-8  may  be  rewritten 


as, 


S      2 


ddsTdds 


ds         ds 
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(3-9) 


*,4{R}Vr>]KH 


(3-10) 
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(3-11) 


Note  that  N'  is  the  differential  of  the  shape  function  with  respect  to  the  generic 
distance,  s,  along  the  cable.  The  analysis  approach  described  in  chapter  2  is  that 
of  incremental  displacement.  It  is  therefore  important  to  define  the  corresponding 
incremental  strain.  Subtracting  the  nodal  coordinates  at  time  t  from  the 
corresponding  coordinates  at  time  t  +  At,  the  incremental  nodal  displacement  is 
obtained.  Similarly,  subtracting  the  expression  for  strain  (Eqn.  3-8)  at  time,  t, 
from  a  corresponding  expression  at  time  t  +  At  yields  the  change  in  strain. 


Ae    =    - 
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{  ds  J 

(3-12) 


67 


coordinate  of  an  arbitrary  point  s  along  the  cable  at  time  t  +  At 


Ads     =    ds~ds 


(3-13) 


Solving  for  dsin  Eqn.  3-13  and  substitute  into  Eqn.  3-12  results  in  Eqn.  3-14. 


Ae 


(dAd^ 
{    ds    J 

T 

+ 

{  ds  ) 

T 

{    ds    ) 

fdAd^fdd,) 
(    ds         {   ds 


(3-14) 


Substituting  the  notation  used  in  Eqn.  3-10  into  3-14,  Eqn.  3-15  is  obtained. 

A£     "     \{AdnT[N'}T["']Adn  +  d/[N>]T[N']Adn+Ad/[N']T[N']dn} 

(3-15) 

The  last  two  terms  in  Eqn.  3-15  are  equal,  therefore  the  above  equation  can  be 
simplified  further  as  shown  in  Eqn.  3-16. 

A£     =     \{*dnT[N']T["']*d»}+dJ[N']T[N']*dn  (3-16) 


The  virtual  work  of  internal  forces  as  given  by  the  last  term  of  Eqn.  3-4  is  given 
below. 
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SU     =     AjfcST  ads 


&8TT*ds  (3-17) 


Applying  the  unit-load  method  [19],  the  virtual  strain  energy  expression  in  Eqn.  3- 
17  can  be  rearranged  to  be  in  terms  of  the  stiffness  for  the  cable  element. 
SU     =     AE^8eTeds  (3-18) 

The  strain  in  the  cable  element  at  time  t  +  At  is  equivalent  to  the  strain  at  time  t 
plus  the  change  in  strain  given  in  Eqn.  3-16.  The  variation  in  strain  at  time  t  +  At 
is  equivalent  to  the  variation  of  the  change  in  strain. 

x-t         sa«t  (3-19) 

oe        =     oAe 

Substituting  Eqn.  3-16  into  3-19  and  then  into  3-18  the  equation  given  in  Eqn.  3- 
20  is  obtained. 

SU=AEl^{Ad/[NY[N']}8Adn+^{sAdJ[N']T[N']}Adn+d/{N']T[N']SAdnJ    e  ds 

(3-20) 
6Adn    =         variation  of  incremental  nodal  displacement  vector. 
Again,  note  that  the  first  two  terms  in  Eqn.  3-20  are  equivalent. 

8U  =  AEj^{Ad/[N']T[N']]8Adn+   d/[N'f [N']8AdnJ  Ids 

(3-21) 

From  Eqn.  3-21,  it  can  be  noted  that  the  term,  6Adn,  is  common  to  all  the  terms  of 

Eqn.  3-4,  and  hence  cancels  out.  Substituting  the  term  for  strain,~e,  Eqn.  3-21 

reduces  to  an  expression  in  terms  of  stiffness. 
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e     +     Ae  (3-22) 


X 


r* 


£  +  {Jr{AdnT[N'V  [N']Adn}+     dnT  lN'V  [N']Adn))d 

(3-23) 


Rearranging  and  collecting  sub-expressions  in  Eqn.  3-23,  Eqn.  3-24  is  obtained. 

fs      =      KEAdn      +      KLAdn      +      RE  (3"24) 


K  E         =        A   E    j     e    [  N  '  ]T   [  N  '  ]ds  (3-25) 


kl    =    ^4({^/[Ar1r[Ar']}+{d/[Ar-f[iV']})rx(I{M/[Ar'f[^]}+{rf/[Ar'f[iV']})A 

(3-26) 

R  E         =        AE\[N']T[N']e     d       ds  (3-27) 

0 


KE  is  the  referred  to  as  the  initial  stress  matrix  and  RE  is  the  force  vector  due  initial 
stress.  Note  that  the  term  AEe  is  equivalent  to  tension  in  the  cable  at  the 
beginning  of  the  time  step.  For  t  =  0,  this  term  is  the  pretension  in  the  element. 
For  subsequent  time  steps  the  tension  term  can  be  computed  as  follows, 


T     =    T     A      + 
t  t-At 


«H-*  (3-28) 


Note  that  KL  is  dependent  on  unknown  displacements.  This  is  the  non-linear 
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component  of  the  element  stiffness,  [K].  This  can  not  be  computed  at  the  start  of 
analysis  and  must  be  substituted  by  an  initial  approximate  stiffness.  Let  this 
approximation  be 


KA      =     AEJ  \dJ[N']T[N'])T  x\dHT[N']T[N'])ds 


(3-29) 


The  expression  in  Eqn.  3-29  was  achieved  by  simply  setting  Adn  in  Eqn.  3-26  to 
zero.  If  the  shape  functions  for  the  2-node  element  is  substituted  into  the  Eqn.  3- 
25,  3-26,  3-27  and  3-29,  the  following  expressions  for  KE,  KL  and  KA  on  the  next 
page  are  obtained.  See  Appendix  A  for  complete  derivation  of  the  expression. 
The  following  notations  are  used. 
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The  expression  for  KE  is  as  follows: 
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(3-30) 
initial  cable  element  length. 
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Eqn.  3-33  can  be  deduced  from  Eqns.  3-25  and  3-27. 
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(3-33) 


(3-34) 


Equation  of  Motion 

The  above  portion  calculated  the  strain  energy  portion  of  Eqn.  3-4.  It  is 
now  necessary  to  consider  the  other  terms  in  Eqn.  3-4.  The  kinetic  energy  term  is 
the  first  term  in  Eqn.  3-4.  This  is  in  terms  of  the  variable  qs  defined  in  Eqn.  3-35 
below. 


is    +   ft 


(3-35) 


gs 
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distributed  load  along  the  element, 
gravity  load  due  to  self  weight, 
inertial  forces. 


Rewriting  the  first  term  in  Eqn.  3-7  yields  Eqn.  3-36. 
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SK    =     fw^Nff^ds    -     jm[N]T[N]{d}ds 
0  0 


(3-36) 


6K 


cable  self  weight  per  unit  length, 
cable  mass  per  unit  length, 
kinetic  energy  due  to  inertia  forces. 


{>.} 


nodal  acceleration. 


h[  is  an  index  array  based  on  the  coordinate  system  adopted.  It  is  used  to 
determine  how  the  weight  of  the  element  is  added  to  the  structure. 
The  self  weight  acts  in  the  negative  z  direction. 
h[  =  {0    0    -1} 
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acceleration  due  to  gravity. 
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The  second  term  in  Eqn.  3-7  and  Eqn.  3-4  is  the  applied  nodal  load. 
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(3-39) 


This  is  the  computed  nodal  load  due  to  dead  weight  of  cladding  and/or  traffic 
signal  heads  and/or  wind  load.  Since  the  cable  element  is  a  two  node  element,  the 
nodal  distribution  of  a  body  load  is  an  even  split  between  the  two  adjacent  nodes. 
The  third  term  in  Eqn.  3-7  is  the  damping  term.  This  is  obtained  as  discussed  in 
chapter  2.  The  equation  of  motion  for  the  cable  element  is  therefore  given  as 
shown  in  Eqn.  3-40. 


76 

[*]ftWc]kW*]K}  =  {'/}*{*.}-{*«}  0-40) 


The  solution  to  Eqn.  3-40  above  can  be  obtained  by  a  myriad  of  numerical 
integration  schemes  mentioned  in  chapter  2.  When  the  cable  element  is  assembled 
with  the  other  element  types,  i.e.  beams  and  connectors  TheRE  term  in  Eqn.  3-40 
above  can  be  replaced  by  a  similar  term  FE.  Where  FE  is  the  initial  stress  vector  for 
all  the  elements  in  the  structure. 


[*3{'.Mc]{J.W*]{".}  -    {',}♦{*,}-{'*}  (3"41> 

It  should  be  noted  that  the  self  weight  vector  remains  the  same  at  the  beginning  of 
the  time  and  at  the  end.  The  change  in  force  due  to  the  self  weight  is  zero.  When 
the  wind  load  is  considered  the  nodal  force  vector  would  include  the  wind  force. 
At  the  beginning  of  each  time  step  the  incremental  displacement  vector  is  zero  and 
the  velocity  and  acceleration  vector  are  as  obtained  from  the  previous  time  step 
analysis.  The  initial  conditions  of  the  dynamic  analysis  is  obtained  from  setting 
both  the  displacement  and  velocity  vectors  to  zero  and  obtaining  the  acceleration 
from  equilibrium  of  the  structure. 

U    \=      iFj}  +  {''hi'*}  (3-42) 

l      ■  J  [M    ] 
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The  use  of  the  incremental  displacement  facilitates  the  use  of  implicit 

numerical  integration  schemes  for  solution  to  Eqn.  3-41.  Chapter  2  discusses  the 

Newmark  Beta  method.  This  method  in  its  modified  form  was  adopted  by  the 

writer.  The  modification  adopted  is  presented  in  the  following  section.  This 

modification  is  referred  to  as  the  alpha  modification. 

Modified  Newmark  Beta  Method 

Due  to  the  inherent  problem  of  "blow  up"  in  the  analytically  unstable 
traffic  signal  cable  system  and  the  fact  that  the  unconditionally  stable  Newmark 
Beta  method  possesses  no  intrinsic  numerical  damping,  Wood  et  al  [20]  proposed 
the  use  of  an  alpha  factor.  This  factor  introduces  a  positive  artificial  damping. 
This  modified  Newmark' s  method  with  positive  artificial  damping  of  higher 
vibrational  modes  is  observed  not  only  to  reduce  the  "blow  up"  effect  but  also 
reduce  the  cumulative  error  which  inevitably  occurs  in  a  numerical  integration 
procedure. 

Applying  the  alpha  modification  to  Eqn.  3-40  results  in  an  averaging  effect 
of  the  inertia  forces  at  the  beginning  and  at  the  end  of  the  time  step. 

(3-43) 
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W,+a,  =  btA,+{4+A,-hi+A, 


(3-44) 


Note  that  if  a  is  zero  Eqn.  3-43  reduces  to  Eqn.  3-40.  Apply  Newmark's  method 
based  on  Eqns.  3-45  and  3-46  to  Eqn.  3-44.  For  convenience  the  nodal 
designation  is  dropped  for  Eqn.  3-45  to  Eqn  3-50.  It  should  however  be 
understood  that  these  are  nodal  displacement,  velocity  and  acceleration. 


t+At 


dl+[(l-y)At]dt  +  (YAt)dt  +  At 


(3-45) 


dt  +  At=dt  +  (A'^t  +  \(lA-ii)(At)2}dt  +  \p^') 


t  +  At 


(3-46) 


Solving  Eqn.  3-46  for  acceleration  at  time  t  +  At  yields  Eqn.  3-47. 


/  ,       "\ 


t+At 


pAt 


2Adt 


pAtd* 


2/J 


(3-47) 


\"r       j 


Substituting  Eqn.  3-47  into  3-45  yields  Eqn.  3-48. 


t+At 


-!—Ad 

PAt  ' 


Of-'K 


•-TjK 


(3-48) 


Following  the  steps  described  in  chapter  2  and  substituting  Eqn.  3-47  and  3-48 
into  Eqn.  3-43,  the  final  expressions  for  the  effective  stiffness  and  load  vectors  are 


obtained. 
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(3-49) 


Note  that 


{*d,}={dl  +  Al}-{dl} 


[K]     =     (l-a)[M] 


jAt: 


+     [C\ 


'-21 ' 

KPAtj 


+      [*] 


(3-50) 


Collecting  terms  and  rearranging  Eqn.  3-49  gives  Eqn.  3-51. 


(3-52) 


(4}  -  W„^t«i{fe!W,^->-  «)W.-»R}.W(HW.^-')K}. 


(3-53) 


When  a  is  zero  the  above  described  Newmark's  method  reduces  to  that  described 
in  chapter  2.  The  values  of  P  and  y  are  as  described  in  chapter  2.  If  the  alpha 
modification  is  to  be  used,  the  values  of  the  constants  are  given  below. 
When  P  =  Va  and  y  =  Vi  for  the  regular  Newmark's  method,  the  constant 
acceleration  method  is  obtained.  This  method  is  a  second  order  solution  with  zero 
artificial  damping  and  unconditional  stability.  The  intent  of  the  alpha  factor  is  to 
insert  positive  artificial  damping  while  keeping  the  property  of  unconditional 
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stability.  For  a  second  order  method  with  positive  artificial  damping  and 

unconditional  stability  the  following  condition  must  be  satisfied. 

p>^>-    &    a<-    &    a  +  y>-  (3-54) 

2      4  2  2 

Wood  et  al  [18]  in  their  paper  suggest  the  use  of  the  following  as  giving  the  best 
results. 

a  =  0.1  P  =  0.3025     y  =  0.6 

a  =  -0.1  p  =  0.5  y  =  0.6 

a  =  -0.1  P  =  0.3025     y  =  0.6 

The  writer  suggests  the  use  of  a  =  0.1      P  =  0.3025     y  =  0.6  for  the  traffic 
signal  cable  system  as  it  yields  the  best  artificial  damping  effect  for  the  system. 
When  the  incremental  displacement  is  obtained  from  Eqn.  3-51,  in  order  to 
maintain  dynamic  equilibrium,  this  computed  incremental  displacement  needs  to 
be  refined.  The  reasons  for  the  need  of  refinement  are  as  follows; 

a)  The  stiffness  matrix  at  time  t  +  At  is  dependent  on  the  recently 
computed  displacements  values. 

b)  The  {FE  }  term  in  the  equation  of  motion  contain  values  dependent 
on  the  nonlinear  cable  stiffness  and  nonlinear  connector  stiffness. 

The  Newton-Raphson  method  is  adopted  to  refine  the  computed 
incremental  displacements.  Application  of  the  Newton  -  Raphson  is  discussed 
next. 
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Modified  Newton  Raphson  Method 

A  graphical  representation  of  the  steps  involved  in  the  modified  Newton 
Raphson  method  is  given  in  Fig.  3.2.  There  are  two  additional  modifications 
made  to  the  regular  non-recursive  Newton  -  Raphson  procedure  [9],  [14].  One 
has  to  do  with  the  nonlinear  loading  and  the  other  the  nonlinear  stiffness  of  the 
cable  elements.  Note  that  both  the  loading  and  the  stiffness  are  dependent  on  the 
incremental  displacements  being  computed.  The  writer  has  taken  these  into 
account  in  modifying  the  Newton  -  Raphson  procedure  adopted. 
Recall  from  Eqn.  3-41  that, 


{Fj }     =         Applied  nodal  loads.  These  consist  of  dead  weight  of  cladding 
and/or  wind  loads  applied  at  the  nodal  points.  Discussion  of 
loading  on  the  traffic  signal  cable  system  is  presented  in  chapter  4. 

{RE}  =  Equivalent  nodal  force  vector  for  member  stress  at  time  t.  This 

accounts  for  the  stresses  in  the  structure  at  the  beginning  of  the  time 
step. 
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Force 


AF. 


Displacement 
Fig.  3.2  Newton  -  Raphson  Iteration 
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The  complete  procedure  to  refine  Ad0  shown  in  Fig.  3.2  on  the  previous 
page  is  given  to  illustrate  the  modified  Newton-Raphson  method  : 
1.         Initialize  data. 


KQ     =    \k\        From  Eqn  (3  -52) 
P    -    Ip\        From  Eqn  (3 -53) 


AR      =     P 
o 


Ad 


-1 


'       =     K  ~l  AR        FromEqn{3-5\) 


In  computing  the  above  initial  values,  AP  is  assumed  to  be  zero  in  Eqn. 
3-53.  It  is  assumed  that  if  At  is  taken  small  enough  the  approximation  error 
introduced  by  this  assumption  is  negligible. 

2.  Compute  the  instantaneous  cable  stiffness  at  a  displacement  Ad0  and  set  it 
equal  to  k± . 

3.  The  incremental  resisting  force: 

A/l      =      KLc    Ado  (3_56) 

Where  the  secant  stiffness  (K,)sec  can  be  approximated  by  Eqn.  3-57. 
(Ki)  =     £o_lA  (3-57) 

V     ' /sec  1 
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AR    +  K.  Ad 
af1    =    — ^ 1 2-  (3-58) 


4.  The  residual  force  at  the  end  of  the  first  iteration  is  as  shown  in  Eqn.  3-59 
ar.    =    p    -    AFj  (3-59) 

5.  Ad,  is  computed  as  in  step  1. 

Adl     =     k^a^  (3-60) 


The  sign  of  the  displacement  computed  in  Eqn.  3-60  is  negative  because 
the 

residual  force  is  negative. 

6.  Repeat  step  2,  compute  k2 . 

7.  Repeat  step  3,  compute  new  resisting  force. 

«o(Ado  +  Adi)+«2(Ado+Adi)  (3-6D 

•••    AF2    =    3 

8.  Repeat  step  4,  compute  the  residual  force. 

a  ^2    =    p    -    a  f2 

9.  Using  Eqn.  3-60  above  compute  Ad2.  The  whole  procedure  is  continued 
until  the  residual  force  is  smaller  than  some  tolerance  value.  The  true 
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displacement,  therefore,  is  as  given  in  Eqn.  3-62. 
Adn    =    ^o+Adl  +  Ad2  (3_62) 


There  are  many  convergence  checks  that  could  be  used.  Some  authors 
suggest  a  check  on  subsequent  displacement  computations;  that  is 

— l-    <     £  i  =  0,1,2  ...  (3-63) 


A  d 
n 


where  £  is  a  prescribed  tolerance.  The  problem  faced  when  convergence  is  strictly 
based  on  the  convergence  of  the  nodal  displacement  is  that  of  different  rates  of 
convergence  for  rotational  displacements  and  translational  displacements.  The 
units  and  range  of  values  for  rotations  differ  greatly  compared  to  translations.  As 
a  result  the  translational  displacements  may  converge  while  the  rotational 
displacements  may  not  simultaneously  converge.  The  writer  has  chosen  a 
convergence  criterion  based  on  the  work  done  by  the  residual  forces. 


{".} 


T 

Ad. 

i 


{AP}T   Adf 


<    I       i- 0.1.2...  (3-64) 


The  use  of  the  modified  Newmark's  method  and  the  above  described  Newton- 
Raphson  method  in  the  entire  solution  procedure  is  described  next. 
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Solution  Algorithm 

Based  on  the  theory  developed  in  the  previous  sections  of  this  chapter  the 
solution  to  the  equation  of  motion  given  in  Eqn.  3-40  was  obtained  as  follows. 

1.  Compute  the  effective  load  vector  as  given  in  Eqn.  3-53.  Note  that  AP 

is  given  by  the  right  hand  side  of  Eqn.  3-41.  For  the  first  time  step  analysis 
AP  is  taken  to  be  zero.  At  the  end  of  each  time  step  Pt+ At  is  computed  and 
the  actual  change  in  load,  AP,  for  that  time  step  is  obtained.  This  change  in 
load  is  used  as  the  estimate  for  the  next  time  step  analysis.  The 
acceleration  at  the  beginning  of  each  time  step  is,  however,  obtained  by 
imposing  total  equilibrium  at  the  end  of  the  previous  time  step  and  using 
the  true  load  in  the  equilibrium  equation. 

2.  Compute  the  effective  stiffness  matrix  using  Eqn.  3-52.  The  initial 
effective  stiffness  matrix  at  time  t  =  0,  uses  [K]  =  [KE]  +  [KA], 

where  [KE]  and  [KA]  are  as  defined  in  Eqn.  3-25  and  3-28.  At  subsequent 
times,  the  effective  stiffness  matrix  uses  [K]  =  [KE]  +  [KJ.  The  nonlinear 
stiffness  [KJ  is  computed  as  given  in  Eqn.  3-26. 
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3.         Solve  for  change  in  displacements. 

{»'.}  ■  [*l-{'},  (3-65) 


4.  The  Newton  Raphson  iteration  scheme  is  used  to  refine  the  incremental 
displacement  obtained  in  Eqn.  3-65.  The  Newton-Raphson  procedure 
gives   a  better  estimate  of  { Adn } .    This  is  done  as  shown  in  steps  1  -  9  on 
pages  83  -  85. 

5.  Compute  the  final  condition  at  the  end  of  the  time  step. 

{'.La,     =     {'.},♦{".}  ^ 


I   "it  +  At 


(3-67) 


(l-a)[Af] 


(3-68) 


6.         Set  the  incremental  displacement  of  the  next  time  step  as  follows  and  repeat 
the  above  steps 

•     K}  -  M 
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The  dynamic  time  step  analysis  is  concluded  when  the  free- vibration  response  dies 

down.  The  complete  analysis  gives  the  full  dynamic  response  history  for  the 

transient  dynamic  loads  applied  and/or  the  equivalent  static  solution  of  the  dead 

load  applied.  For  the  static  equivalent  solution  the  "die-down"  criteria  is  as 

follows. 


![*]{a^}2     ![k]K} 


p-i 


q  -  p 


s     £ 


(3-69) 


This  means  that  the  strain  energy  of  the  structure  is  averaged  over  several  time 

steps  and  compared  to  the  average  of  the  next  several  time  steps.  When  the 

difference  between  the  two  strain  energies  computed  is  negligible,  die-down  is 

achieved. 

/,  p  and  q  are  beginning  and  ending  time  step  cycles  for  strain  energy  computation. 


Solution  Flowchart 
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[      Start       J 


Read  in  material  and  geometric  properties 


Form  [M],  [C]  and  [K]  for  each  member 


Read  in  load  and  assemble  structure  properties 


Obtain  effective  stiffness.  Eqtn  3-52 
Obtain  effective  load.  Eqtn  3-53 


Solve  Eqtn  3-51  for  displacement 


Compute  new  effective  stiffness 
using  instanteous  tangent  stiffness 


Newton-Raphson  iteration  loop: 


Compute  resisting  force.  Eqtn  3-56 


Compute  residual  force.  Eqtn  3-59 


Continued  on  next  page 
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Next  time  step 


Compute  final  condition  at  the  end 
of  the  time  step.  Eqtn  3-66  -  Eqtn  3-68 


Fig  3.3  Flowchart  of  Solution  Algorithm 


CHAPTER  4 
LOADING 


Introduction 


The  loads  on  the  traffic  signal  system  can  be  split  into  two  categories.  The 
dead  load  and  the  wind  load.  The  dead  load  is  the  summation  of  the  weight  of  the 
cable,  weight  of  the  signal. heads.  The  dead  load  assumed  to  be  acting  on  the 
traffic  light  system  is  as  given  in  Ref.  [21].  The  wind  load  acting  on  the  structure 
is  discussed  in  this  chapter. 

Two  models  of  wind  loading  were  developed.  The  variable  area  wind  load 
model  and  the  variable  velocity  dynamic  wind  load  model.  The  variable  wind 
load  area  model  was  largely  developed  in  Ref.  [21].  The  application  of  this  model 
is  slightly  different  in  this  research.  The  variable  velocity  model  is  a  simplified 
true  wind  simulation.  Both  models  and  their  applications  are  discussed. 

Variable  Area  Wind  Load 


The  wind  speed  is  assumed  to  be  a  constant  for  this  model.  The  wind  load 
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is  taken  to  act  at  the  centroid  of  the  traffic  signal  heads  and  the  traffic  signs.  As 
the  wind  impacts  the  surfaces  of  these  objects  the  wind  pressure  is  taken  to  act 
normal  to  the  these  surfaces.  The  resultant  normal  force  is  resolved  into  the 
vertical  uplift  force  and  the  horizontal  drag  force. 


Wind 


Light  Signal  Head 


Wind  Uplift  Force 
(WL) 


Wind  Drag  Force 
(WD) 


Fig.  4.1    Resolved  Wind  Load  on  Signals 


The  wind  pressure  is  obtained  directly  from  Bernoulli's  general  law  for  an 
ideal  fluid. 


- =      0.002558VZ      lb  s  I  ft1 


(4-1) 


V 


mass  density  of  air. 
along-wind  velocity.  (ft/sec} 


p         =         wind  pressure  (lbs/ft). 
If  wind  gust  is  to  be  taken  into  account  we  introduce  a  factor  Fz.  The  gust  factor 
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varies  with  height.  The  gust  factor  at  any  height  is  given  by  Eqn.  4-2  [24]. 

(4-2) 


Ft       =       ^30 


,  3  Qx  0.0625 


\  z  J 


F30       =         gust  factor  at  30  ft,  recommended  to  be  1 .30. 

z  =         height  of  the  centroid  of  the  signal  heads  above  ground. 

For  the  structures  being  investigated,  the  gust  factor  is  approximately  1 .3  for 
height,  z,  of  17  ft.  The  final  pressure  equation  becomes, 

p    =    0.002558  (1.3  V)2  ^4"3^ 

The  resolved  forces  are  computed  from  Eqns.  4-4  and  4-5  respectively.  These 
Eqns.  introduces  two  coefficients,  the  lift  and  drag  factors.  These  are  used  to 
account  for  the  aspect  ratio  and  orientation  of  the  wind  surface. 
WD    =    CDpA  (4-4) 

WL      =      CL  P  A  (4-5) 

CD       =         drag  coefficient 

CL       =         lift  coefficient 

A         =         area  of  the  wind  surface 

The  drag  and  lift  factors  account  for  the  reduction  in  the  area  impacted  by 
wind  as  the  signal  head  rotates.  The  drag  and  lift  coefficients  in  Eqns  4-4  and  4-5 
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are  given  in  Ref.  [21]  and  [22]  and  are  shown  in  Fig.  4.2. 
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Fig.  4.2  Lift  and  Drag  Coefficients 
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The  expressions  for  the  drag  and  lift  coefficients  shown  in  Fig.  4.2  is  given  by 

Eqns  4-7  and  4-8  [23]. 

CD     =     1.2-1.08xl0-2a-4.16xl0"5«2  (4-7) 

CL=232xlO~2a   -  1.87x  10_4a2  -  2.23x  10-6a3  +  8.61  x  10"8a4  -  7.79x  10~10a5 

(4-8) 

Due  to  the  fact  that  the  signal  heads  are  free  swinging,  the  surface  area  impacted 

by  the  wind  at  any  given  time  is  dependent  on  the  tensions  in  the  cables  and  the 

corresponding  deformations  in  the  cables.  For  dynamic  analysis,  the  time  step 

method  facilitates  the  computation  of  the  angle  a  at  every  time  step  and  hence  the 

calculation  of  the  new  area  impacted  by  the  wind  in  each  time  step. 

Variable  Velocity  Dynamic  Wind  Load 

The  wind  speed  in  this  model  is  taken  to  be  the  summation  of  a  mean  wind 
speed  value  and  a  fluctuating  (or  gust)  component. 

V(t)     =     V+v(t)  (4-9) 

The  mean  wind  speed  value  is  the  same  as  in  the  previous  wind  load  model. 
The  fluctuating  component  can  be  obtained  from  experimental  data  or  be 
artificially  generated.  The  generation  of  artificial  wind  gust  data  is  often  the 
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easier  choice.  A  vast  body  of  literature  exists  on  the  methods  of  artificially 

generating  the  gust  component  of  wind.  Only  a  brief  presentation  of  the  method 

adopted  for  this  research  is  made.  Further  information  can  be  obtained  from  the 

following  references;  Davenport  [24],  Chen  and  Ahamdi  [25],  Iannuzzi  and 

Spinelli  [26],  Solari  [27]  [28]  and  Haritos  [29]. 

The  gust  component  of  wind  can  be  represented  by  a  random  process  with 
a  substantially  Gaussian  distribution.  This  means  that  the  process  has  zero  mean 
and  a  unit  variance.  The  process  is  taken  as  stationary  (i.e.  invariant  with  time) 
and  ergodic  (i.e.  the  statistics  of  the  sample  are  the  same  as  that  for  the  ensemble). 
These  assumptions  are  important  because  the  theoretical  derivation  of  the  power 
spectrum  from  time  records  is  only  possible  for  a  stationary  process  and  for  the 
wind  sample  generated  to  be  realistic,  the  ensemble  of  random  wind  signals  must 
be  assumed  spatially  homogenous. 

The  longitudinal  field  of  turbulence  for  a  wind  trace  must  conform  to  the 
following  spectral  density  suggested  by  Davenport  [24]: 

AkVJ       v2 


S»    =  ">        A     %  (4-10) 

(l+X2) 


CO 


Sv        =         power  spectrum  density  of  horizontal  gustiness  at  height  z10. 

z,0  is  standard  reference  height  and  is  set  equal  to  10  m.  (33  ft.). 


97 


io       - 


CO 
K 


mean  hourly  wind  velocity  at  the  reference  height  (ft/sec), 
frequency  (cycles/sec). 

drag  coefficient  for  ground  disturbance  at  reference  height. 
See  table  from  Ref.  [24]  below. 


X    = 


4000  co 


V, 


10 


CO 


— —    =  wave  number  in  cycles/  ft 
k  10 


Type  of  surface 


Open  terrain  with  very  few  obstacles:  open  grass  or 
farmland  with  few  trees,  hedgerows  etc. 


Terrain  uniformly  covered  with  obstacles  30-50  ft 
in  height :  residential  suburbs,  small  towns,  small 
fields. 


Terrain  with  large  and  irregular  obstacles  e.g.  large 
cities,  very  broken  country  with  many  tall  trees. 


Drag  coefficient 


0.005 


0.015 


0.050 


There  are  two  methods  of  simulating  the  necessary  random  fields 
required,  one  is  based  on  the  superposition  of  sine  and  cosine  functions.  This  is 
known  as  the  wave  superposition  method  [26]  and  the  other  is  based  on  the  use  of 
numerical  filters  on  sequences  of  random  numbers,  this  is  commonly  referred  to  as 
linear  filtering.  For  this  research  the  later  method  is  adopted.  The  random 
number  sample  generated  must  have  a  normal  distribution  given  in  Eqn  4-1 1  and 
shown  in  Fig.  4.3. 
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/w  rfc* 


-Xx2 


(4-11) 


Normal  distribution 
o  =  1 
x  =  0 


Distribution  of  all 
Sample  values 


Frequency  Density 


Fig  4.3  Distribution  of  sample  generated 


99 
The  random  number  generator  similar  to  that  proposed  by  Leva  [30]  is 

utilized  here.  This  is  a  normal  pseudorandom  number  generator.  The  numbers 

are  generated  using  an  acceptance-rejection  technique.  Each  number  generated  is 

checked  to  see  if  it  lies  within  an  acceptable  Gaussian  distribution  field. 

The  input  for  the  normal  random  number  generator  is  another  random 

number  generator,  usually  a  flat  spectral  random  number  generator.  The  random 

number  routine  used  could  be  the  implicit  rand()  routine  in  Fortran  77  library  or 

any  other  random  number  generator.  The  writer  has  chosen  to  use  a  random 

number  generator  provided  in  the  standard  C  library.  This  random  number 

routines  generate  numbers  in  the  range  [0,1]  with  a  flat  spectral  density 

distribution.    The  routine,  drand(),  generates  a  sequence  of  48-bit  integer  values, 

Xj,  according  to  the  linear  congruential  formula. 

Xi+l  =  (a  X,  +  c)modm 

The  default  values  for  the  Unix  Sun  system  for  the  constants  a,  c,  and  m  are 

a  =  5DEECE66D16  =  2736731631558 

c  =  B16=138 

m  =  248 
The  routine  drand  must  be  initialized  with  a  seed  value,  X;.  A  default  value  exists 
for  this  but  it  is  often  advised  not  to  use  this.  The  seed  value  is  provided  to  the 
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computer  by  the  user  with  the  routine  srand  (long  seedval).  This  routine  is  called 

once.  Subsequent  calls  of  drand  uses  the  value  last  number  returned  as  seed.  This 
ensures  the  data  is  truly  random  and  no  seed  is  used  twice. 

The  random  data  generated  represents  white  noise.  This  must  be  filtered  to 
conform  to  the  spectral  density  function  given  in  Eqn.  4-10.  To  do  this  another 
property  of  the  spectral  density,  the  covariance  function  must  be  considered. 
Iannuzzi  and  Spinelli  [26]  defined  the  covariance  function  using  the  Weiner- 
Khinchin  co-spectrum  equation.  Their  formulation  accounted  for  certain  along- 
wind  fluctuations  as  would  be  evident  if  two  objects  in  the  path  of  the  along-wind 
turbulence  were  separated  by  some  distance.  For  the  purpose  of  this  research  one 
object  in  the  stream  of  the  wind  is  considered,  hence  wind  analysis  of  one  signal 
process  is  sufficient  as  opposed  to  the  correlation  of  two  signal  processes  provided 
in  Ref.  [26].  The  covariance  function,  in  this  case  becomes  the  auto-covariance 
function  (Eqn.  4-12).  For  the  purpose  of  applying  numerical  filters  the  auto- 
covariance  function  shall  be  evaluated  at  five  different  lag  times  (t).  Values  of  x 
chosen  are  given  in  Eqn.  4-13. 

CXX  (*)=ls(°>)e~i2KOir  da>  =  ls((o)cos(2n(or)d(o  (4"12) 

0  0 


T     =     n  At         n--2,  -1,  0,  1,  2 


(4-13) 
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To  be  able  to  evaluate  the  integral  given  in  Eqn.  4-12  it  is  important  to 
understand  another  property  of  the  spectral  density  function  given  in  Eqn  4-10. 
This  has  to  do  with  the  frequency  content  of  a  typical  gust  wind.  Using  as  an 
example  a  mean  wind  velocity  of  60  mph  and  a  drag  coefficient  of  0.05,  the 
spectral  density  using  Eqn.  4-10  yields  Fig.  4.4  shown  below. 
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Fig.  4.4    Spectral  density  of  fluctuating  part  of  wind  (V=60  mph) 
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It  is  observed  that  the  spectral  amplitude  increases  rapidly  at  small 

frequencies  reaching  a  maximum  at  about  0.03  Hz.  After  which  there  is  a  rapid 

decay  of  spectral  energy.  By  0.2  Hz,  90%  of  the  spectral  energy  i*  passed  and  by 

1  Hz,  the  percentage  is  99%  .  Therefore,  for  the  integral  given  in  Eqn.  4-12  a 

frequency  integration  limit  of  0  -  1  Hz,  in  this  case,  represents  the  entirety  of  the 

spectrum.  The  expression,  Eqn  4-12,  could  be  solved  exactly  by  using  the  theory 

of  complex  contour  integration.  However  for  this  research  numerical  integration 

was  chosen  to  solve  the  integral.  Simpson's  method  is  used.  The  largest  positive 

frequency  that  would  be  the  limit  of  integration  is  determined  from  the  length  of 

record  of  wind  speed  to  be  generated  and  the  time  interval  of  the  wind  record 

required.  The  smallest  frequency  that  can  be  detected,  referred  to  as  the  frequency 

resolution  (or  spectral  resolution)  f0,  is  called  the  fundamental  frequency.  The 

largest  frequency,  fq,  is  called  the  Nyquist  frequency  or  the  folding  frequency. 

/      =    -  (4-14) 

J  O  T 

I 


f      =    HL.  (4-15) 

q  2 


T,        =         Length  of  record  (sec) 

N        =         Number  of  wind  speed  data  (Must  be  an  even  number). 

With  the  definition  of  the  auto-covariance  it  is  now  possible  to  make  the 
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white  noise  generated  by  the  random  generators  conform  to  the  spectral  density 

function  of  given  in  Eqn.  4-12.  Several  numerical  filters  have  been  used  in  the 

available  body  of  literature  to  do  this.  Auto-regressive  filters  were  chosen  for  this 

study  because  these  have  been  proven  to  need  the  least  computational  time  [26]. 

The  only  caveat  is  that  the  series  generated  using  a  time  step  of  less  than  0.1  sec 

tend  to  be  divergent.  For  situations  where  the  time  step  for  structural  response  is 

smaller  than  0.1  sec  (as  is  often  the  case),  linear  interpolation  schemes  with  wind 

data  generated  for  a  time  step  of  0.1  sec  should  be  used. 

Auto-regressive  Filters 

An  auto-regressive  filter  of  order  p  is  given  in  Eqn.  4-16 

p 

v(0  -  5Xv('-wA0+e(0  (4-i6) 

m=\ 


v(t)      =         new  stream  of  filtered  data  with  the  spectral  distribution 

in  Eqn.  4-12.  This  represents  the  fluctuating  portion  of  the  wind, 
a,,,  auto-regressive  parameters  computed  from  the  covariance  functions, 

t,  At     =         time  and  time  increment  of  wind  history, 
p         =         order  of  regression.  For  this  research  p  =  3  was  chosen. 
e(0     =      <Jxx(t)  (4-17) 


104 
ox        =         standard  deviation  of  the  filtered  random  process.  This  can  be 
computed  from  Eqn.  4-19. 
The  regression  parameters  can  be  computed  from  the  system  of  linear 
equations  obtained  from  Eqn.  4-18. 

Cxx{jAt)=%Cxx[(j-m)At]am      j  =  l,...P  (4-18) 

m  =  l 

P 

°>      =     Cxx(0)  +  ^amCxx(mAt)  (4.19) 

m=l 


Substituting  the  previously  selected  value  of  p  equal  to  3  into  Eqn.  4-16  we  obtain 

the  following  set  of  linear  equations. 

Cx x  ( At)  =      Cxx  (0)a1  +  Cx  x  (-At)a2  +  C^  ^  (-2Ar>3  (4.20) 

CZX(2A0=      Cxx(At)al+  Cxx(0)a2+Cxx(rAt)a3  (4-21) 

Czx  (3A/)=      CXX(2A0«1  +  CXX(At)a2  +  CXx(°>3  <4"22) 

The  parameters  a„  a2,  a3  can  be  evaluated  form  Eqns  4-20,  4-21  and  4-22. 
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White  Noise  Data 


Wind  Gust  Data 


Numerical  Filter 


Psf  -  Probablity  density  function 

SD  -  Spectral  density  function 

Fig.  4.5  Numerical  filtering  Process 


This  procedure  should  result  in  a  realistic  variable  wind  trace.  The  final 
step  is  to  check  if  the  spectral  density  function  of  the  variable  wind  speeds 
generated  in  Eqn.  4-16  is  the  same  as  that  given  in  Eqn.  4-12.  This  is 
accomplished  by  computing  the  variance  of  the  filtered  data  and  compare  this  with 
the  variance  of  the  spectral  density  function  in  Eqn.  4-12.  Both  variances  should 


be  identical. 


Since  dynamic  wind  generated  accounts  for  gust,  the  velocity 


computed  in  Eqn  4-9  is  used  in  Eqn  4-1  instead  of  4-3.  The  variable  wind  model 
equation,  therefore,  becomes 


p    =    0.002558  (V  +v(f)) 


(4-22) 


The  pressure  value  computed  from  Eqn.  4-22  is  therefore  used  in  the  lift  force  and 


drag  force  equations  as  in  the  previous  wind  load  model.  A  sample  of  an  actual 
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Test  Data 

A  wind  load  test  was  performed  on  a  traffic  signal  system.  Details  of  the 
tests  are  given  in  Ref.  [31].  The  results  from  these  experiments  are  compared  with 
predictions  from  the  program  developed  in  this  dissertation.  Most  of  the 
comparisons  of  results  are  discussed  in  chapter  6.  One  comparison  germane  to 
this  chapter  is  the  comparison  of  the  statistics  of  the  raw  wind  velocities  measured 
during  these  experiments  to  the  spectral  density  function  given  in  Eqn.  4-12. 
Before  comparison  can  be  made  it  should  be  noted  that  random  data  obtained 
were  measured  using  low  frequency  anemometers  of  which  examples  are  pitot 
static  tubes,  propellers  and  vane  anemometers.  In  this  case,  Gill  3-Cup 
anemometers,  Model  12102,  were  used.  The  readings  obtained  are  biased  by  the 
response  time  of  the  instrument  used.  In  addition,  wind  loading  recorded  were  not 
natural  wind  events  but  were  generated  by  propellers  of  air  boats.  The  wind  gusts 
generated  by  the  air  boats  were  therefore,  atypical  of  the  natural  wind  gusts  in 
terms  of  randomness  and  duration  of  the  events.  Three  different  tests  set-ups  were 
considered  (Fig.  4.7).  The  two-point  connection  signal  system,  the  one-point 
connection  signal  system  and  another  two-point  connection  signal  system  with 
near  zero  catenary  sag.  Only  results  from  the  first  test  set-up  are  considered  in 
chapter  6  to  verify  the  theories  presented  in  this  dissertation. 
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The  Results  for  test  set-up  #  3  are  used  in  this  chapter.  The  essence  of 
comparisons  in  this  chapter  is  to  determine  if  the  wind  loads  from  the  air  boat  wind 
simulation  could  be  generated  by  the  numerical  wind  simulation  presented  in  this 
chapter. 

Wind  load  was  simulated  in  stages.  The  first  stage  is  the  low  wind  speed 
range,  then  a  medium  wind  range  and  finally  a  high  wind  speed  range.  When  the 
wind  speed  ensemble  is  taken  the  spectral  density  distribution  is  reflective  of  that 
given  in  Eqn.  4-10  as  is  apparent  in  Fig.  4-8  below. 
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Fig.  4-8    Spectral  Density  (All  wind  speed  ranges) 
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The  wind  samples  collected,  however,  were  observed  not  to  be  ergodic. 

This  is  as  to  be  expected  because  wind  mean  velocity  was  forced  to  be  varying  in 

such  a  way  as  to  allow  for  the  partitioning  of  the  data  into  the  wind  speed  ranges 

mentioned  above.  The  spectral  density  plots  for  the  ranges  are  shown  in  Fig.  4.9  to 

Fig.4.11. 


Plot  of  the  Spectral  Density  Function  for  Low  wind  speeds  range 
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Fig.  4.9  Spectral  Density  (Low  Wind  Speeds) 
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Plot  of  the  Spectral  Density  Function  for  Medium  Wind  Speeds 
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Fig.  4.10  Spectral  Density  (Medium  Wind  Speeds) 


The  wind  speed  trace  generated  by  the  air  boats  in  the  low  and  medium  wind  speed 
range  obviously  does  not  conform  to  the  spectral  distribution  of  the  Fig.  4.4.  This 
would  have  to  be  taken  into  consideration  when  comparing  the  results  in  chapter  6, 
since  the  results  are  compared  in  ranges.  For  the  high  wind  speeds  some 
correlation  was  observed  with  the  spectral  density  distribution  given  in  Fig.  4.4  as 
can  be  seen  in  Fig.  4.11. 
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Plot  of  the  Spectral  Density  Function  for  High  Wind  Speeds 
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Fig.  4.11  Spectral  Density  (High  Wind  Speeds) 


There  apparent  conclusion  therefore  is  that  numerical  simulation  can  be  used  with 
some  certainty  of  accuracy  to  predict  the  air  boat  wind  fields  for  the  entire  wind 
speed  ensemble  but  would  not  simulate  the  gradation  in  wind  speed. 


CHAPTER  5 
PROGRAM  VERIFICATION  FOR  STATIC  LOADS 


Static  Problem  Definition 

Full  scale  static  load  tests  were  conducted  on  a  Two-point  traffic  signal 
system.  These  tests,  performed  at  the  University  of  Florida  Structural  Laboratory, 
are  used  to  verify  the  theories  presented  in  this  dissertation.  Details  of  actual 
testing  and  instrumentation  can  be  seen  in  Ref.  [31].  For  the  purpose  of 
subsequent  discussion  only  the  test  system  configuration,  the  sequence  of  load 
application  and  the  pertinent  test  results  are  of  importance.  Prior  to  the  static  load 
test  discussion  a  simpler  single  cable  structure  is  used  to  verify  the  cable  element 
developed. 

This  single  cable  structure  was  presented  in  a  paper  [32]  by  Baron  and 
Vendatesan.  The  results  from  this  paper  was  compared  in  Ref.  [7,  pg  113]  to 
verify  another  method  of  analysis.  Results  from  both  references  are  compared 
here  with  results  from  DYNASS,  the  program  developed  in  this  dissertation. 
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Single  Cable  Structure  Example 

The  single  cable  example  is  shown  in  Fig.  5.1  below.  The  area  of 
the  cable  is  0.227  in2  and  E  =  12  000  ksi.  The  corresponding  self  weight  is  0.7945 
lb/ft.  Though  the  self  weight  was  not  considered  in  the  static  methods  of  analysis 
presented  in  Ref.  [7]  and  [32],  the  weight  is  required  for  dynamic  analysis  and  is 
hence  included  in  DYNASS. 


50  ft. 


4  kips  4  kips      Initial  dead  load 

I  Excludes  self  wt. 

3  kips      Additional  Applied  Load 


Fig  5.1  Single  Cable  Problem 

For  comparison  in  DYNASS,  analysis  is  started  from  a  slack  no  load  state.  The  4 
kip  loads  are  then  applied  for  the  prestressed  state.  Finally,  the  additional  3  kip 
load  is  added.  The  comparison  of  results  is  only  done  for  the  additional 
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deflections  due  to  the  added  3  kip  load  and  the  stresses  in  the  cable  at  final 
equilibrium.  Deflection  is  compared  only  at  nodes  1  and  2. 


Table  5-1  Displacement  Comparison  for  Suspended  Cable 


Node 

# 

Component 

Ref.  [7] 
inches 

Ref.  [32] 
inches 

DYNASS 

inches 

%ofdiff. 

DYNAS- 

Ref.  [7] 

%ofdiff. 
DYNAS- 
Ref.  [32] 

1 

X 

26.00 

26.28 

25.47 

-2.04 

-3.08 

z 

52.36 

52.08 

49.86 

-4.77 

-4.26 

2 

X 

23.80 

23.40 

22.93 

-3.66 

-2.01 

z 

-49.11 

-48.48 

-47.22 

-3.85 

-2.60 

Table  5-2  Stress  Comparison  for  Suspended  Cable 


Member  # 

Ref.  [7] 
ksi 

Ref.  [32] 
ksi 

DYNASS 
ksi 

%ofdiff. 

DYNAS- 

Ref.  [7] 

%ofdiff. 
DYNAS- 
Ref.  [32] 

1 

52.96 

53.32 

53.50 

0.65 

0.34 

2 

48.52 

48.91 

49.04 

1.06 

0.27 

3 

55.18 

55.52 

55.93 

1.34 

0.73 

The  cable  element  in  DYNASS  appears  to  be  slightly  stiffer  than  that  assumed  in 
the  Ref.  [7]  and  [32].  It  should  however,  be  noted  that  inherent  in  the 
assumptions  in  static  solutions  in  Ref.  [7]  and  [32]  is  the  principle  of 
superposition.  This  principle,  according  to  some  schools  of  thought,  has  been 
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observed  to  overestimate  structural  response  for  non-linear  problems.  The 
analysis  in  DYNASS  applies  all  load  concurrently,  hence  does  not  use  the 
principle  of  superposition.  In  addition,  initial  pre-tensioning  of  the  cable  due  to 
its  self  weight  is  accounted  for  in  DYNASS  and  not  in  the  Ref.  [7]  and  Ref.  [32]. 

Static  Test  System  Configuration 

The  traffic  signal  system  tested  is  given  in  Fig.  5.2  also  available  in  Ref. 
[31  -  p  16].  The  poles  were  made  of  2  -  W  12x26  each  braced  in  the  plane  of 
loading  by  a  two  double  angle  section  (2L  -  3Vix3V6xV4  )  and  also  in  the  out-of- 
plane,  by  a  single  angle  section  (L  -  3x4x14).  A  lA  inch  zinc  coated  steel  wire  was 
used  for  the  messenger  cable  and  a  %  inch  zinc  coated  steel  wire,  for  the  catenary 
cable.  The  hanger/connectors  were  steel  pipes  \Vi  inch  in  diameter.  Actual  signal 
heads  were  not  used  but  the  weight  of  a  typical  3  light  (12"  lens),  polycarbonate 
material  plus  disconnect  box  signal  head  was  simulated.  Properties  of  each 
member  of  the  set-up  is  given  in  Table  5-3. 
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Table  5-3    Section  Properties  for  Members  shown  Fig.  5-2 


Member 

A 

(in2) 

12 

(in4) 

13 

(in4) 

J 

(in4) 

Strain  Pole  (W  12x26) 

7.65 

204 

17.3 

0.300 

In-plane  bracing  (2L  -  3V4x3V£xV4) 

3.38 

4.02 

1.59 

0.0386 

Out-of-plane  bracing  (L  -  3x4x1/4) 

1.69 

2.77 

1.36 

0.0386 

Pipe  Connector  (IV2  0) 

0.799 

0.31 

0.31 

0.620 

The  modulus  of  elasticity  (E)  used  for  all  sections  except  for  cables  was  29000 
ksi.  For  both  cables,  a  young's  modulus  of  23000  ksi  was  deemed  appropriate. 
The  coordinates  for  each  node  fed  into  the  program  is  as  given  in  Table  5-4. 
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Table  5-4     Dimensions  and  node  coordinates  for  Fig.  5.2 


Node 

X 

Y 

Z 

(in) 

(in) 

(in) 

1 

0 

0 

0 

2 

0 

0 

59.75 

3 

0 

0 

80.75 

4 

0 

0 

122.75 

5 

149 

0 

c=l* 

6 

298 

0 

c=r 

7 

449 

0 

C=l* 

8 

600 

0 

122.75 

9 

149 

0 

80.75 

10 

298 

0 

80.75 

11 

449 

0 

80.75 

12 

600 

0 

80.75 

13 

298 

0 

varies** 

14 

600 

0 

59.75 

15 

600 

0 

0 

16 

96 

0 

0 

17 

504 

0 

0 

18 

0 

72 

-10 

1      .9 

600 

72 

-10 

This  initial  coordinate  is  computed  by  the  program. 
Load  was  applied  at  varying  heights  for  different  tests. 
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Loading  Sequence  of  Tests  Performed 


There  were  in  all  3  static  tests  performed.  The  loading  sequence  and  initial 
stress  conditions  for  each  test  are  given  in  Tables  5-5,  5-6  and  5-7. 

Table  5-5    Loading  Sequence  for  Test  #  1 


Test  #  1 

Load  Stages 

Load  Applied  (lbs) 

Comments 

Vertical 

Lateral 

Initial  conditions 

Initial  messenger  cable 
tension  was  1.141  kips. 
Initial  catenary  cable 
tension  was  0.165  kips 
after  the  hanger  weight 
alone  was  applied. 

Stage  1 

50.4 

0 

Hanger  weight  plus 
concentrated  load.  Sag 
in  catenary  cable  was 
3.633%  after  dead  load 
application.  Messenger 
cable  was  attached  after 
dead  had  been  applied. 

Stage  2 

50.4 

29.0 

Point  of  application  of 
load  was  46.375  in. 
from  node  6. 

Stage  3 

50.4 

56.3 

Stage  4 

50.4 

83.6 

Stage  5 

50.4 

101.2 

Last  load  sequence. 
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Table  5-6    Loading  Sequence  for  Test  #  2 


Test  #2 

Load  Stages 

Load  Applied  (lbs) 

Comments 

Vertical 

Lateral 

Initial  conditions 

Initial  messenger  cable 
tension  was  1.127  kips. 
Initial  catenary  cable 
tension  was  0.167  kips 
after  hanger  weight 
only. 

Stage  1 

60.4 

0 

Hanger  weight  plus 

concentrated  load.  Sag 

in  catenary  cable  was 

3.743%  after  the  dead 

load  application. 
Messenger  cable  was 
attached  after  dead  load 
was  applied. 

Stage  2 

60.4 

29.0 

Point  of  application  of 

load  was  46.375  in. 
from  node  6. 

Stage  3 

60.4 

56.3 

Stage  4 

60.4 

83.6 

Stage  5 

60.4 

119.6 

Last  load  sequence. 
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Table  5-7    Loading  Sequence  for  Test  #  3 


Test  #  3 

Load  Stages 

Load  Applied  (lbs) 

Comments 

Vertical 

Lateral 

Initial  conditions 

Initial  messenger  cable 
tension  was  1.046  kips. 
Initial  catenary  cable 
tension  was  0.522  kips 
after  hanger  weight  only. 

Stage  1 

70.4 

0 

The  messenger  cable 
was  attached  before 
application  of  the 
concentrated  load.  Sag 
in  catenary  cable  was 
1,1 12%  after  dead  load 
application. 

Stage  2 

70.4 

29.0 

Point  of  application  of 
load  was  63.94  in.  from 
node  6. 

Stage  3 

70.4 

56.3 

Stage  4 

70.4 

83.6 

Stage  5 

70.4 

110.9 

Last  load  sequence. 

For  test  3  the  effect  of  sag  was  being  investigated  hence  the  drastic  change 
in  sag.  The  results  from  the  above  tests  and  their  comparison  to  the  predicted 
values  from  DYNASS  is  presented  next. 
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Test  Results  and  Predicted  Values 

The  test  results  from  the  tests  are  tabulated  against  the  values  predicted 
from  the  DYNASS  program.  Results  compared  are  displacements  of  nodes  6 
(catenary  cable)  and  10  (messenger  cable)  and  tension  in  the  cables  Since  a 
dynamic  theory  has  been  used  to  approximate  a  static  load  case,  plots  showing  the 
complete  structural  response  including  die-down  is  included  for  the  dead  load 
cases. 
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Discussion  of  Results 

To  best  model  the  actual  testing  procedure  as  presented  in  Ref.  [31].  The 
following  analysis  steps  were  adopted  by  the  writer.  For  Test  1  and  2  the  hanger 
was  attached  to  the  catenary  cable  alone  prior  the  additional  gravity  load  being 
attached.  Both  the  hanger  and  the  catenary  cable  were  then  attached  to  the 
messenger  cable.  In  the  analysis  of  the  Test  1  and  2  set-up.  The  hanger  self 
weight  was  applied  to  the  catenary  cable.  The  stresses  and  displacements  for  this 
load  case  was  obtained  and  served  as  the  base  reference  point.  When  the 
additional  gravity  load  is  applied  the  change  in  the  configuration  from  the 
reference  configuration  is  compared  with  what  was  measured.  For  the  application 
of  the  lateral  loads  the  prestressed  configuration  from  the  gravity  load  case  served 
as  input.  A  pre-tensioning  of  the  connector  element  is  done  to  ensure  that  the 
gravity  load  pick-up  by  the  messenger  cable  is  minimized.  This  procedure  is 
observed  to  follow  the  actual  testing  procedure  within  limits  of  experimental  error. 

Comparing  tension  computations  for  the  static  tests,  a  good  match  is 
observed  for  tests  1  and  2  with  the  maximum  percent  difference  of  3.82%.  For 
test  3,  however,  there  is  a  larger  difference  with  a  percent  difference  of  7.67%.  It 
should  be  noted  that  loading  procedure  for  test  3  was  changed  in  the  middle  of 
testing.  The  hanger  was  connected  before  the  additional  gravity  load  was  applied. 
This  caused  some  relaxation  in  the  catenary  cable.  The  near  zero  sag  for  test  3 
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could  also  be  a  factor  in  the  resulting  difference  between  the  measured  and  the 
computed  tensions,  this  is  because  absolute  zero  sag  in  the  catenary  cable  causes 
numerical  floating  errors  in  the  analysis  method  adopted.  Fig.  5.3  to  Fig.  5.5 
show  graphically  the  change  in  tension  relative  to  the  change  in  lateral  load. 

The  differences  in  displacements  measured  to  that  computed  was  more 
varied  and  less  consistent  than  with  the  tension  differences.  This  is  probably  due 
to  the  fact  that  the  relative  difference  in  displacement  is  an  average  of  half  an 
inch.  The  readings  in  the  lab  were  observed  to  fluctuate  more  than  half  an  inch. 
Hence,  the  difficulty  in  establishing  if  the  variations  are  due  to  non-conforming 
results  or  experimental  errors.  The  results,  nevertheless  validate  the  theories 
espoused  by  this  research.  Displacement  comparisons  are  seen  shown  in  Tables  5- 
17  to  Table  5-19.  The  dynamic  response  of  the  structure  is  best  seen  with  a  plot  of 
a  nodal  displacement  with  time.  The  structural  response  of  the  starting 
configuration,  i.e.  hanger  weight  load  case,  is  shown  in  Fig.  5.6  to  Fig.  5.8.  The 
regions  of  structural  instability  are  identified  in  each  plot  for  test  1  and  2.  The 
behavior  of  the  test  3  was  decisively  different  but  is  included  for  completeness. 
The  regions  identified  are  points  where  the  structural  response  had  been  observed 
to  diverge  when  a  larger  time  step  is  used  in  the  analysis. 
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CHAPTER  6 
PROGRAM  VERIFICATION  FOR  DYNAMIC  LOADS 


The  Structure's  Configuration 


A  full  scale  two-point  traffic  signal  system  was  tested  at  the  University  of 
Florida.  The  span  of  the  structure  was  typical  of  a  four  lane  road.  Details  of  the 
actual  tests  can  be  found  in  Ref.  [31].  The  system  set-up  and  node  numbering  is 
shown  in  Fig.  6.1.  The  strain  poles  used  were  type  N-V  prestressed  poles  with 
concrete  strength,  fc\  of  6000  psi.  The  catenary  and  messenger  cables  were  both 
%  in.  0  cables.  Cables  were  connected  using  strap  connectors  and/or  144  in.  0 
steel  pipes.  The  traffic  signals  mounted  were  2  -  5  head  light  signals,  12"  lens, 
Aluminum,  no  backplate  and  one-way  bracket.  The  wind  load  applied  was 
predominantly  in  the  out-of-plane  direction.  In-plane  wind  forces  are  neglected. 
The  devised  wind  loading  method  involved  the  use  of  air  boats  to  create  the  wind 
fields.  Wind  velocity  was  varied  between  10  mph  and  100  mph  to  simulate  gust. 
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Coordinates  of  Wind  Load  Test  System  Configuration: 

Table  6-1     Dimensions  and  node  coordinates  for  Fig.  6.1 
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Node 

Y 

z 

(in) 

(in) 

(in) 

1 

0 

0 

16 

2 

0 

0 

112 

3 

0 

0 

164 

4 

0 

0 

166 

5 

180 

0 

c=V 

6 

360 

0 

c=r 

7 

504 

0 

c=r 

8 

684 

0 

c=V 

9 

864 

0 

164 

10 

864 

0 

142 

11 

180 

0 

112 

12 

360 

0 

112 

13 

504 

0 

112 

14 

684 

0 

112 

15 

864 

0 

112 

16 

360 

0 

112-27.3 

17 

504 

0 

112-27.3 

18 

864 

72 

0 

This  initial  coordinate  is  computed  by  the  program. 
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Wind  Test  Set-up 

The  poles  were  tapered  beams  but  are  modeled  as  piece  wise  prismatic 
beam.  This,  of  course,  introduces  some  approximations  in  the  analysis.  It  is 
considered,  however,  not  to  be  significant.  The  poles  were  also  prestressed.  To 
adequately  model  this,  the  young  modulus  of  the  poles  were  selected  such  as  to 
yield  the  1  inch  displacement  at  the  top.  Each  pole  was  pre-loaded  prior  to  testing 
to  obtain  the  force  that  would  produce  1  inch  deflection  at  the  top.  Ignoring  shear 
deformation,  the  required  EI  was  then  computed.  This  also  accounted  for 
possible  movement  in  the  foundation.  The  foundation  for  the  poles  required  an 
embedment  length  of  five  to  six  feet.  Back  fill  earthwork  was  used  around  the 
pole  footing  instead  of  concrete  for  ease  of  removal.  The  properties  of  the  pole 
were  chosen  such  that  the  EI  value  of  pole  1  is  3,550,000  kip-in2  and  the  EI  value 
for  pole  2  is  3,280,000  kip-in2. 
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Table  6-2    Section  Properties  for  Members  shown  Fig.  6.1 


Section 

Properties 

Young's 
Modulus  (ksi) 

Area 

(in2) 

Moment  of  Inertia  (in4) 

Poles 

4000 

93.4 

I2=888 
I3=888 
J=484 

4000 

93.4 

I2=820 
I3=820 
J=484 

Cables 
Mess.  &  Cat. 

23000 

0.079 

— 

Connectors 
1-Viinch  pipe 

29000 

0.799 

I2=0.310 
I3=0.310 
J=0.620 

strap 

29000 

0.5324 

I2=0.167 

I3=0.0153 

J=0.500 

Description  of  Wind  Load  Test 


Three  Wind  load  tests  systems  were  considered  in  Ref.  [31].  A  two-point 
connection  system,  refer  to  Fig.  4.7,  a  one-point  connection  system  where  both  the 
messenger  and  catenary  cables  are  both  connected  to  the  same  point  on  the  poles, 
and  a  two-point  connection  system  with  drastically  reduced  sag.  For  the  purpose 
of  this  dissertation  only  the  first  test  system  is  considered.  This  is  because  the  aim 
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is  to  verify  the  theories  presented  in  this  dissertation. 

Comparison  of  Results 

For  the  dynamic  tests,  the  author  has  chosen  to  concentrate  only  on  the 
variation  in  cable  tensions  in  relation  to  the  variation  in  wind  velocities.  The 
loading  pattern  adopted  is  described  next.  The  dead  loads  were  applied  first.  The 
signals  were  hung  directly  from  the  catenary  cable  using  disconnect  boxes  and 
hangers.  The  pipe  hanger  was  used  for  the  East  signal  which  lie  along  nodes  6,  14 
and  17  in  Fig.  6.1 .  This  was  the  only  signal  instrumented  for  lateral 
displacements.  The  strap  hanger  was  used  for  the  West  signal  along  nodes  7,  15 
and  18  in  Fig.  6.1.  Wind  load  was  applied  using  the  air  boats  at  three  stages:  a 
low  speed  range  between  33  to  46  mph,  a  medium  speed  range  of  47  mph  to  64 
mph  and  a  high  speed  range  of  73  mph  to  99  mph.  The  wind  load  generated  on 
the  signals  were  kept  similar  as  best  as  possible.  The  range  of  difference  of  wind 
speeds  applied  to  the  signals  vary  from  about  7  mph  to  a  maximum  of  1 1  mph. 
The  duration  of  each  test  was  1 3  sec.  The  average  wind  speed  for  the  low  wind 
speed  range  is  41.684  mph  (West  Signal)  and  38.198  mph  (East  Signal).  The 
average  wind  speed  for  the  medium  wind  speed  range  is  53.146  mph  (West 
Signal)  and  51.819  mph  (East  Signal).  The  average  wind  speed  for  the  high  wind 
speed  range  is  88.173  mph  (West  Signal)  and  104.178  mph  (East  Signal). 
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Test  Data 


The  initial  tensions  in  the  cables  prior  to  application  of  wind  loading  are 
843  lbs  in  the  catenary  cable  and  725  lbs  in  the  messenger  cable.  For  Low  wind 
speed  range  the  data  obtained  from  the  2  point  connection  test  set-up  is  shown  in 
Table  6-2. 


Table  6-3     Wind  Data  for  Low  Wind  Speed  Test 


Time 
(hr:min:sec) 

Catenary 

Cable 

Tension  (lb) 

Messenger 

Cable 
Tension  (lb) 

West  Air  boat 

Wind  Speed 

(mph) 

East  Air  boat 
Wind  Speed 

(mph) 

10:52:18 

741.4 

993.6 

38.5 

38.8 

10:52:19 

773.6 

975.4 

38.7 

40.7 

10:52:20 

750.2 

1014.3 

41.0 

39.1 

10:52:21 

770.6 

997.7 

44.1 

36.4 

10:52:22 

722.1 

1036.6 

42.1 

38.0 

10:52:23 

762.9 

1000.1 

39.6 

38.4 

10:52:24 

736.7 

1033.2 

41.6 

39.2 

10:52:25 

713.8 

1051.4 

43.3 

35.9 

10:52:26 

744.3 

1021.8 

43.2 

39.6 

10:52:27 

712.7 

1052.7 

43.0 

36.0 

10:52:28 

760.5 

1005.8 

41.9 

39.1 

10:52:29 

725.5 

1005.3 

42.6 

36.9 

10:52:30 

724.8 

1021.7 

42.4 

38.5 

Averages 

741.5 

1016.1 

41.7 

38.2 
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For  medium  wind  speed  range  the  data  obtained  from  the  2  point-  connection  test 
set-up  is  shown  in  Table  6-3. 


Table  6-4     Wind  Data  for  Medium  Wind  Speed  Test 


Time 
(hr:min:sec) 

Catenary 

Cable 

Tension  (lb) 

Messenger 

Cable 
Tension  (lb) 

West  Air  boat 

Wind  Speed 

(mph) 

East  Air  boat 

Wind  Speed 

(mph) 

10:52:49 

702.3 

1155.2 

51.2 

54.7 

10:52:50 

709.1 

1202.6 

51.4 

56.6 

10:52:51 

734.6 

1182.2 

52.2 

48.7 

10:52:52 

750.6 

1174.0 

52.3 

53.1 

10:52:53 

727.6 

1178.5 

49.7 

54.7 

10:52:54 

744.8 

1178.1 

54.2 

57.3 

10:52:55 

735.9 

1203.9 

54.3 

54.1 

10:52:56 

734.3 

1196,1 

52.1 

54.4 

10:52:57 

737.0 

1197.9 

55.0 

47.6 

10:52:58 

740.9 

1164.6 

52.4 

53.1 

10:52:59 

734.0 

1152.3 

52.3 

43.1 

10:53:00 

705.1 

1168.1 

55.7 

44.8 

10:53:01 

691.9 

1158.0 

56.9 

50.9 

10:53:02 

689.8 

1208.6 

54.3 

52.7 

Averages 

724.1 

1178.0 

53.1 

51.8 

152 


For  high  wind  speed  range  the  data  obtained  from  the  2  point-  connection  test  set- 
up is  shown  in  Table  6-4. 


Table  6-5     Wind  Data  for  High  Wind  Speed  Test 


Time 
(hr:min:sec) 

Catenary 

Cable 

Tension  (lb) 

Messenger 

Cable 
Tension  (lb) 

West  Air  boat 

Wind  Speed 

(mph) 

East  Air  boat 

Wind  Speed 

(mph) 

10:53:42 

662.6 

1540.4 

88.9 

107.5 

10:53:43 

674.8 

1602.1 

88.5 

107.5 

10:53:44 

660.8 

1591.4 

88.7 

101.4 

10:53:45 

656.2 

1586.5 

93.1 

102.4 

10:53:46 

663.6 

1595.6 

86.2 

103.7 

10:53:47 

632.5 

1577.5 

83.5 

102.1 

10:53:48 

661.0 

1651.6 

85.4 

102.6 

10:53:49 

670.2 

1591.6 

86.0 

107.3 

10:53:50 

653.2 

1625.4 

92.7 

103.4 

10:53:51 

661.5 

1636.9 

87.1 

105.4 

10:53:52 

637.7 

1635.2 

89.3 

108.7 

10:53:53 

671.5 

1626.1 

87.3 

101.0 

10:53:54 

639.8 

1620.7 

89.5 

101.2 

Averages 

657.3 

1606.2 

88.2 

104.2 

153 


Wind  Load  Modeling 


Comparisons  are  made  between  the  test  results  and  two  simulated  wind 
load  models.  It  is  important  to  note  that  though  it  could  have  been  easier  and 
maybe  more  accurate  to  input  the  entire  wind  history  from  testing  into  the 
computer  program  for  comparison  of  analysis,  the  writer  has  chosen  to  test  the 
wind  generation  routines  incorporated  into  the  program,  hence  only  simulated 
wind  history  is  used.  The  wind  loading  models  are  discussed  in  chapter  4.  The 
first  wind  model  is  based  on  the  averages  mentioned  in  the  preceding  section  held 
constant  for  the  duration  of  the  test.  The  second  wind  loading  is  also  based  on  the 
averages  but  with  a  fluctuating  component  about  the  mean  values.  The  average 
wind  speed  used  as  input  in  the  wind  analysis  was  the  average  of  the  averages  on 
the  East  and  West  air  boat  signal  heads.  These  values  are  shown  in  Table  6-6 
below. 


Table  6-6  Wind  Speeds  and  Resulting  Drag  Forces  Used  in  Dynamic  Analysis 


Wind 
Speed 

Average  Input 
Force  on  Signal  (lbs) 

Input  Wind  Speed 
(mph) 

Constant* 
Velocity 

Variable** 
Velocity 

Test 
East/West 

Constant* 
Velocity 

Variable** 
Velocity 

Low 

45.54 

40.19 

38.2/41.7 

39.9 

40.2*** 

Medium 

43.60 

46.43 

51.8/53.1 

52.5 

53.2*** 

High 

108.5 

97.88 

104.2/88.2 

88.2 

88.8*** 
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Table  6-7  Cable  Tensions  for  Wind  Load  Test  #  1 


Wind 
Speed 

Catenary  Cable  Tensions 
(lbs) 

Messenger  Cable  Tensions 
(lbs) 

Test 

Constant* 
Velocity 

Variable** 
Velocity 

Test 

Constant* 
Velocity 

Variable** 
Velocity 

Low 

741 

628.0 

752.6 

1016 

964.9 

670.0 

Medium 

724 

436.0 

614.2 

1180 

1181.0 

991.0 

DYNASS  Constant  Wind  Speed  Model  (See  chapter  4  for  details) 
DYNASS  Variable  Wind  Speed  Model  (See  chapter  4  for  details) 
Average  velocity  from  fluctuating  wind  speed  generated. 
High  wind  ranges  produced  non-conclusive  results  due  to  divergence  of 
analysis.  There  are  many  variables  in  dynamic  wind  load  computation  that  may 
account  for  the  differences  between  the  measured  and  computed  results.  Both 
constant  and  variable  speed  wind  models  under  predict  the  tensions  in  the  catenary 
cable.  The  percent  difference  for  low  wind  speed  is  15.25%  and  ~  0%  for  the 
constant  velocity  and  variable  velocity  models  respectively.  For  the  medium  wind 
speed  range  the  percent  difference  is  39.8%  and  15.17%.  The  discrepancy  in  the 
results  could  be  as  a  result  of  the  inability  to  more  adequately  model  the  loading 
procedure.  The  comparisons  for  the  messenger  cable  are  the  same  with  a  percent 
difference  of  ~  0%  and  34.1%  respectively  for  the  low  wind  range  and  ■  0%  and 
16.0%  respectively  for  the  medium  wind  range.  At  this  time  the  analysis  for  the 
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high  wind  speed  range  for  winds  of  about  100  mph  gives  a  non  converging  result. 
This  is  because  at  that  wind  speed  the  analysis  predicts  an  overturning  of  the  signal 
light  head.  This  effect  of  this  overturning  was  observed  to  cause  drastic  variations 
of  the  drag  and  lift  forces  applied  at  the  signal  heads.  The  range  of  the  drag  force 
was  300  lbs  to  0  lbs.    Mathematical  recovery  after  such  an  overturning  is  being 
investigated  and  is  recommended  for  further  research.  Other  factors  that  could 
account  for  the  differences  in  the  results  include  discrepancies  in  initial  tension 
readings  between  tests  and  an  inadequate  modeling  of  wind  pressure  sustained  at 
the  signal  head.  The  low  wind  range,  medium  wind  range  and  high  wind  range 
were  tested  continuously  without  stopping.  The  computer  model,  however,  only 
accounts  for  the  change  in  initial  conditions  between  tests  by  changing  the  initial 
tension  in  the  cables.  The  dynamic  effects  in  the  transition  were  lost  and  were  not 
modeled.  Both  wind  models  developed  attempts  to  take  account  of  the  rotation  of 
signal  head,  in  the  wind  stream,  by  the  use  of  drag  and  lift  coefficients.  Inherent  in 
the  use  of  these  coefficients  is  the  assumption  that  there  is  an  uncoupling  of  the 
along  the  wind  response  of  the  lights  and  the  cross  wind  response.  And  a  further 
assumption  that  the  along  wind  response  controls.  Both  these  assumptions  are 
presently  unverifiable.  The  dynamic  solution  approach,  barring  the  assumptions  in 
loading,  is  seen  to  give  a  reasonable  lower  bound  solution  to  the  cable  -  frame 
problem. 


CHAPTER  7 
CONCLUSION 


Summary 

The  purpose  of  this  dissertation  was  to  present  an  alternative  analysis 
approach  to  the  force  density  method  used  in  the  ATLAS  program  developed  in 
the  Civil  Engineering  Department  at  University  of  Florida.  This  new  approach 
was  to  increase  flexibility  of  analysis,  solve  the  static  analysis  of  a  cable  -  frame 
assemble  and  be  able  to  model  true  dynamic  loading  on  the  traffic  signal  system. 
The  new  approach  espoused  uses  the  principles  of  dynamic  analysis.  This 
solution  approach  uses  finite  element  techniques  to  model  the  traffic  signal  system 
as  made  up  of  frame  elements,  cable  elements,  connector  elements  and  truss 
elements.  The  elastic  and  inertia  properties  of  each  element  in  the  structure  is 
assembled  to  form  the  structures'  stiffness,  mass  and  damping  properties.  The 
nodal  equations  of  motion  for  the  structure  are  then  set-up  in  terms  of  the 
structures'  assembled  global  properties.  These  nodal  equations  of  motion  are 
solved  using  the  Newmark's  Beta  numerical  integration  and  Newton  -  Raphson 
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iteration  scheme.  The  analysis  technique  presented  was  compared  to  test  results 

from  tests  performed  at  the  University  of  Florida  and  to  other  analysis  methods. 

Review  of  Analysis  Method  Developed 

A  2-  node  cable  element  was  explicitly  derived  for  the  finite  element 
analysis  of  the  cable  -  frame  system.  Built  into  the  formulation  of  this  cable 
element  are  the  assumptions  that  the  element  can  not  resist  bending  or 
compression.  The  cable  element  is  hence  rendered  inactive  when  subjected  to  any 
system  of  compressive  forces.  The  inertia  properties  of  the  cable  element, 
however  remain  active  even  under  compressive  forces  hence  restoring  the  system 
and  keeping  the  system  in  equilibrium.  The  cable  element  properties  are  derived 
directly  in  terms  of  the  deformed  global  coordinates  of  the  nodes  as  opposed  to  the 
classical  local  coordinate  system  with  a  global  -  local  transformation.  The  finite 
element  method  presented  is  based  upon  the  total  Lagrangian  version  of  the 
incremental  equation  of  motion.  Solution  to  the  equation  of  motion  is  obtained 
using  an  improved  Newmark's  Beta  method.  The  improvement  is  in  the  form  of 
the  introduction  of  numerical  damping  factor.  The  numerical  damping  uses  novel 
manipulations  of  the  inertia  forces  to  stabilize  an  otherwise  statically  precarious 
system.  This  is  used  in  addition  to  the  conventional  dissipative  damping  matrix. 
The  main  difference  between  both  damping  processes  is  that  using  numerical 
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damping  retains  the  "true"  dynamic  response  of  system,  while  arbitrarily 

increasing  the  conventional  damping  can  result  in  over  damping  that  can  mask  the 

"true"  response  of  the  system.  The  improved  Newmark's  Beta  method  reduces 

the  nonlinear  set  of  equations  to  corresponding  set  of  linear  equations  that  can 

then  be  solved  using  Newton  -  Raphson.  A  modified  Newton  -  Raphson  was 

presented.  In  addition  to  the  regular  modifications  to  the  method  to  improve 

convergence,  the  author  makes  use  of  a  technique  of  averaging  the  tangent 

stiffness  of  the  structure  in  successive  Newton  -  Raphson  iterations.  This  is 

observed  to  curb  the  divergent  tendencies  of  the  algebraic  equations  in  an  unstable 

cable  analysis. 

The  solution  method  is  completed  with  the  definition  of  two  dynamic  wind 

load  models.  The  first  uses  a  constant  average  wind  speed  but  accounts  for  the 

rotation  of  traffic  signal  head  by  using  lift  and  drag  coefficients.  The  second 

model  uses  the  average  wind  speed  plus  a  gust  component.  This  model  also 

accounts  for  the  rotation  of  the  traffic  signal  head.  The  gust  component  of  the 

second  wind  model  is  generated  using  auto  regressive  numerical  filters.  Both 

wind  models  generate  along  wind  forces.  Cross  wind  effects  are  ignored.  The 

solution  method  presented  adequately  handles  a  nonlinear  analysis  with  the 

sources  of  nonlinearity  being  large  displacement  in  cable  and  connector  beam 

elements  and  nonlinear  wind  load. 
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Verification  of  Analysis  Method 

A  program  was  developed  based  on  the  solution  method  presented.  This 
program,  DYNASS,  was  applied  to  the  traffic  signal  system.  The  computed 
results  from  this  program  was  compared  to  test  results  of  experiments  done  at  the 
University  of  Florida.    Good  correlation  was  observed  for  the  static  load  case 
results.  For  the  dynamic  load  case,  computed  results  are  considered  comparable  to 
the  measured  result  but  with  appreciable  difference  due  assumptions  in  loading. 

Advantages  and  Disadvantages  of  Solution  Method 

The  solution  method  was  verified  for  the  traffic  signal  system.  A  major 
advantage  of  this  solution  technique,  however,  is  it's  flexibility.  The  method  of 
analysis  presented  can  be  used  for  any  system  of  cables  and  frames.  This  retains 
the  power  of  the  finite  element  technique  of  being  valid  for  any  structure. 
Another  advantage  of  this  solution  approach  is  the  additional  information 
obtained  from  a  typical  analysis.  Though  this  solution  approach  can  be  used  for 
static  analysis,  a  lot  more  information  is  gathered  from  the  analysis  than  just  the 
static  solution.  Since  the  solution  method  uses  dynamic  theories,  an  advantage  of 
this  solution  approach  is  the  capability  to  model  dynamic  loads. 

A  major  disadvantage  is  the  time  it  takes  for  a  typical  solution.  A  typical 
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solution  can  be  obtained  in  3  sec  analysis  time.  However,  due  to  the  nonlinearity 

introduced  by  the  cables  and  connector  elements,  the  solution  time  step  required 

range  between  0.001  and  0.0001  sec.  This  implies  that  a  typical  solution  may 

require  30000  iterations.  Each  iteration  takes  about  0.01  sec.  in  real  time.  This 

means  a  typical  nonlinear  analysis  could  take  as  long  as  50  mins. 


Recommendations  for  Further  Research 


There  are  many  ways  the  solution  method  presented  could  be  improved. 
One  way  is  the  implementation  of  a  variable  time  step.  This  would  make  the 
analysis  faster.  A  variable  time  step  allows  for  the  use  of  a  large  time  step  in 
portions  of  the  analysis  not  sensitive  to  the  time  step  chosen.  It  would,  however, 
allow  for  a  reduction  in  the  time  step  if  divergence  of  solution  is  observed. 
Another  area  recommended  for  further  research  is  the  wind  load  computation. 
The  effects  of  cross  wind  in  structural  building  design  is  increasingly  gaining 
attention.  These  cross  wind  effects  are  presently  ignored  in  the  analysis  as  is  the 
norm.  It  is  the  author's  opinion  that  ignoring  this  phenomenon  under  states  the 
forces  acting  on  the  traffic  signal  system.  A  final  recommendation  is  an 
implementation  of  other  load  types.  Since  this  solution  method  can  be  applied  to 
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mooring  lines,  it  is  recommended  that  an  additional  load  type  be  considered,  i.e. 
hydrostatic  pressure.  Finally,  extensions  to  the  theory  could  be  made  to  handle 
base  excitation  as  would  be  experienced  in  an  earthquake. 


APPENDIX  A 
CABLE  ELEMENT  DERIVATION 
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Nonlinear  2-node  cable  element. 


6 

A 
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♦  1-1 


+  2-1 


+  3=1 


u=: 


*5S 


♦6°; 


The  above  are  the  shape  functions  for  the  2  node  element  shown. 
The  equation  for  generic  displacement  of  any  point  along  the  element 
with  respect  to  the  nodal  displacements  can  be  given  as  shown  below. 


U=i|)1-u1  +  +4^2 
V=(j»2v1  •  <t>5v2 

W=(j>3-W  1  r +  g-W2 

Note  that  the  notations  and  coordinate  systems  adopted  in  this  appendix  is  as 

described  in  the  chapter  3. 

s  is  the  arc  length  of  the  cable  from  node  1  to  any  generic  point  along  the  cable. 

L  is  the  initial  arc  length  of  the  cable. 

{xj  yj  Zj  x2  y2  Z2  }  are  me  n°dal  coordinates. 

{u,  v,  Wj  u2  v2  w2  }  are  the  nodal  displacement  vector. 
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In  matrix  form,  N  is  the  above  linear  displacement  field.  N'  is  the 
first  differential  of  the  shape  function  with  respect  to  s.  i.e.  the  arc 
length. 


To  obtain  the  equivalent  nodal  loads  due  to  the  weight  of  the  cable  we  have 

•L 


9= 


0       1 


WN'I 

0 
0 


The  h  vector  is  a  coordinate  dependent  matrix.  Since  the 
load  is  acting  in  the  negative  z  direction  to  our  assumed 
coordinate  system.  The  h  vector  is  as  shown  below. 


o 

0  |i 
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o 
o 

(1-8) 

L 

0 

0 


=NTh 


WN'hdsa 


■L 


■L 


W 


Consider  the  initial  tension  in  the  cable  as  a  distributed  load. 


Straight  cable 


fCf 


T0N'hds 


166 


o      1 


0       1 


1 

•  0 
0 


■NTh 


(1-8) 

L 

0 
0 
s 

L 
0 
0 


«NT-h 


Note  that  because  of  the  assumed  tension 
distribution  along  the  member  this  f0  is  in  the 
element  local  axis. 


fo- 


4T, 


1  + 


2s 


[(L-s)j 

L 

T'. 

0 

0 

0 
s 

L 

ds» 

0 
3T» 

0 

0 

0 

0 

In  structure's  global  axis: 


dx=X2     x.| 
dy=*y  2    y  1 


dz=z< 


1 


xl-Vdx2  t  dy2  t  dz2 


dx 

xl 

dy 

xl 

dz 

'o-3'To- 

xl 
dx 

xl 

dy 

xl 

dz 

xl 

The  mass  matrix  is  computed  from: 


W       T 

— -N-Nds 
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Note  that  W  is  the  total  weight  of  the  element. 

L  is  the  length  along  the  undeformed  arc  of  the  element. 

g  is  the  acceleration  due  to  gravity.  This  is  used  to  convert  the  weight  into  mass. 


0       1 


0      1 

0         0 

s        0 
L 

0        • 
L 


1°       0 
L 


0    0 


0       1  0       0         0  =NTN 

L  L 


0  0       1     S    0    0    S 


(Us)' 


(-Us). 


NTN= 


(Us)- 


(-Us)' 
L2 

0 


(Us)' 
L2 

0 


(-Us)- 


(  Us)- 


(L  +  s)- 


(L+s)- 


0 


s2 
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m= 


.  W 

'eg 


2   0   0   10   0 

0   2   0   0   10 

0   0   2   0   0   1 

10   0   2   0   0 

0   10   0   2   0 

0   0   10   0   2 

This  the  mass  matrix  in  global  coordinates. 


From  chapter  3  the  elastic  stiffness  of  the  above  element  can  be  obtained  from 
the  following.  Eqtns  3-10  and  3-22  give  the  expression  for  strain. 


E.(tt  At)»6_{t)  -hAe, 


dnT(N')-(N')-(dn)-1 


Ass=dnT-(N')T(N')-(Adn)  +  i 


AdnT(N')T(N')-(Adn 


(3-16) 


Change  in  strain. 


AE 


n')t(n' 


d  nT-  (n')T  (n')-  (d  n)    1  ds       (3-25)  KE  is  the  linear  part 

of  the  cable  stiffness. 


N' 


'AE- 


dnl 


AdnT(N')T(N 


dnMN')T-(N%i 


AdnT(N')T 


ds  (3-26) 


KL  is  the  non-linear  part  of  the  cable  stiffness  formulation. 

Substituting  the  shape  functions  into  the  above  expressions  for  KE  we  have  the 

following. 
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Compute  the  term:  d  n ' •  (N  j  •  (N  j  •  (d  n)  -  1 


0      0 

1    0 

L 


0     0 


0 


1 
L 

0 

0 

1 

L 


- 1                1 

—     0     0         0    0 

L                  L 

1                1 
0    —    00-0 

L                L 

■(n 

)'(n 

1              1 
0      0           0    0 

L              L 

1      0     o      1    0     0 


1  1 

0—0     0    —    0 


1  1 

0     0  0     0    — 


1  1 

—    0000 


1  1 

0    —    0     0    —    0 


1  1 

0      0-00 


.(n')t-(n) 
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x  1    Y1    z1    x2   V2   z2 


1  1 

0      0  0      0 


1  1 

0-0      0-0 


1  1 

0      0  0      0 


1  1 

0      0  0      0 


1  1 

0  0      0  0 


1  1 

0      0  0      0 


x1 

Y1 

z1 

x2 

y2 

z2 


■dnT(N')T(N')-(dn) 


2x1x2-2y1y2  +  x12  +  y12  +  z12-  2zrz2^x22  ^  y22    z22 


■dn'-lNj  -INJ-  dn 


2-xrx?-  2y1y2  +  x12 ry/tz/    2zrz2    x22    y22  +  z22)  T 


T  /    • 


2xrx2    2y1y2  +  x12  +  y12  +  z12    2zrz2  ,  x22  +  y22  t  z22 


N')  -(N')ds 


Multiply  the  above  equation  with  each  term  in  the  NHSf  matrix  and  factor  expression 
we  have  the  following: 
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2xrx2    2yry2-x12fy12  +  z12     2zrz2 i  x22 1  y22  +  z22 


ds 


Note  that  x,  y,  and  z  coodinates  are  for  any  time  (t). 


■ir 


2LJ 


X!    x2)2+  (yr  y2)2  + 


K22aJ4{(xi-x2)2+(yi-y2)2+ 


•33" 


'44" 


K55= 


2L3 
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14" 


2LJ 


K25= 


2LJ 


(xi  -  x2)2^  (vi-y2)2+- 


(x1-x2)2+(y1-y2)5 


K36sr-J3{(xi-x2)2^(yry2: 


'41' 


2-L 


2L° 


x1-x2)2^(y1-y2)2  + 


K52=^{(xi-x2)2+(yr-y2)5 


'63- 


2L° 


2LJ 


x1-x2)2+(y1-y2)2 


Zl-z/     L2 


z1     z2^     L2 


Zl-z2l2     L2 


-z2l2     L2 


z^z/     L2 


zrz2^     L2 


Z1     z2)2     L2 


z1     ^     L2 


Z1     z2)2     L2 


Z1     z2,2     L2 


Z1-Z2^2     L2 


Z1-Z2)2     L2 
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KE"AE- 


11 

0        0      K14      0 

0 

0 

K  oo      0        0         25 

0 

0 

0      K33      0        0 

K36 

41 

0        0      K44      0 

0 

0 

K52      °        °      K55 

0 

0        0      K63      0        OK 


66 


KE  is  referred  t  as  the  initial  stress  matrix. 

It  is  computed  at  the  beginning  of  each  time 
step.  It  is  remains  constant  within  the  time 
step. 


KE« 


AE 
2L3' 


i-x2  +  (y  1  - y 2)  +■  zi  z2 


2      ,2 


10  0  10  0 
0  10  0  10 
0  0  10  0  1 
10  0  10  0 
0  10  0  1  0 
0     0      10     0     1 


Note  that  the  tension  in  the  cable  at  time  t  is 


h'f{^-*2)2-^--y2)2-^--^  l2 


It  is  more  convenient  to  cast  this  tension  expression  in  terms  of  incremental 
strain  in  order  to  account  for  pre-tension.  In  which  case  this  expression  becomes: 

T  t»T(   M  t  AE  At  Change  in  strain  can  be  computed  using  Eqn.  3-16. 


The  initial  tension  is  taken  from  either  a  given  value  or  computed  from  the  slack 
configuration. 


AT-EA- 


dn-lNj    N 


AdnT-(N')TN' 


•(Adn) 


Compute  the  RE  as  shown  in  Eqn.  3-27. 
•L 


AE 


(n')t(n 


dnT(N')T(N')-(dn' 


dnds 


0 
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The  coresponding  nodal  force  vector  for  the  stress  at  any  time  t  is  as  follows. 
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0        0      K 
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0      K 


0        0 
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.A-E 
2-L3 


11x  1  14x2 

K22-yi  K25y2 

K33'z1  '  K36'z2 

^ 41  "x  1  '  ^44x2 

K52-yi  K55y2 

K63'z1  K66'z2 


■Ri 


RE= 


A-E 
2L3 


x2)2Wyi-y2)2 


i,\2-\.2 


x1 

x2 

V1 

V2 

zr 

z2 

x2 

x1 

y2 

yi 

z2 

zi 

Or  alternatively, 


RE=: 


AE      A 

'-*     2-L2 


x1 

x2 

u 

V2 

z1 

z2 

x2" 

x1 

V2 

V1 

Z2 

z1 

Similarly  we  can  compute  the  KL  matrix  as  follows. 

•L 


K,   -AE- 


dnT-(N')T(N')  +  AdnT(N')T(N 


1  [.,   W„-\T 


dn'-lNj  •iNj  +  ^-|Adn' 


ds 


0 


The  term     d nT- (n)  •  (n)      is  equal  to  the  following. 
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x1    V1    z1    x2   V2   z2 
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=dnT-t 


i-x2    yi-y2    zi-z2   -  xi-x2   ■  yi-y2j    zi-z2 


L2  L2 


■dn-  N    •  N 


The  term       Ad  /-(n)  ■  (n)     is  nonlinear  and  dependent  on  Adn.  Note  that  Adn  is 

unknown  at  the  start  of  an  analysis.  Therefore,  we 
replace  the  above  expression  for  KL  with  KA 

as  an  initial  approximation. 
■L 


KA=A-E. 


<y-u 


dn'-lNj 


ds 


(3-28) 


This  expression  is  used  at  the  beginning  of  the  time  step  analysis.  Within  each  time 
step,  i.e.  during  the  Newton-Raphson  iteration,  the  true  expression  will  be  used. 
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Compute  KA: 


Let  Q* 


dnM 


dn'lNj 


x1     x2 


Y1-Y2 


z1     z2 


I/ 
x1  +  x2 


1/ 


I/ 
z1  +  z2 


xi  x2j  (vi-y2j  (zi"z2j    (xi-x2j    (y  1  y 2;  izi-z2 


=Q 


For  convenience,  let  the  above  expressions  be  replaced  as  follows. 


1~x2       yi-*2       z1-z2     ix1-x2     -yi-V2         z1z2     _1 


L2  L2  L2  L2  L2  L2  L2 


(  dx    dy    dz  dx  dy  dz) 


u1     u2)     (v1"v2)     (w1"w2)      (u1-u2)      (v1-v2)      (w1"w2 


L2  L2  L2  L2  L2  L2 


•(  du     dv    dw  du  dv  dw) 


The  matrix  KA  is  therefore  as  given  below: 
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(  dy)-(dx)     (-dz)-(dx) 
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(  dx) 
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2 


(  dx)-(-dy)   (-dx)-(-dz) 

2 


dz)' 


dx)(  dz)    (  dy)(  dz)       (-dz)' 


The  above  expression  for  KA  is  used  at  the  beginning  of  the  time  step. 
Within  the  time  step  the  following  derived  KL  is  used. 
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T/V 
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The  term  in  the  integral  therefore  becomes  the  following. 
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See  the  next  page  for  the  KL  matrix. 


At  the  beginning  of  the  time  step  analysis:       K=K  E  -t-  K  A 
For  the  rest  of  the  analysis:  K=K  EiKL 
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APPENDIX  B 
USER'S  MANUAL  FOR  DYNASS 
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DYNASS  USER'S  GUIDE 


INTRODUCTION 


DYNASS  is  a  dynamic  analysis  program,  it  does  a  nonlinear  analysis  of  a  cable  structure 
including  but  not  limited  to  the  traffic  signal  systems.  The  method  of  analysis  used  is  the 
implicit  numerical  integration  of  a  piece- wise  linear  equation  of  motion.  The  default 
implicit  numerical  scheme  is  the  Newmark's  Beta  method.  The  program  can,  however, 
simulate  other  numeric  integration  schemes  e.g.  central  difference,  backward  difference 
and  Euler  integration  method.  For  the  Newmark's  Beta  method,  the  program  provides  a 
numerical  damping  feature  called  the  alpha  modification.  This  enables  a  faster  die-down 
of  the  structure's  vibration  while  still  maintaining  the  true  dynamic  properties  of  the 
system.  Table  B-l  shows  the  factors  to  used  for  the  different  methods. 


Table  B-l    Integration  < 

constants  foi 

use  in  DYNASS 
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Method  of  Analysis 

a 

P 

Y 

Newmark's  Beta  method 

Constant  acceleration 

Linear  acceleration 

Variable  acceleration 

0 
0 
0 

0.25 

0.167 

0.125 

0.5 
0.5 
0.5 

Alpha  Modification 
Option  1 
Option  2 
Option  3* 

0.1 
-0.1 
-0.1 

0.3025 

0.5 
0.3025 

0.6 
0.6 
0.6 

Central  Difference  method 

0 

0.5 

0.5 

Backward  Difference  method 

0 

1 

0.5 

Euler  Integration  method 

0 

0 

0.5 

This  option  is  suggested  for  traffic  signal  system  analysis. 
To  run  DYNASS  it  is  required  that  the  user  prepare  an  input  file,  the  following  section 
discusses  the  terminology  in  use  for  an  adequate  input  to  the  program  and  the  general 
form  of  the  input  file.  The  program  requires  that  the  file  serving  as  input  to  the  program 
be  called  INPUT.    The  results  from  the  program  is  written  to  a  file  called  OUTPUT.  In- 
addition,  two  plot  files  are  generated  by  the  program.  The  first  is  the  time  -  displacement 
history  for  previously  chosen  nodal  points.  The  second  is  time-velocity  listing  for  the 
wind  loading  generated  by  the  program.  The  first  plot  file  is  called  PLOT  and  the  second 
is  SVEL.  These  can  be  plotted  by  regular  plotting  routine.  A  plotting  routine  is  currently 
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not  included  in  the  latest  version  of  DYNASS.  The  terminology  used  in  the  input  file 
refers  to  the  traffic  signal  cable  structure  as  shown  in  Fig.  B.l . 


Catenary  Cable 


Preslressed  concrete  pole 


Fig.  B.l  Definition  of  terms  used  in  the  INPUT  file 
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NOMENCLATURE 


CATENARY-  This  indicates  the  primary  dead  load  supporting  cable.  The  word  catenary 
is  derived  from  the  root  word  "catenate"  which  implies  to  link  or  tie  strings  together.  The 
initial  configuration  for  this  cable  is  assumed  as  a  parabolic  shape  when  there  is  no  pre- 
tension. 

MESSENGER  -  This  indicates  the  lower  cable.  It's  name  was  derived  from  the  electrical 
messenger  wires  that  it  supports. 

BEAMS  -  This  refers  to  the  support  strain  poles  at  each  end  of  the  span.  In  the  case  of  the 
traffic  signal  system,  these  are  prestressed  concrete  poles. 

CONNECTOR  -  This  refers  to  the  element  connecting  the  catenary  cable,  the  messenger 
cable  and,  often  times,  a  traffic  signal  head  or  auxiliary  sign  or  some  other  cladding  or 
load  attachment. 

INITIAL  SAG  -  This  is  the  sag  in  the  cable  due  only  to  cable  self  weight  and  tension.  Sag 
due  to  signal  lights  is  not  included.  This  is  used  as  a  starting  position  for  the  solution 
procedure. 

WIND  -  This  refers  to  wind  loads  acting  on  the  signal  light,  accessories,  cables,  and  poles. 
This  load  is  specified  by  providing  an  average  wind  speed  and  the  direction  of  wind.  To 
model  gust  the  length  of  wind  record  and  the  sampling  interval  for  the  wind  data  must  be 
provided. 


LOADS  -  This  refers  to  the  additional  gravity  and  lateral  loads  acting  on  the  signal  light 
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and  is  used  in  combination  with  the  wind  loads. 


SIGNAL  -  This  is  used  to  complete  the  wind  load  definition.  A  detailed  description  of  the 
signal  or  sign  gives  the  area  impacted  by  the  wind  pressure  and  the  gravity  load  of  the 
signal  light. 

NONLINEAR  -  The  nonlinear  parameters  pertaining  to  the  cable  geometric  properties  are 
the  tension  in  the  cable  and  the  original  sag  in  the  catenary  cable. 

DYNAMICS  -  The  die-down  tolerance  and  the  integration  constants  given  in  Table  B-l 
are  part  of  the  dynamic  parameters  used  by  the  program.  Other  parameters  specified  as 
dynamic  parameters  are  discussed  later  in  this  section. 

INPUT  FILE  SYNTAX 

The  following  is  syntax  rules  that  should  observed  when  attempting  to  write  the  input  for 
the  DYNASS  program.  The  input  file  is  free-format  hence  only  blank  spaces  or  commas 
are  needed  to  separate  data.  The  following  are  the  conventions  used  for  input: 

1 .         A  "C"  in  column  1  of  any  line  will  cause  the  line  to  be  echoed  as  a  comment  on 
the  console. 


A  back  slash  "\"  at  the  end  of  information  on  the  first  line  will  allow  the  next  line 
to  be  interpreted  as  a  continuation  of  the  first  line:  therefore,  a  160  character 
record  is  possible. 
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3.  A  colon  ":"  indicates  the  end  of  information  on  a  line.  Any  information  entered  to 
the  right  of  the  colon  is  ignored  by  the  program.  This  space  can  be  used  to  provide 
additional  comments  within  the  input  file. 

4.  If  a  blank  identifier  is  specified,  the  data  string  is  assumed  to  be  the  first  data 
string  of  the  record. 

5.  If  fewer  data  exist  than  are  specified,  the  values  returned  will  be  either  a  zero  or 
blank  according  to  the  type  of  data  expected. 

6.  Real  numbers  do  not  require  decimal  points:  E  formats  with  +  or  -  exponents  are 
accepted. 

7.  Simple  arithmetic  statements  may  be  used  within  the  input.  The  functions  that  can 
be  used  are  +,-,*,/•  The  order  of  evaluation  is  sequential,  not  hierarchical  as  in  the 
FORTRAN  language. 

INPUT  DATA 


The  following  is  a  guide  to  the  user  for  preparing  an  INPUT  file.  The  data,  as  shown 
below,  is  entered  in  sections.  Each  section  has  a  designation  and  is  separated  from  the 
others  by  a  colon  ":".  The  first  two  lines  have  no  apparent  section  designation  but  is 
referred  to  as  the  control  section.  The  first  line  is  the  80  character  title  line.  The  second 
line  contains  the  general  information  parameters.  The  positioning  of  the  first  two  lines  is 
critical  to  the  running  of  the  program  and  hence  can  not  be  omitted  or  repositioned.  The 
order  of  the  rest  of  the  sections  is  arbitrary.  It  is  not  necessary  to  separate  the  first  two 
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control  lines  from  the  next  section  with  the  colon  mentioned  above. 

The  information  necessary  in  the  second  line  control  section  are 

1 .  The  number  of  nodes  in  the  structure.  The  syntax  for  this  is  nodes=  #.  Where  the  # 
sign  is  a  numeric  value.  Note  that  there  should  not  be  any  white  space  between  the 
letter  's'  and  '='. 

2.  The  number  of  elements  types  in  the  structure.  DYNASS  is  capable  of  handling 
only 

beam,  truss,  cable  and  connector  elements.  The  syntax  for  this  is  elem_no=  #.  Or 
simply,  No=  #. 

3.  The  number  of  cables  in  the  structure  must  be  specified  next.  The  syntax  for  this 
is  cable=  #.  If  no  cable  present.  This  could  be  omitted  or  set  to  zero. 

4.  The  final  control  section  parameter  is  the  istat.  This  is  used  to  identify  the  kind  of 
analysis  to  run.  The  program  is  capable  of  three  types  of  analysis  depending  on 
the  wind  load  model  selected. 

A)  Wind  load  applied  as  a  lateral  load  not  altered  in  the  time  step 
solution.  For  this  istat=0.  The  computation  of  the  lateral  load  is 
left  to  the  user.  The  lateral  load  must  be  entered  in  the  section 
designated  as  loads. 

B)  Variable  area  wind  model.  Since  the  wind  is  impacting  a  free- 
swinging  traffic  signal.  The  rotation  of  the  traffic  signal  in  the 
along-wind  direction  is  taken  into  consideration  in  the  computation 
of  the  wind  force.  The  non-fluctuating  average  wind  speed  based 
on  the  ASCE  7,  3  sec  wind  gust  criteria  is  used  for  the  computation 
of  wind  pressure.  For  this  istat=l. 

C)  Variable  velocity  wind  model.  The  fluctuating  part  of  wind 
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velocity  is  generated  in  the  program  using  auto-regressive 
numerical  filters.  The  fluctuating  portion  of  the  wind  velocity  is 
added  to  the  average  wind  velocity  similar  to  that  in  paragraph  B 
above  to  compute  the  required  wind  pressure.  See  chapter  4  for 
details  of  the  different  wind  models  and  the  equations  used  in  wind 
load  computation.  For  this  istat=-l. 


Mx      Fx 


Fig  B.2  Coordinate  system  used  in  DYNASS 


Using  the  above  coordinate  system.  A  typical  input  could  be  generated  for  a  traffic  signal 
system  using  the  following  sequence  of  section  definitions.  It  is  important  to  note  that 
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the  default  plane  for  a  structure  in  DYNASS  is  the  x-z  plane.  For  a  3-D  structure,  the  z 
direction  as  the  elevation  of  the  poles  must  be  maintained.  For  a  structure  defined  in  the 
x-  y  plane  an  error  statement  of  not  finding  the  structure's  plane  will  be  printed  and  the 
program  terminated. 

NON-LINEAR  PARAMETERS 
This  section  describes  the  nonlinear  cable  elements.  The  syntax  rules  are  as  follows: 

NONLINEAR  -  (Section  Designation) 

NND   ND=N„N2,N3,...,NNND    SG=S1     TEN=T1     W=W1     P=P1 

NND    =  Number  of  nodes  in  the  cable  being  specified.  When  cables  are 

connected  to  poles,  this  number  includes  the  nodes  at  the  cable- 
pole  supports.  This  number  must  be  greater  or  equal  to  3. 
The  more  elements  per  cable  used,  the  better  the 
accuracy  of  the  solution. 

N,        =  First  node  in  the  cable.  Preferably  the  left  cable-pole  node. 

N2        =  Second  node  in  the  cable,  moving  from  left  to  right. 

N3        =  Third  node  in  the  cable  and  so  on 

Nnnd      =  Last  node  in  the  cable.  Preferably  the  right  cable-pole  node. 

S 1         =  The  original  sag  in  the  catenary  cable. 

Tl        =  The  pretension  in  the  catenary  cable.  This  is  optional.  If  none  is 

specified,  pre-tensioning  is  done  due  to  self  weight  of  the  cable 
only. 
If  one  is  specified.  The  value  must  be  specified  in  kips. 


Wl 


PI 
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weight  of  the  cable  (lbs/ft).  This  must  be  specified.  Note  the  units. 
Weight  must  be  specified  in  lbs/ft. 

This  is  the  type  of  cable.  P=  0  indicates  a  catenary  cable.  P=  1 
indicates  a  messenger  cable.  This  also  must  be  specified.  This 
Information  is  used  for  accurate  printing  of  results. 


NOTE:  The  prefix  ND=  SG=,  TEN=  W=  and  P=  must  precede  the  specification  of  the 
numerical  values  mentioned  above.  This  command  line  is  specified  for  each  cable  in  the 
structure.  Both  messenger  and  catenary  cables  are  specified  separately  if  both  present  in 
the  structure. 

This  section  MUST  end  with  a  blank  line. 

DYNAMIC  PARAMETERS 

This  section  describes  the  dynamic  parameters  necessary  for  analysis. 

DYNAMICS  -  (Section  Designation) 

TOL     A=  ALPHA    B=BETA    G=  GAMMA     P=P1    N=N1 


TOL     =  This  is  the  die-down  criteria.  This  must  always  be  less  than  or 

equal  to  0.001 .  The  smaller  the  value  the  better  the  equivalent 
static  solution.  The  only  caveat  is  that  the  smaller  the  tolerance 
value  the  longer  the  program  runs  to  achieve  die-down. 


190 

ALPHA  =        This  is  the  alpha  term  selected  from  Table  B- 1 . 

BETA  =  This  is  the  beta  term  selected  from  Table  B-l . 

GAMMA=      This  is  the  gamma  term  selected  from  Table  B-l . 

PI        =  In  order  to  be  sure  of  die-down,  the  structure's  response  of  a 

critical  node  must  be  plotted.  The  node  to  be  plotted  is  PI. 
Nl        =  The  displacement  degree  of  freedom  (dof)  for  the  node  to  be 

plotted  must  also  be  specified.  X  dof=l.  Y  dof=2.  Z  dof=3. 
This  section  MUST  end  with  a  blank  line. 

JOINT  COORDINATE  DATA 

This  section  defines  the  nodal  coordinates  in  a  positive  Cartesian  coordinate  system.  The 
following  lines  will  define  the  coordinates; 

COORDINATES  -  (Section  Designation) 

ND  X=X,  Y=Y,  Z=Z,  or  C=C, 


ND       =  Node  number  for  the  given  coordinates, 

X,         =  X  coordinate  for  the  node  ND, 

Y,         =  Y  coordinate  for  node  ND, 

Z,         =  Z  coordinate  for  node  ND, 

C,         =  Cable  number  that  the  node  belongs. 

NOTE:  Any  node  that  is  not  specified  will  default  to  the  origin. 

This  section  MUST  end  with  a  blank  line. 
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BOUNDARY  CONDITION  DATA 

This  section  defines  the  active  degrees  freedom  for  the  structure  to  be  analyzed.  All  the 
degrees  of  freedom  must  be  specified  as  fixed  or  released,  according  to  the  boundary 
conditions  and  the  types  of  elements  connected  to  the  node.  The  following  lines  of  input 
describe  the  state  of  the  DOF. 


BOUNDARY  -  (Section  Designation) 
NF,NL,NI  DOF=F„F2,F3,F4,F5,F6 

NF       =  First  node  in  a  generation  sequence  for  which  the  degree  of 

freedom  is  specified. 
NL       =  Last  node  in  a  generation  sequence  for  which  the  degree  of 

freedom  is  specified. 
NI        =  Increment  for  generating  node  numbers  between  NF  and  NL  for 

which  the  DOF  specification  is  used.  NL  and  NI  can  be  left  blank 

if  no  generation  is  desired. 
F|         =  State  ofthei-th  DOF,  either  fixed  or  released. 

Therefore  F(  can  have  ONLY  the  two  following  values;  F;  —  F  is 

for  fixed  or  F,  =  R  is  for  released. 

The  boundary  specification  line  may  be  repeated  as  many  times  as  is  necessary  to  define 
the  boundary  conditions  of  all  the  nodes.  Any  node  that  is  not  defined  will  default  to  the 
fully  fixed  condition. 
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This  section  MUST  end  with  a  blank  line. 


WIND  PARAMETERS 


This  section  defines  the  wind  loading  parameters.  The  syntax  rules  are  as  follows: 


WIND  -  (Section  Designation) 


SPEED=S1    DIRECTION=Dl    GUST=G1,G2    K=K1     S=STD 


SI 
Dl 


Average  wind  speed  (mph). 

Wind  approach  angle.  If  wind  impact  angle  is  perpendicular  to  the 

default  x-z  plane,  Dl=  90°.  See  Fig.  B3. 


Dl  =0' 


Dl  =90' 


Fig.  B-3  Wind  Impact  Angle 


Gl 


G2 


Kl 


STD     - 
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If  istat  is  specified  as  -1  in  the  control  section  then  the  gust 

parameters  must  be  specified.  Otherwise  this  must  be  omitted.  Gl 

is  the  duration  of  the  wind  data  to  generate.  Duration  could  range 

from  3  seconds  to  90  seconds.  There  is  a  default  limit  of  90 

seconds  incorporated  into  the  program. 

Sampling  interval  for  the  gust  wind  generation.  This  must  not  be 

less  than  0. 1  seconds.  The  numerical  filters  used  to  generate  the 

wind  load  tends  to  be  divergent  if  G2  is  less  than  0. 1 .  Note  also 

that  G2  can  not  be  greater  than  Gl . 

roughness  factor.  This  is  depended  on  the  environment 

surrounding  the  structure.  Kl  values  can  be  obtained  from  chapter 

four. 

This  the  standard  deviation  from  the  mean  velocity  desired  in  the 

random  velocity  component.  If  omitted,  a  default  value  is 

computed  based  on  the  spectral  density  distribution  in  chapter  4.  If 

the  default  STD  value  is  not  desired  a  value  of  one  implies 

standard  deviation  is  ignored. 


This  section  MUST  end  with  a  blank  line. 


ELEMENT  DATA 


One  section  each  must  be  supplied  for  the  different  element  types  that  makes  up  the 
structure.  The  units  of  element  properties  unless  otherwise  specified  must  be  in  kips  and 
inches. 
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CABLE  -  (Section  Designation) 

NM,NP 

NM      =  Number  of  cable  members. 

NP       =  Number  of  material  properties. 

Material  Property  Lines  -  NP  lines 

MID  A=A  E=E  D=D 

MID     =  Material  property  set  number. 

A         =  Cross  sectional  area,  A. 

E         =  Modulus  of  elasticity,  E. 

D         =  Damping  ratio,  D. 

Cable  Element  Lines  -  NM  lines 

NN,NI,NJ  M=MP  C=C, 

NN  =  Member  identification  number. 

NI  =  Joint  number  I  (Starting  node). 

NJ  =  Joint  number  j  (Ending  node). 

MP  =  Material  property  number. 

C,  =  Number  of  cable  that  the  element  belongs. 


This  section  MUST  be  terminated  with  a  blank  line. 
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TRUSS  -  (Section  Designation) 

NM,NS 

NM      =  Number  of  truss  members. 

NS       =  Number  of  material  properties. 

Material  property  lines  -  (NS  lines). 

MID  A=A  E=E  D=D  M=M 


MID 

A 
E 
D 
M 


Material  property  set  number. 
Cross  Sectional  Area,  A. 
Modulus  of  elasticity,  E. 
Damping  ratio 
mass  density  of  material,  (lbs/ft3). 


TRUSS  element  lines  -  NM  lines 


NN,NI,NJ  M=MP 


NN       =  Member  identification  number. 

NI        =  Joint  number  I  (Starting  node). 

NJ        =  Joint  number  J  (Ending  node). 

MP      =  Material  property  number. 

This  section  MUST  be  terminated  with  a  blank  line. 
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THREE  DIMENSIONAL  BEAM,  AND  CONNECTOR  ELEMENTS 


There  are  two  types  of  three  dimensional  beam  elements  used  in  DYNASS.  Both  types 
of  beam  elements  are  used  to  model  members  which  have  bending  stiffness  about  the  two 
principal  axes  and  axial  stiffness.  Shearing  deformations  are  neglected  for  both  elements. 
The  members  are  prismatic  and  connected  to  nodes  i  and  j.  A  dummy  node,  k,  is  used  to 
specify  the  direction  of  the  principal  bending  axis.  Axis  1  is  the  i-j  direction;  Axis  3  is 
normal  to  axis  1  and  is  in  the  i-j-k  plane;  Axis  2  is  normal  to  the  plane  3  formed  by  the  1 
and  3  axes  and  completes  the  right  hand  rule.  The  orientation  of  the  axes  is  shown  in  Fig. 
B.3. 


M1 


ef_ 


A 


V3 
M3 


Fig.  B.3  Axis  Orientation 


The  information  for  each  type  must  be  supplied  in  a  separate  block.  The  Beam  element  is 
a  full  3-D  beam  which  serves  as  a  pole  in  the  structure.  The  Connector  element  is  a  large 
displacement  3-D  beam  element  and  it  connects  the  catenary  cable  with  the  messenger 
cable  and  the  signal  light  or  sign. 
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BEAM  -  (Section  Designation) 

NM,NP 

NM      ■  Number  of  beam  members. 

NP       =         Number  of  material  types. 

Material  Property  Lines  -  NP  lines 

MID  I=I2,I3  J=J  A=A  E=E  G=G  T=TYPE  M=M  D=D 

MID  =  Material  property  number. 

A  =  Cross  sectional  area. 

J  =  Torsional  moment  of  inertia  about  the  1  axis. 

12  =  Moment  of  inertia  about  the  2  axis,  I22.  Strong  Axis. 

13  =  Moment  of  inertia  about  the  3  axis,  I33.  Weak  Axis. 
E  =  Modulus  of  elasticity  (Young's  Modulus),  E. 

G         =  Shear  Modulus,  G,  only  used  to  calculate  the  torsional  stiffness. 

TYPE  =  1  (beam  used  as  part  of  the  strain  pole).  2  (other  beam  usage 

excluding  connector.  Necessary  for  other  types  of  structures.) 
D         =  Damping  ratio. 

M         =  Mass  density  of  the  beam  material  (lbs/ft3). 

Beam  element  Connectivity  -  NM  lines 

NN,NI,NJ,NK,NPL  M=MP 
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NN  =  Member  number. 

NI  =  Node  number  I  (Staring  node). 

NJ  =  Node  number  j  (Ending  node). 

NK  =  Node  number  k  (Dummy  node  for  orientation). 

NPL  =  Pole  number  to  which  the  beam  belongs. 

If  T=2,  this  value  must  be  omitted. 

MP  =  Material  property  for  this  member.  (DEFAULT=  previous). 

This  section  MUST  be  terminated  with  a  blank  line. 


CONNECTOR  -  (Section  Designation) 


NM,NP 


NM 
NP 


Number  of  connector  members. 
Number  of  material  types. 


Material  Property  Lines  -  NP  lines 


MID  I=I2,I3  J=J  A=A  E=E  M=M  D=D 


MID  =  Material  property  number. 

A  =  Cross  Sectional  area. 

12  =  Moment  of  inertia  about  the  2  axis,  I22. 

13  =  Moment  of  inertia  about  the  3  axis,  I33. 

E  =  Modulus  of  elasticity  (Young's  Modulus),  E. 
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J 
M 

D 


Torsional  moment  of  inertia  about  the  1  axis. 
Mass  density  of  the  connector  material,  (lbs/ft3). 
Default=490  lbs/ft3. 
Damping  ratio. 


Connector  element  Connectivity  -  NM  lines 


NN,NI,NJ,NK  M=MP 


NN       =  Member  number. 

NI        =  Node  number  I. 

NJ        =  Node  number  J. 

NK      =         Node  number  K. 
MP       =  Material  property  set  for  this  member. 

(DEFAULT^  previous). 
This  section  MUST  be  terminated  with  a  blank  line. 

SIGNAL/SIGN  ELEMENTS 


These  two  elements  are  not  like  the  other  elements  because  their  stiffness  properties  are 
never  calculated.  The  signals  and  sign  are  specified  to  compute  wind  loading.  If  the 
typical  signal  heads  or  auxiliary  signs  are  used  [appendix  C,  ref  1 9]  the  program  is 
equipped  to  compute  the  necessary  areas  and  weights.  If  a  new  atypical  cladding  is  used, 
the  sign  type  OTHER  should  be  selected. 
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SIGNAL  -  (Section  Designation') 

NM,NP 

NM      =  Number  of  signals  and  /or  signs  nodes. 

NP       =  Number  of  different  signal  or  sign  specification. 

Property  Specification  Lines  -  NP  lines 

MID  S=S  D=D„..,DS  E=E„...,ES  L=L„...,LSG=G„...,GS 

or 

MID  N=N  A=A 

MID     =  Material  property  number. 

S  ■  Type  of  bracket.  S=l  indicates  a  one-way  bracket  hence  one  signal 

head  attached  to  the  node.  S=4  is  the  largest  number  of  signal 

heads  that  can  be  attached  to  any  one  node. 

This  indicates  a  4-way  light. 
Dj        =  Diameter  of  the  lens.  There  are  only  two  choices.  8"  or  12"  lens. 

A  choice  of  lens  diameter  should  be  specified  for  each  face  of  light 

depending  on  the  value  of  S  specified.  1=1  ..S 
Ej         =  Border  preference  for  the  signal  light.  The  choices  are  without, 

WO,  with  5"  border  backplate,  W5,  and  with  8"  border  backplate, 

W8.. 
Lj         =  Material  type  for  the  signals.  There  are  only  two  choices, 

Aluminum,  AL,  or  Poly  Carbonate,  PC. 
Gj         =  Geometry  of  the  signal  face.  See  ref  [19]  and  Table  B-4  for  details. 

The  order  the  signal  faces  must  be  specified  is  clock-wise  starting 
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from  the  face  parallel  to  the  x-z  plane.  So  if  Fig.  B-3  is  referred  to. 
The  face  first  specified  corresponds  to  the  wind  impact  angle  Dl  of 
90°,  the  next  face  corresponds  to  an  impact  angle  Dl  of  0°  and  the 
next  to  an  impact  angle 

of  270°  and  finally  the  face  corresponding  to  an  impact  angle  of 
180°.  If  the  structure  exist  in  the  y-z  plane,  i.e.  facing  the  impact 
wind  direction  of  0°,  a  zero  specification  (e.g.  D=  0,D2 )  must 
precede  the  lens  diameter  list  for  a  single  bracket  light. 


Alternatively  an  auxiliary  sign  could  be  specified. 


N 


Typical  Auxiliary  sign  notation.  See  ref  [19]  for  details. 
If  N=OTHER  is  specified.  Then  specify  area. 
Surface  area  of  the  sign  or  atypical  signal.  Note  this  is  only 
specified  in  conjunction  with  the  sign  type  OTHER. 
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Table  B-4  Light  face  arrangement  for  G  specification 


Geometry  Denotation 


Light  Configuration 


Description 


G=l 


1  section  light 


Used  as  blinking  caution  light 


G=2 


2  linear  section 
Commonly  at  railway 
crossing. 


G=3 


3  linear  section 


G=4 


4  linear  section 


8 


G=5 


5  linear  section 


G=6 


5  special  (dog  house 
type) 


. 


o 


G=7 


5T  special 


• 
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Signal/Sign  Node  Assignment  -  NM  lines 

NN,NM   M=MP 

NN      =  Signal  or  Sign  number. 

NM      =  Signal  or  Sign  node. 

MP      -  Signal  Light  or  Traffic  Sign  type  to  assign  to  the  node. 

This  section  MUST  be  terminated  with  a  blank  line. 

LOADS  -  (Section  Designation) 

NFF=FX,FY,FZ 


NF 

Fx 
FY 
Fz 


Node  on  which  the  load  will  act. 

Load  corresponding  to  the  X  transnational  degree  of  freedom. 
Load  corresponding  to  the  Y  transnational  degree  of  freedom. 
Load  corresponding  to  the  Z  transnational  degree  of  freedom. 


The  load  line  specification  line  may  be  repeated  as  many  times  as  is  necessary  to  define 
all  the  loads  acting  on  the  structure.  Any  loads  not  specified  will  default  to  zero. 
This  section  MUST  end  with  a  blank  line. 
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MEMORY  MANAGEMENT  AND  SPECIAL  FEATURES  IN  DYNASS 

It  is  necessary  to  discuss  special  features  in  the  program.  One  of  such  feature  is  tied  to 
the  dynamic  memory  management  employed  in  the  program.  For  a  more  efficient 
memory  management  the  program  uses  the  profile  equation  solver.  The  profile  form  of 
the  mass,  stiffness  and  damping  matrices  are  stored  because  of  the  skyline  form  and  the 
sparseness  of  these  matrices.  The  inverse  of  the  assembled  mass  and  effective  stiffness 
matrices  was  required  at  various  stages  in  the  analysis  process,  matrix  inversion  was 
performed  using  a  typical  inversion  technique.  The  commonly  used  algorithms  used  to 
perform  a  matrix  inversion  include  the  Gauss  elimination,  standard  Cholesky  (A=m) 
and  the  modified  Cholesky  (A=LDLT).  The  modified  Cholesky  is  adopted  in  this 
program.  The  profile  solver  programmed  into  DYNASS  is  a  based  on  Carlos  Fellipa's 
Profile  solver  given  in  ref  [33].  The  Algorithm  for  matrix  inversion  is  simply; 

Ax  =  b;  A  is  then  x  n  matrix,  x  is  the  unknown  displacement 

vector,  b  is  the  force  vector. 

A    =  LDU;      Factorize  A.  Lis  the  unit  lower  triangular  matrix, 
D  is  a  non-singular  diagonal  matrix, 
U  is  the  transpose  ofL. 

Lz  =  b;  Forward  substitution. 

y   =z/D;       Scaling. 

Ux=y;  Back  substitution. 
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The  above  procedure  is  typical  of  all  modified  Cholesky  solution  procedure.  What  is  not 
typical  is  the  storage  of  the  coefficients  of  matrix  [A].  If  [A]  is  defined  below  a  9  x  9 
matrix  with  zero  terms. 


A  = 


a 


11 


a 


a 


12 

22 


SYMM. 


a 


23 


a 


a. 


33       34 


a 


44 


a 


a 


a 


a 


15 
0 

a28 

35 

a38 

45 

"46 

a47 

*48 

55 

"56 

0 

fl58 

fl59 

"66 

%1 

0 

"69 

all 

"78 

"79 
0 

88 


a 


99 


Matrix  A  has  a  skyline  profile.  The  profile  matrix  stored  is  as  shown  below. 

H"ll     a12    "22     "23    "33    "34    "44    "l5     °    "35 
"45    "55    "46    "56    "66    "47     °    "67    "77    "28 


"38    "48    "58     ° 


"78    "88    "59     "69    "79     °      "99} 
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A  pointer  array  was  developed  to  accurately  store  and  recover  the  above  profile  matrix. 

The  pointer  array,  referred  to  as  the  NA  array  in  the  program,  for  the  above  matrix  is 

therefore 

NA  =  {\    3    5    7    12    15    19    26    31} 


This  locates  the  terms  along  the  diagonal  of  the  original  [A]  matrix.  To  conform  with 
Carlos  Fellipa's  profile  solver  [33].  There  is  a  negative  unit  regressive  shift  of  the  above 
pointers,  such  that  the  pointer  now  pointing  to  the  first  diagonal  term  in  the  [A]  matrix  is 
zero.  The  equivalent  Fellipa's  NA  is 
NA  =  {0    2    4    6    11    14    18    25    30} 

One  of  the  limitations  in  the  program  has  to  do  with  the  maximum  size  of  the  NA  array. 
The  default  maximum  permissible  size  is  2000.  The  above  example  has  an  NA  size  of  9, 
corresponding  to  the  number  of  degrees  of  freedom.  Another  memory  restriction  is  the 
size  of  the  profile  matrix.  For  the  above  example  the  profile  matrix  size  is  3 1 .  For  the 
DYNASS  program  the  default  maximum  profile  size  is  6000.  Any  attempt  to  exceed 
these  limits  results  in  a  memory  restriction  error  to  be  printed  and  an  abrupt  termination 
of  the  program.  Included  is  a  sample  INPUT  and  OUTPUT  files. 
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SAMPLE  INPUT  AND  OUTPUT  FILES  FOR  DYNASS 

The  sample  input  file  included  is  for  the  wind  load  test  #1  discussed  in  chapter  6.  All 

coordinates,  member  properties  and  connectivity  information  related  to  Fig.  6. 1 .  The 

modified  Newmark's  Beta  method  is  adopted  for  analysis  with  a  tolerance  of  0.000133 

and  numerical  parameters  of -0.1  for  alpha,  0.3025  for  beta  and  0.6  for  gamma,  he 

average  wind  speed  is  39.941  and  wind  model  2  (i.e.  istat=  -1)  is  used.  The  wind  is 

generated  for  2.5  sec,  same  as  for  total  time  of  analysis.  Time  step  of  analysis  is  0.0003 

but  time  step  for  wind  speed  generation  is  0.1.  Dynamic  response  for  node  16,  y 

direction  is  printed  as  output  plot  file  SVEL. 

First  Wind  Test  --  Low  Wind  Speed  load  test 

nodes=18,  elem_no=3,  cable=2,  istat=-l 

nonlinear 

6  nd=3, 5, 6, 7, 8, 9  sg=5.00  ten=0.6  w=0.27  p=0 

6  nd=2,ll,12,13,14,15  ten=0.829  w=0.27  p=l 

dynamics 

t=0. 000133  a=-0.1  b=0.3025  g=0.6  time=2 . 5, 0 . 0003  p=16  n=2 

wind 

speed=39.941  direction=90  gust=2 . 5  ,0.1  k=0.05  s=l 


coordinates 

1 

x= 

0 

y= 

0 

.000 

z  = 

16 

2 

x= 

0 

y= 

0 

.000 

z  = 

112 

3 

x= 

0 

y= 

0 

.000 

z  = 

164 

4 

x= 

0 

y= 

0 

.000 

z  = 

166 

5 

x= 

15*12 

y= 

0 

.000 

z  = 

145. 

28 

6 

x= 

30*12 

y= 

0 

.000 

z  = 

126. 

55 

7 

x= 

42*12 

y= 

0 

,000 

z= 

126. 

55 

8 

x= 

57*12 

y= 

0. 

,000 

z  = 

145. 

28 

9 

x= 

72*12 

y= 

0. 

000 

z  = 

164 

10 

x= 

72*12 

y= 

0. 

000 

z  = 

142 

11 

x= 

15*12 

y= 

0. 

000 

z  = 

112 

12 

x= 

30*12 

y= 

0. 

000 

z  = 

112 
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13  x=  42*12  y=  0.000   z  =  112 

14  X=  57*12  y=  0.000   Z  =  112 

15  X=  72*12  y=  0.000   z  =  112 

16  x=  30*12  y=  0.000   z=  112-27.3 

17  x=  42*12  y=  0.000   z=  112-27.3 

18  x=  72*12  y=  0.000   z=  0 

boundary 

1,18,1  dof=r  r  r  r  r  r 

I  dof=f  f  f  f  f  f 
5  dof=r  r  r  f  f  f 
8  dof=r  r  r   f   f   f 

II  dof=r  r  r  f  f  f 
14  dof=r  r  r  f  f  f 
18  dof=f  f  f  f  f  f 


Cable 

10,1 

1   A=0.079  E=23000  D=0 . 05 


1   3,5   M=l  C=l 

2   5,6   M=l  C=l 

3   6,7   M=l  C=l 

4   7,8   M=l  C=l 

5   8,9   M=l  C=l 

6   2,11  M=l  C=2 

7   11,12  M=l  C=2 

8   12,13  M=l  C=2 

9   13,14  M=l  C=2 

10   14,15  M=l  C=2 

beam 

6,2 

1  A=93.4  1=888,888 

J=484  E= 

4000 

M=195.0  T=l 

D=0.05 

:EI  =  3  550  000 

2  A=93.4  1=820,820 

J=484  E= 

4000 

M=195.0  T=l 

D=0.05 

:EI  =  3  280  000 

1    1,2,18,1   M=l 

2    2,3,18,1   M=l 

3    3,4,18,1   M=l 

4    9,10,1,2   M=2 

5    10,15,1,2   M=2 

6    15,18,1,2   M=2 
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connector 
4,2 

1  A=0.799  1=0.310,0.310  J=0.620  E=29000  M=490.2  D=0.05  :  pipe  connector 

2  A=0.5324  1=0.167,0.01534  J=0.50  E=29000  M=490.2  D=0.05  :  strap 

1  6,12,18   M=l  P=0. 07140 

2  12,16,18   M=l  P=0. 07140 

3  7,13,18   M=2  P=0. 07140 

4  13,17,18   M=2  P=0. 07140 

signal 

2,1 

1  R=1.5  S  =  l  D=12  E=WO  L=AL  G=6 

1  16  m=l 

2  17  m=l 

loads 

12  F=0,0, -8.15/1000,0,0,0 

13  F=0, 0,-8. 15/1000, 0,0,0 

6  F=0,0, -4.35/1000,0,0,0 

7  F=0, 0,-4. 35/1000, 0,0,0 


The  corresponding  OUTPUT  file  for  the  previous  INPUT  file  is  given  next.  Note  the 
units  used  as  output  is  kips  and  inches  unless  otherwise  stated.  The  OUTPUT  includes  an 
echo  of  the  input  parameters  and  the  results  of  computations.  There  also  two  different 
output  generated  by  the  program,  a  displacement  plot  in  the  file  PLOT  and  the  wind 
generation  information  given  in  the  plot  file  SVEL.    The  header  for  the  PLOT  file  is 
time,  nodal  displacement,  nodal  velocity  and  nodal  acceleration.  The  header  for  the 
SVEL  file  is  time,  wind  speed,  nodal  wind  drag  force,  tension  in  the  catenary  cable(s)  and 
finally  the  tension  in  the  messenger  cable(s). 
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+ 

##### 

+ 

#       # 

+ 

#       # 

+ 

#       # 

+ 

#       # 

+ 

##### 

+ 

+ 

++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

#     #    #       #       ##        ####       ####  + 

########  #  + 

#  #  #    #    #       #     ####       ####  + 

#  #  #  #  ######     #     #  + 

#  ####       ##       ##       # 

#  #       #    #       #     ####       #### 

* * 

Developed 
By 

Dr.  Adeola  K.  Adediran 

In  Conjuction 
With 

Dr.  Marc  Hoit 
Dr.  Ronald  Cook 
Dr.  Petros  Christou 

* * 


+ 
+ 
+ 


+ 
+ 
+ 
+ 
+ 
+ 


+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 


++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


Title  and  General  Information 


First  Wind  Test  --  Low  Wind  Speed  load  test 
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Number  of  joints 
Number  of  element  types 
Number  of  cables 


18 
3 
2 


Units  are  in  Kips  and  Inches  unless  otherwise  stated 


ECHO  OF  INPUT  DATA 


A.  Boundary  conditions 


B.  Input  nodal  point  coordinates 


Node  X 


XX   YY 


zz 


X 


1 

F 

F 

F 

F 

F 

F 

0.000 

0.000 

16.000 

2 

R 

R 

R 

R 

R 

R 

0.000 

0.000 

112.000 

3 

R 

R 

R 

R 

R 

R 

0.000 

0.000 

164.000 

4 

R 

R 

R 

R 

R 

R 

0.000 

0.000 

166.000 

5 

R 

R 

R 

F 

F 

F 

180.000 

0.000 

145.280 

6 

R 

R 

R 

R 

R 

R 

360.000 

0.000 

126.550 

7 

R 

R 

R 

R 

R 

R 

504.000 

0.000 

126.550 

8 

R 

R 

R 

F 

F 

F 

684.000 

0.000 

145.280 

9 

R 

R 

R 

R 

R 

R 

864.000 

0.000 

164.000 

10 

R 

R 

R 

R 

R 

R 

864.000 

0.000 

142.000 

11 

R 

R 

R 

F 

F 

F 

180.000 

0.000 

112.000 

12 

R 

R 

R 

R 

R 

R 

360.000 

0.000 

112.000 

13 

R 

R 

R 

R 

R 

R 

504.000 

0.000 

112.000 

14 

R 

R 

R 

F 

F 

F 

684.000 

0.000 

112.000 

IB 

R 

R 

R 

R 

R 

R 

864.000 

0.000 

112.000 

16 

R 

R 

R 

R 

R 

R 

360.000 

0.000 

84.700 

17 

R 

R 

R 

R 

R 

R 

504.000 

0.000 

84.700 
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18     F     F     F     F     F     F 


864.000 


0.000 


0.000 


C.  Element  information 


Beam  Element  information 

Property  Number       =    1 

Modulus  of  Elasticity  =    4000.00 

Area  =   93.400000 

Moment  of  Inertia  (Strong  axis)=  888.000000 

Moment  of  Inertia  (Weak   axis)=  888.000000 

Torsional  Moment  of  Inertia    =   484.0000 

Shear  Modulus  =    173  9.13 

Damping  Ratio  =   0.050000 

Material  density  =  195.000000   lbs/ft*3 

Property  Number       =   2 

Modulus  of  Elasticity  =    4000.00 

Area  =   93.400000 

Moment  of  Inertia  (Strong  axis)=  820.000000 

Moment  of  Inertia  (Weak   axis)=  820.000000 

Torsional  Moment  of  Inertia    =   484.0000 

Shear  Modulus  =    1739.13 

Damping  Ratio  =    0.050000 

Material  density  =  195.000000   lbs/ft*3 

-  Number  of  Elements  =   6 


Member  #   Beam  Type   Node  I   Node  J   Node  K   Material  # 
1  1        1        2       18        1 
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2 

3 
4 
5 
6 


1 

2 

3 

18 

1 

1 

3 

4 

18 

1 

1 

9 

10 

1 

2 

1 

10 

15 

1 

2 

1 

15 

18 

1 

2 

Beam  type  1  -  Pole 

Beam  type  2  -  Other  Beam  Usage 


Cable  element  information 

Property  Number       =  1 

Modulus  of  Elasticity  =  23000.00 

Area  =  0.079000 

Damping  Ratio         =  0.050000 


Number  of  Elements 


10 


Member  # 

Node  I 

Node  J 

Ma 

terial  # 

Cable 

# 

Initial  Tension 

(kips) 

l 

3 

5 

1 

1 

0.652482 

2 

5 

6 

1 

1 

0.652018 

3 

6 

7 

1 

l 

0.648627 

4 

7 

8 

1 

l 

0.652018 

5 

8 

9 

1 

1 

0.652482 

6 

2 

11 

1 

2 

0.829000 

7 

11 

12 

1 

2 

0.829000 

8 

12 

13 

1 

2 

0.829000 

9 

13 

14 

1 

2 

0.829000 

10 

14 

15 

1 

2 

0.829000 

Cable  initial  sag 

Cable   1:   0.050000 
Cable   2:   0.000000 


Connector  Element  information 
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Property  Number       =   1 

Modulus  of  Elasticity 

Area  = 

Moment  of  Inertia  (Strong  axis)= 

Moment  of  Inertia  (Weak   axis)= 

Torsional  Moment  of  Inertia    = 

Shear  Modulus  = 

Damping  Ratio 

Material  density  = 


29000.00 
0.799000 
0.310000 
0.310000 
0.62 
11153.85 
0.050000 
490.200000  lb/ftA3 


Property  Number       =   2 

Modulus  of  Elasticity  = 

Area 

Moment  of  Inertia  (Strong  axis)= 

Moment  of  Inertia  (Weak   axis) = 

Torsional  Moment  of  Inertia 

Shear  Modulus 

Damping  Ratio 

Material  density 


29000.00 
0.532400 
0.167000 
0.015340 
0.50 
11153.85 
0.050000 
490.200000  lb/ft*3 


Number  of  Elements  = 


Member  #   Node  I   Node  J   Node  K   Material  # 
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i 

2 
3 
4 
M=  1  Single  Light 


6 

12 

ie 

12 

16 

18 

7 

13 

18 

13 

17 

18 

D=12 

E  = 

WO 

ii- 

G=** 

1 
1 

2 

2 


Signal  and  Sign  information 

Sign(al)  #   Node  #   Material  #   Weight  (kips)    Area  (inA2;  1-3,  1-2 
planes) 

1  16       1  0.71E-01        0.11E+04        0.00E+00 

2  17       1  0.71E-01        0.11E+04        0.00E+00 


***  Load  Analysis  Started  *** 


G.  Concentrated  Nodal  Loads 

Node  Load         x         y         z        xx        yy  zz 

12  1   0.00E+00   0.00E+00  -0.82E-02  0.00E+00  0.00E+00  0.00E+00 

13  1   0.00E+00   0.00E+00  -0.82E-02  0.00E+00  0.00E+00  0.00E+00 

6  1   0.00E+00   0.00E+00  -0.43E-02  0.00E+00  0.00E+00  0.00E+00 

7  1   0.00E+00   0.00E+00  -0.43E-02  0.00E+00  0.00E+00  0.00E+00 


***  Formation  of  Stiffness  Matrix  *** 


E.  Time  Step  Analysis: 
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Time  Step  Interval   =     0.000300  sees. 
Max.  Time  Allowed   =    90.000000  sees, 


F.  Displacements  After  Die-Down 


Node 
l 
2 
3 
4 
5 
6- 
7- 
8- 
9- 
10- 
11 
12- 
13- 
14- 
15- 
16- 
17 
18 


0.00000E+00 
0.16734E+00 
0.32096E+00 
0.32701E+00 
0.38091E-01 
0.84318E-01 
0.92674E-01 
0.21992E+00 
0.49211E+00 
0.40003E+00 
0.52082E-01 
0.60736E-01 
0.48410E-01 
0.16262E+00 
0.27777E+00 
0.16572E-01 
0.30688E-01 
0.00000E+00 


0.00000E+00 
0.43437E-02 
0.81251E-02 
0.82726E-02 
0.25155E+01 
0.50543E+01 
0.51359E+01 
0.26392E+01 
0.12680E-01 
0.10382E-01 
0.67560E+01 
0.13502E+02 
0.13499E+02 
0.67573E+01 
0.72958E-02 
0.29365E+02 
0.29437E+02 
0.00000E+00 


0.00000E+00  0 

•0.18611E-04-0 

0.28003E-04-0 

■0.28003E-04-0 

■0.22159E+01    0 

■0.30377E+01    0 

•0.29750E+01    0 

•0.21046E+01    0 

■0.30868E-04-0 

■0.26925E-04-0 

■0.52327E+00    0 

■0.33418E+00    0 

■0.33101E+00    0 

■0.50794E+00    0 

■0.21548E-04-0 

0.47473E+01  0 

0.48043E+01  0 

0.00000E+00  0 


XX 

00000E+00 

70688E-04 

73741E-04 

73741E-04 

00000E+00 

61937E+00 

61044E+00 

00000E+00 

10466E-03 

.10406E-03 

00000E+00 

,  61966E+00 

.  61621E+00 

.  00000E+00 

10130E-03 

62021E+00 

.  62699E+00 

00000E+00 


YY 

0.00000E+00  0 

0.28118E-02  0 

0.30254E-02  0 

0.30254E-02  0 

O.OOO0OE+00  0 

•0.13009E-02-0 

•0.37933E-02    0 

0.00000E+00  0 

-0.41988E-02  0 

•0.41575E-02    0 

0.00000E+00  0 

•0.12998E-02-0 

-0.37534E-02  0 

0.00000E+00  0 

•0.39682E-02    0 

•0.12980E-02-0 

-0.36755E-02    0 

0.00000E+00  0 


zz 

.  00000E+00 
.00000E+00 
.OOOOOE+00 
. 00000E+00 
.00000E+00 
.96823E-03 
.72197E-04 
. 00000E+00 
. 00000E+00 
.00000E+00 
. O00OOE+O0 
.96690E-03 
.12170E-03 
. 00000E+00 
.  00000E+00 
. 96452E-03 
.21618E-03 
. 00000E+00 


H.  Final  Displaced  Configuration 


Node 

1 


X 
0.000 


Y 
0.000 


z 

16.000 
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2 

0. 

167 

0. 

004 

112. 

000 

3 

0. 

321 

0. 

008 

164. 

000 

4 

0. 

327 

0, 

008 

166. 

,000 

5 

180. 

038 

2. 

,516 

143. 

,064 

6 

359. 

,916 

5. 

.054 

123, 

,512 

7 

503. 

.907 

5 

.136 

123 

.575 

8 

683, 

.780 

2 

.639 

143 

.175 

9 

863 

.508 

0 

.013 

164 

.000 

10 

863 

.600 

0 

.010 

142 

.000 

11 

180 

.052 

6 

.756 

111 

.477 

12 

359 

.939 

13 

.502 

111 

.666 

13 

503 

.952 

13 

.499 

111 

.669 

14 

683 

.837 

6 

.757 

111 

.492 

15 

863 

.722 

0 

.007 

112 

.000 

16 

359 

.983 

29 

.365 

89 

.447 

17 

504 

.031 

29 

.437 

89 

.504 

18 

864 

.000 

0 

.000 

0 

.000 

***  Force  Recovery  Started  *** 


I .  Element  Forces 


1.  Primary (Catenary)  cables 


Cable  Member  Forces 


Cable      Member       Force      Stress 


1  0.5648E+00  0.7149E+01 

2  0.5643E+00  0.7143E+01 

3  0.5608E+00  0.7098E+01 

4  0.5637E+00  0.7136E+01 

5  0.5642E+00  0.7141E+01 


2.  Secondary (Messenger)  cables 


Cable  Member  Forces 
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Cable 


Member 


Force 


Stress 


6  0.9992E+00  0.1265E+02 

7  0.9991E+00  0.1265E+02 

8  0.9986E+00  0.1264E+02 

9  0.9997E+00  0.1265E+02 
10  0.9997E+00  0.1265E+02 


Beam  Member  Forces 


1-2  Plane  1-3  Plane 

Member  Node       Moment       Shear      Moment       Shear  Axial  Force 

1     1    0.1789E+03   0.4580E-01   0.4814E+01   -.1560E+01  0.7243E-01 

2    -.2918E+02   -.4580E-01   -.4170E+00   0.1560E+01  -.7243E-01 

Axial  Torque  =   0.00000E+00 


2  0.2918E+02   0.7972E-02   0.4158E+00   - . 5611E+00      0.6748E-01 

3  0.4148E-03   -.7972E-02   -.1288E-02   0.5611E+00      -.6748E-01 

Axial  Torque  =   0.00000E+00 
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3  0.3470E-05   0.6444E-05   0.1571E-04   -.1175E-04      0.4226E-07 

4  0.2003E-04   -.6444E-05   -.2825E-05   0.1175E-04      -.4226E-07 

Axial  Torque  =   0.00000E+00 


9 

10 


1937E-02   -.8211E-02   -.2193E-03   0.5600E+00      0.6696E-01 
1232E+02   0.8211E-02   -.1804E+00   -.5600E+00      -.6696E-01 

Axial  Torque  =   0.00000E+00 


10    0.1232E+02   -.8089E-02   0.1803E+00   0.5590E+00  0.6696E-01 

15    -.2909E+02   0.8089E-02   -.4230E+00   -.5590E+00  -.6696E-01 

Axial  Torque  =   0.00000E+00 

15    0.2906E+02   -.4547E-01   0.4205E+00   0 . 1556E+01  0.7188E-01 

18    -.2034E+03   0.4547E-01   -.5513E+01   -.1556E+01  -.7188E-01 

Axial  Torque  =   0.00000E+00 


Connector  (Hanger)  Forces 


1-2  Plane  1-3  Plane 

Member  Node       Moment       Shear      Moment       Shear  Axial  Force 

1     6    -.7777E-05   0.2472E-01   0.2399E-05   -.1774E-03  -.4768E-01 

12    0.2588E-02   -.2472E-01   0.3597E+00   0.1774E-03  0.4768E-01 

Axial  Torque  =  0.15738E-06 


12    -.2532E-02   -.1317E-01   -.3597E+00   0.9571E-04      -.5883E-01 
16    -.8104E-04   0.1317E-01   0.1136E-03   -.9571E-04      0.5883E-01 

Axial  Torque  =   0.14089E-06 
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3  7    0.1534E-03   0.2419E-01   -.2255E-03   -.5635E-04  -.4790E-01 
13    0.6665E-03   -.2419E-01   0.3522E+00   0.5635E-04  0.4790E-01 

Axial  Torque  =  0.19264E-05 

4  13    -.2507E-03   -.1291E-01   -.3519E+00   0.8344E-05  -.5920E-01 
17    0.2288E-04   0.1291E-01   -.5080E-03   -.8344E-05  0.5920E-01 

Axial  Torque  =  0.40491E-05 
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